The Dark Universe 



On 
O 
O 
(N 

C 
>—> 

(N 

o 

u 

43 

Oh 
i 

o 



> 

m 
o 
in 

*sD 
O 
OS 
O 

> 

X 

C3 



Matthias BartelmanrQ 

Zentrum fur Astronomie der Universitat Heidelberg, Institut fur Theoretische Astrophysik, Albert-Uberle-Str. 2, 
69120 Heidelberg, Germany 

(Dated: draft version of June 27, 2009) 

For a few years now, cosmology has a standard model. By this term, we mean a consistent theoretical background 
which is at the same time simple and broad enough to offer coherent explanations for the vast majority of cos- 
mological phenomena. This review will briefly summarise the cosmological model, then proceed to discuss what 
we know from observations about the evolution of the Universe and its contents, and what we conclude about the 
origin and the future of the Universe and the structures it contains. 



I. THE COSMOLOGICAL STANDARD MODEL 
A. Historical outline 

This review will not follow the historical sequence of dis- 
coveries but rather present the arguments in a perhaps more 
logical order than history has allowed. Yet, let us begin with a 
brief and necessarily fragmental historical outline. 

Modern cosmology began in 1915 with Einstein's theory 
of General Relativity, from which Friedmann in 1922 derived 
the homogeneous and isotropic class of world models which 
to this date form the foundation of the cosmological standard 
model. His discovery that static world models were impossi- 
ble confirmed Einstein in his introduction of the cosmological 
constant, which he abandoned in the 1930s after Slipher, Hub- 
ble and Humason had discovered the cosmic expansion. Ear- 
lier, Lemaitre had speculated about a primordial fireball which 
seemed the natural starting point of an expanding universe. 

First evidence for dark matter was found by Zwicky in 
galaxy clusters, and by Babcock and Oort in galaxies, between 
1933 and 1940. In 1946, Lifshitz worked out relativistic per- 
turbation theory and started applying it to the linear growth 
of cosmic structures. In the late 1940s, Gamow, Alpher and 
Herman worked out how nuclear fusion may have proceeded 
in the early universe and predicted a cosmic radiation back- 
ground with a temperature of a few Kelvin. 

The cosmic microwave background (CMB) was discovered 
by Penzias and Wilson and explained by Dicke and collab- 
orators in 1965. In 1970, Peebles and Yu as well as Sun- 
yaev and Zel'dovich independently predicted structures in the 
CMB and found that they should have relative amplitudes near 
10~ 4 , where subsequent searches did not find them. Peebles 
suggested in 1982 that a lower CMB fluctuation level could 
be explained if the dark matter could not interact with light. 
It was quickly realised by means of the numerical simula- 
tions by Davis, Efstathiou, Frenk and White that cosmolog- 
ical structure formation could be explained within the emerg- 
ing paradigm of cold dark matter. This paradigm experienced 
further strong support when temperature fluctuations at the re- 
vised level were finally found in the CMB with the COBE 



satellite in 1992. 

It was suggested by Guth in 1981 that a phase of inflation- 
ary expansion could solve the horizon and flatness problems 
of Friedmann cosmology. Mukhanov and Chibisov immedi- 
ately saw that inflated quantum fluctuations could be the ori- 
gin of cosmic structures. At least in simple scenarios, inflation 
requires a spatially flat universe and thus a total energy den- 
sity equal to its critical value, but it became obvious that mat- 
ter alone cannot be as dense. The difference was tentatively 
attributed to the cosmological constant, which was confirmed 
by several groups in the late 1990s through the accelerated 
cosmic expansion inferred from type-la supernovae. The cos- 
mological standard model could be considered complete when 
the recent measurements of the CMB temperature fluctuations 
confirmed that the universe is spatially flat. 

Modern cosmology is thoroughly described in many text- 
books. To mention only a few, se e (|Coles and Lucchin||2002| 



Dodelson 



2003; 



Padmanabhanl |1993[ |Peacock[ |1999| |PeeH 



bles] [l9"9"3~t |Weinberg| [2008] l and the textbook on the CMB 
( |Durrer||2008l ~ 



This review is structured as follows. The introduction pro- 
ceeds by summarising the framework of the Friedmann mod- 
els in LB and the growth of cosmic structures in I.C Section|H| 
describes what we know about the age of the Universe from 



its ingredients through nuclear cosmo-chronology (II. A i, stars 
( II. B i and the cooling of white dwarfs ( II. C I. In Sect. Ill mea- 



surements of the Hubble constant from Hubble's law ( III.A[ ) 
and gravitational lensing (III.Bi are reviewed. Big-bang nu- 



cleosynthesis is the subject of Sect. IV with subsections on 
the production of light elements (IV.Ai and the observed el- 
ement abundances (TV.B) . Section~|V| gives an overview of 
the matter density in the Universe, as constrained by galax- 



ies (V.Ai, galaxy clusters as individual objects (V.Bi and as 



an evolving population (V.Cl. Section VI on the cosmic mi- 
crowave background begins with the isotropic CMB ( |VLA| ) 
and continues with its fluctuations ( |VI.B[ >. Cosmic structures 
are reviewed in Sect. VII which starts with the quantifica- 
tion of structures ( |VII.A i and continues with their measure- 
ment and results obtained ( |VII.B) . Section |VIII| on cosmo- 
logical weak lensing follows, in which gravitational light de- 
flection is first summarised ( |VIII.A| i before measurements 
are describe d (|VIII.B I. Supernovae of type la are the sub- 
ject of Sect IX) where the principles of their cosmological 
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application (IX.Ai, the reasons thereof and the observational 
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results are described ( |IX.B| Section [X] on the normalisation 
of the matter-fluctuation power spectrum begins with an in- 



troduction (X.Ai and continues describing constraints from 



CMB fluctuations (X.Bi, cosmological weak lensing (X.Ci 



galaxy clusters (X.D i and the Lyman-ff forest (X.Ei. The final 



Sect.|Xl]discusses the motivation and the evidence for cosmo 



logical inflation (XI. A i and ends with remarks on dark energy 
(jXTB). 



space-time. The scale factor is commonly normalised to 
Aq = 1 at the present time. As usual, the line element ds 
gives the proper time measured by an observer moving by 
{dw, ffc{w)d9, /k{w) sin 8d<p) within the coordinate time inter- 
val dt . For light, ds = 0. 

Observers attached to the coordinate grid (w, 0, (p) on the 
spatial hypersurfaces for constant cosmological time t are 
called comoving. The cosmic expansion does not change the 
comoving coordinates x, but the physical coordinates r = ax. 



B. Friedmann models 

1 . The metric 

Cosmology studies the physical properties of the Universe 
as a whole. The only of the four known interactions which 
can play a role on cosmic length scales is gravity. Electro- 
magnetism, the only other interaction with infinite range, has 
sources of opposite charge which tend to shield each other on 
comparatively very small scales. Cosmic magnetic fields can 
perhaps reach coherence lengths on the order of > 10 Mpc, but 
they are far too weak for them to be important for the cosmic 
evolution. The weak and strong interactions, of course, have 
microscopic range and must thus be unimportant for cosmol- 
ogy as a whole. 

The best current theory of gravity is Einstein's theory 
of general relativity, which relates the geometry of a four- 
dimensional space-time manifold to its material and energy 
content. Cosmological models must thus be constructed as 
solutions of Einstein's field equations. Symmetry assump- 
tions greatly simplify this process. Guided by observations 
to be specified later, we assume that the Universe appears ap- 
proximately identical in all directions of observation; in other 
words, it is assumed to be isotropic on average. While this 
assumption is obviously incorrect in our cosmological neigh- 
bourhood, it holds with increasing precision if observations 
are averaged on increasingly large scales. The assumption 
of isotropy can only be valid in a preferred reference frame 
which is at rest with respect to the mean cosmic motion. The 
motion of the Earth within this rest frame must be subtracted 
before any observation can be expected to appear isotropic. 

The second assumption asserts that the Universe should 
appear equally isotropic about any of its points. Then, it 
is homogeneous. Isotropic and homogeneous solutions for 
Einstein's field equations admit the Robertson-Walker metric 
whose line element is 

ds 2 = -c 2 dt 2 + a 2 (t) [dw 2 + f 2 (w) (d6> 2 + sin 2 6»d0 2 )] , (1) 

with the radial coordinate w and the radial function 



K- ]/2 sm(K l/2 w) (K>0) 
w (K = 0) 

\K\- 1/2 sinh(\K\ 1/2 w) (K<0) 



(2) 



The curvature parameter K, which can be positive, nega- 
tive or zero, has the dimension of an inverse squared length. 
The scale factor a(t) isotropically stretches or shrinks the 
three-dimensional spatial sections of the four-dimensional 



2. Redshift and expansion 

The changing scale of the Universe causes the cosmologi- 
cal redshift z. The wavelength of light from a distant comov- 
ing source seen by a comoving observer changes by the same 
amount as the Universe changes its scale while the light is 
travelling. Thus, the emitted and observed wavelengths A and 
/to, respectively, are related by 



A 



a 



(3) 



where a is the scale factor at the time of emission. The relative 
wavelength change is the redshift, 



z = 



Aq-A 



1 



thus 



l+z 



(l+zr 



(4) 



(5) 



When the metric is inserted into Einstein's field equations, 
two ordinary differential equations result, 



c-f - 



8ttG 
4nG 



Kc 1 



i a 



P + 



* + 3 ' 
3p\ A 



(6) 
(7) 



relating a and a to the density p, the pressure p, the cosmo- 
logical constant A and the curvature K. Equation (|7]i can be 
eliminated and replaced by 



— (pc a ) + p— - — 
dt dt 



= 0, 



(8) 



from which the change of p with a can be inferred once p is 
given. Relativistic matter ("r", "radiation") has p = pc 2 /3 
while nonrelativistic matter ("m", "dust") is characterised by 
pc 2 » p « 0. Thus 



p r oc a 



-3 



(9) 



For a model universe containing only matter, radiation and 
the cosmological constant, Eq. ([6]) can be brought into the 
form 



H 1 



Hi 



a 4 a? 



1 - Q 



mO ' 



AO 



= H 2 Q E 2 (a) 



(10) 
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This is commonly called Friedmann's equation, in which the 
relative expansion rate a/a is replaced by the Hubble func- 
tion H(a) whose present value is the Hubble constant Hq = 
H(ao) = H(l), and the matter-energy content is described by 
the three density parameters Q^), QmO and Qao. 
Hq is often expressed in dimension-less form by 



h = 



H 



100 km s-' Mpc- 



(11) 



Since lengths in the Universe are typically measured with re- 
spect to the Hubble length, they are often given in units of 
/z^'Mpc. Similarly, luminosities are typically obtained by 
multiplying fluxes with squared luminosity distances and are 
thus often given in units of h~ 2 L B . We avoid this notation 
where possible in the following and insert h = 0.7 where ap- 
propriate. 

The dimension-less parameters O m o and Q.^ describe the 
densities of matter and radiation today in units of the critical 
density 



_ 3H 2 _ 3flg 

Pa ' ~ 8;rG ' Pcr ° " SttG 



(12) 



A Robertson-Walker metric whose scale factor satis- 
fies Friedmann's equation is called a Friedmann-Lemaitre- 
Robertson- Walker metric. The cosmological standard model 
asserts that the Universe at large is described by such a metric, 
and is thus characterised by the four parameters Q m o, Q-o, Q A o 
and //(). Since the critical density and the densities themselves 
evolve in time, so do the density parameters. 

The remainder of this review article is devoted to answering 
two essential questions: (1) What are the values of the param- 
eters characterising Friedmann's equation? (2) How can we 
understand the deviations of the real universe from a purely 
homogeneous and isotropic space-time? 

Table [I] adapted from (Komatsu et al. 20081, summarises 
the most important cosmological parameters and the values 
adopted below, unless stated otherwise. 



3. The radiation-dominated phase 

It is an empirical fact that the Universe is expanding. Earlier 
in time, therefore, the scale factor must have been smaller than 
today, a < 1. In principle, it is possible for Friedmann models 
to never reach a vanishing scale, a — 0, within finite time 
into the past. However, a few crucial observational results 
suffice to rule out such "bouncing" models. This implies that a 
Universe like ours which is expanding today must have started 
from a — a finite time ago; thus, there must have been a Big 
Bang. 

Equation ( [T0| shows that the radiation density increases like 
gT 4 as the scale factor decreases, while the matter density in- 
creases with one power of a less. Even though the radiation 
density is very much smaller today than the matter density, 
there must have been a period in the early evolution of the 
Universe when radiation dominated the energy density. The 
scale factor when both densities were equal is called a eq . This 



radiation-dominated era is very important for several obser- 
vational aspects of the cosmological standard model. Since 
the radiation retains the Planckian spectrum which it acquired 
in the very early Universe in the intense interactions with 
charged particles, its energy density is fully characterised by 
its temperature T. As the energy density is both proportional 
to T 4 and cT 4 , its temperature falls like T oc a -1 , 

At the times relevant for our purposes, only photons and 
neutrinos need to be considered as relativistic species. Since 
photons are heated by electron-proton annihilation after neu- 
trinos decoupled from the cosmic fluid, the photons are heated 
above the neutrinos by a factor (H/4) 1 ^ 3 , which follows di- 
rectly from entropy conservation. If T is the temperature of 
the cosmic microwave background after electron-positron an- 
nihilation, the energy densities in the photons and the three 
neutrino species are 



Pr.CMB 



n 2 (kT) 4 
15 (he? 



4/3 

:|— | Pr.CMB (13) 



or, as long as all neutrino species are relativistic, 



: Pr.CMB 



21 / 4 



8 111 



4/3 



1.68p r ,cMB ■ (14) 



4. Age, distances and horizons 

Since Friedmann's equation gives the relative expansion 
rate a I a, we can use it to infer the cosmic time, 



r' _ r w da_ r w da 1 r w 

Jo Jo a Jo aH(a) H J 



da 



aE(a) 
(15) 

which illustrates that the age scale is the inverse Hubble con- 
stant H Q 1 . A simple example is given by the Einstein-de Sitter 
model, which (unrealistically) assumes O m o = 1, Q. T o = and 
Q A = 0. Then, from Eq. ([ToJ, E(a) = a~ 3/2 and 

1 P 1 2 
1 = TT \ V^dfl = — . (16) 

"0 Jo J«0 

The cosmic age is to = t(ao) = t(V>- The cosmic time and the 
lookback time to — t are shown in Fig. [TJ as functions of the 
redshift. 

Distances can be defined in many ways which typically lead 
to different expressions. The proper distance D pmp is the dis- 
tance measured by the light-travel time, thus 



dD prop = cdt 



D„ 



c 

Ho 



I 



da 
aE(a) 



(17) 



where the integral has to be evaluated between the scale fac- 
tors of emission and observation of the light signal. The co- 
moving radial coordinate w is the comoving distance mea- 
sured along a radial light ray. Since light rays propagate with 
zero proper time, ds = 0, this gives 



dw 



cdt 



w 



/da c f 
ad Ho J 



da 



i 2 E(a) 



(18) 
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parameter 


symbol 






WMAP-5 


comment 






alone 


+ BAO + SNe 














(reference) 




CMB temperature 


TcMB 


2.728 


+ 


0.004 K 




from dFixsen et al. 1996b 
l ii i 


total energy density 


^tot 


1.099 


+ 


0.100 
0.085 


1.0052 ± 0.0064 


matter density 


£2m0 


0.258 


+ 


0.03 


0.279 ± 0.015 


assuming spatial flatness 


baryon density 


^bO 


0.0441 


+ 


0.0030 


0.0462 ± 0.0015 


here and below 


cosmological constant 




0.742 


+ 


0.03 


0.721 ± 0.015 




Hubble constant 


h 


0.719 


+ 


0.026 
0.027 


0.701 ± 0.013 




power-spectrum normalisation 


0"8 


0.796 


+ 


0.036 


0.817 ± 0.026 




age of the Universe in Gyr 


k 


13.69 


± 


0.13 


13.73 ± 0.12 




decoupling redshift 


£dec 


1087.9 


± 


1.2 


1088.2 ± 1.1 




reionisation optical depth 


T 


0.087 


± 


0.017 


0.084 ± 0.016 




spectral index 


n s 


0.963 


+ 


0.014 
0.015 


960 + 0014 

u ,uu - 0.013 





TABLE I Cosmological parameters obtained from the 5-year data release of WMAP (Komatsu et al. , 2008 I, without and with the additional 
constraints imposed by baryonic acoustic oscillations (BAO, § |VII.A.6[ ) and type-la supernovae (SNe, § |fX) . The dimension-less Hubble 
constant h is defined in the normalisation parameter erg in J266[ l and the spectral index in j28}. Reionisation is briefly discussed in 
§ |X.B.2| Note that spatial flatness (K = 0) was assumed in deriving most of these values, as mentioned in last column. 
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FIG. 1 Cosmic time (solid curve) and lookback time (present age of 
the universe minus cosmic time) in Gyr as functions of redshift. Here 
and below, the parameters labelled as "reference" in Tab.|l]are used. 
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FIG. 2 Angular-diameter distance (solid curve) and proper distance 
in Mpc as functions of redshift. 



The angular-diameter distance D ang is defined such that the 
same relation as in Euclidean space holds between the physi- 
cal size of an object and its angular size. It turns out to be 

D mg {a) = af K [w(a)] , (19) 

with /k(w) given by The luminosity distance D\ um is anal- 
ogously defined to reproduce the Euclidean relation between 
the bolometric luminosity of an object and its observed flux. 
This results in 

Amg(a) f K [w(a)] 
Aum(a) = -, — = . (20) 

a L a 

Figure [2] shows the angular-diameter and the proper distance 
as functions of redshift. 



These distance measures rapidly diverge for scale factors 
a < 1. For small distances, i.e. for a « 1, they all reproduce 
the linear relation 



D(z) = . (21) 



Since time is finite in a universe with a Big Bang, any par- 
ticle can only be influenced by, and can only influence, events 
within finite regions, called horizons. Several different def- 
initions of horizons exist. They are typically characterised 
by some speed, e.g. the light speed, times the inverse Hub- 
ble function which sets the time scale. 
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C. Structures 

1 . Structure growth 

The hierarchy of cosmic structures is assumed to have 
grown from primordial seed fluctuations in the process of 
gravitational collapse: overdense regions attract material and 
grow. They are described by the density contrast 6(x, t) as a 
function of the comoving coordinates x, which is the density 
fluctuation relative to the mean density p(t), 



6(x, t) = 



p(x, t) - pi t) 
Pit) 



(22) 



Linear perturbation theory shows that, during the matter- 
dominated era, the density contrast 5 of sub-horizon pertur- 
bations is described by the second-order differential equation 
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FIG. 3 Fractional cosmic age t/t (solid curve) and structure growth 
as functions of redshift. 



S + 2H6 - AnGpS = 



(23) 



if the dark matter is cold, i.e. if dark-matter flows have neg- 
ligible velocity dispersion. Equation ( 23 1 has two solutions, 
one growing and one decaying. While the latter is irrelevant 
for structure growth, the growing mode is described by the 
growth factor D+(a), defined such that the density contrast at 
the scale factor a is related to the density contrast today 6q 
by 6(a) = 6oD + (a). In most cases of practical relevance, the 
growth factor is accurately described by 



D + (a) = 



G(a) 

cm 



(24) 



with the fitting formula 



G(a) = a£l n 



Qi /7 - n A + i + 



On 



1 + ^) 

70 / 



(25) 



where the density parameters have to be evaluated at the scale 
factor a. For a standard cosmological model, (25i deviates 
from the accurate solution by < 0.2 % for a e [0.01, 1]. The 
linear growth factor is shown in Fig. [3] together with the frac- 
tional cosmic age t/to. 

A very important length scale for cosmic structure growth 
is set by the horizon size at the end of the radiation-dominated 
phase. Structures smaller than that became causally connected 
while radiation was still dominating. The fast expansion due 
to the radiation density inhibited further growth of such struc- 
tures until the matter density started dominating. Small struc- 
tures are thus suppressed compared to large structures, which 
became causally connected only after radiation domination. 
The horizon radius u; eq at the end of the radiation-dominated 
era thus separates larger structures which could grow with- 
out inhibition, from smaller structures which were suppressed 
during radiation domination. In comoving coordinates, 



^ = 2(V2~1)^( 



a, 

Om0 / 



1/2 



(26) 



2. The power spectrum 

It is physically plausible that the density contrast in the Uni- 
verse is a Gaussian random field, i.e. that the probability for 
finding a density contrast between 6 and 6 + dd is given by 
a Gaussian distribution. The principal reason for this is the 
central limit theorem. A Gaussian random process is charac- 
terised by two numbers, the mean and the variance. By con- 
struction, the mean of the density contrast vanishes, such that 
the variance defines it completely. 

In linear approximation, density perturbations grow in 



place, as Eq. ( 23 1 shows because the density contrast at one 
comoving position x does not depend on the density contrast 
at another. As long as structures evolve linearly, their scale 
will be preserved, which implies that it is advantageous to 
study structure growth in Fourier rather than real space. The 
variance of the density contrast 6(k) at the comoving wave 
vector k in Fourier space is called the power spectrum 

{6(k)6*(k')) = (2nfP s (k) 6 D (k -1c'), (27) 

where the Dirac function £ D ensures that modes with different 
wave vectors are independent. 

Once the power spectrum is known, the statistical proper- 
ties of the linear density contrast are completely specified. It 
is a remarkable fact that two simple assumptions about the 
nature of the cosmic structures and the dark matter constrain 
the shape of P$(k) completely. If the rms mass fluctuation en- 
closed by the horizon is independent of time, and if the dark 
matter is cold, the power spectrum will behave as 



Ps(k) 



I*"- (fc«£ eq ) 
Ik"*- 4 (fc»£ eq ) 



(28) 



with the spectral index n s = 1 (Harrison 1970[|Peebles 1982 



Peebles and Yu|[197^|Wernber g 2008 ; Zeldovich|[1972[ >. The 
comoving wave number £ eq = 0.01 Mpc -1 of the peak location 
in the power spectrum is set by the comoving horizon radius 



at matter-radiation equality, w eq , given in (26 1. The steep de- 
cline for structures smaller than w eq reflects the suppression of 



6 



structure growth during radiation domination. Figure|4]shows 
the linearly and the nonlinearly evolved power spectrum of 
cold dark matter (CDM). 
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FIG. 4 Linearly and non-linearly evolved cold dark matter (CDM) 
power spectra. The linear CDM spectrum is shown for a = 1 and 
scales like D 2 + (a). 



3. Non-linear evolution 

As the density contrast approaches unity, its evolution be- 
comes non-linear. In the course of non-linear evolution, 
overdensities contract, causing matter to flow from larger to 
smaller scales. Power in the density-fluctuation field is thus 
transported towards smaller modes, or towards larger wave 
numbers k. This mode coupling process deforms the power 
spectrum on small scales, i.e. for large k. Detailed studies 
of the non-linear evolution of cosmic structures require nu- 
merical simulations, which need to cover large scales and to 
resolve small scales well at the same time. Much progress has 
been achieved in this field within the last two decades due to 
the fortunate combination of increasing computer power with 
highly sophisticated numerical algorithms, such as particle- 
mesh and tree codes, and adaptive mesh refinement techniques 



Barnes and Hut, 


1986; Efstathiou et al. 1985| Hockney and 


Eastwood| 1988; 


Monaghan 2005; Springel and Hernquist 


2002). 



The onset of non-linear evolution can be described by the 
so-called Zel'dovich approximation, which approximately de- 



scribes particle trajectories (Zel'dovich 1970| l. Although the 
Zel'dovich approximation breaks down as the non-linear evo- 
lution proceeds, it is remarkable for two applications. First, it 
allows a computation of the shapes of collapsing dark-matter 
structures and arrives at the conclusion that the collapse must 
be anisotropic, leading to the formation of sheets and fila- 



ments (Doroshkevich 19701. Filamentary structures thus ap- 
pear as a natural consequence of gravitational collapse in a 
Gaussian random field. Second, it provides an explanation for 
the origin of the angular momentum of cosmic structures. 



II. THE AGE OF THE UNIVERSE 
A. Nuclear cosmo-chronology 

1 . The age of the Earth 

How old is the Universe? We have no direct way to mea- 
sure how long ago the Big Bang happened, but there are vari- 
ous ways to set lower limits to the age of the Universe. They 
are all based on the same principle: since the Universe cannot 
be younger than any of its parts, it must be older than the old- 
est objects it contains. Three methods for age determination 
have been developed. One is based on the radioactive decay 
of long-lived isotopes, another constrains the age of globular 
clusters, and the third is based on the age of white dwarfs. 

Nuclear cosmo-chronology compares the measured abun- 
dances of certain long-lived radioactive isotopes with their 
initial abundances, which are eliminated by comparing abun- 
dances in different probes. 

To give a specific example, consider the two uranium iso- 
topes 235 U and 238 U. They both decay into stable lead iso- 
topes, 235 U — > 207 Pb through the actinium series and 238 U — > 
206 Pb through the radium series. The abundance of any of 
these two lead isotopes is the sum of the initial abundance, 
plus the amount produced by the uranium decay, e.g. 



N 2 m = N 2 m,o+N235 (e**'-l) 



(29) 



for 207 Pb, where /V235 is the abundance of 235 U nuclei today. 
A similar equation with 235 replaced by 238 and 207 replaced 
by 206 holds for the decay of 238 U to 206 Pb. The decay con- 
stants for the two uranium isotopes are measured as 



A235 = (1.015 Gyr)" 



A 2iS = (6.45 Gyr)- 



(30) 



The idea is now that ores on Earth or meteorites formed dur- 
ing a period very short compared to the age of the Earth f e , so 
that their abundances can be assumed to have been locked up 
instantaneously and simultaneously a time f e ago. Chemical 
fractionation has given different abundances to different sam- 
ples, but could not distinguish between different isotopes of 
the same element. Thus, we expect different samples to show 
different isotope abundances, but identical abundance ratios 
between different isotopes. 

The unstable lead isotope 204 Pb has no long-lived parents 
and is therefore a measure for the primordial lead abundance. 
Thus, the abundance ratios between 206 Pb and 207 Pb to 204 Pb 
calibrate the abundances in different samples. Suppose we 
have two independent samples a and b, in which the abun- 
dance ratios 



D _ N2O6 An- -^207 

"206 = — — ana « 2 07 = 77 — 

^V204 ^V204 



are measured. They are given by 



N 



R,. = j? w + _i-(e**-l) , 



(31) 



(32) 



204 



with (/, j) = (206, 238) or (207, 235). The present lead abun- 
dance ratios R206 and #207 should be the same in the two 
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samples and cancel when the difference between the ratios in 
the two samples is taken. Then, the ratio of differences can be 
written as 



R a — R b 
"207 "207 



7V 23 5 e A 



1 



"206 



"206 



/V 238 qIksU - i 



(33) 



Once the lead abundance ratios have been measured in the two 
samples, and the present uranium isotope ratio 



#235 
7V 2 ; 



= 0.00725 



(34) 
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is known, the age of the Earth t e is the only unknown in ( [33] l. 
This method yields ( |Patterson| 1 1 956| [Wasserburg et al. \ \\ 911) 



t e = (4.6 + 0.1) Gyr . 



(35) 



2. The age of the Galaxy 

A variant of this method can be used to estimate the age 
of the Galaxy, but this requires a model for how the radioac- 
tive elements were formed during the lifetime of the Galaxy 
until they were locked up in samples where we can measure 
their abundances today. Again, we can assume that the Galaxy 
formed quickly compared to its age f g . 

Suppose there was an instantaneous burst of star formation 
and subsequent supernova explosions a time f g ago and no 
further production thereafter. Then, the radioactive elements 
found on Earth today decayed for the time f g - f e until they 
were locked up when the Solar System formed. If we can 
infer from supernova theory the initially produced abundance 
ratio N235/N23&, we can conclude from its present value (34i 
and the age of the Earth what the age of the Galaxy must be. 

The situation is slightly more complicated because element 
production did not stop after the initial burst. Suppose that a 
fraction / of the total number 7V p of the heavy nuclei locked 
up in the Solar System was produced in a burst at t = 0, and 
the remaining fraction 1 — / was added at a steady rate until 
t — f g - f e when the Earth was formed. Given the constant 
production rate p, the abundance of a radioactive element with 
decay constant A is 



n = ar M + 



(36) 



before f g — t e , with a constant C to be suitably chosen, and 

N = N e- Al '- ( 's-'< )] (37) 

thereafter, where No is the abundance of the element locked 
up in the Solar System, as before. The initial conditions then 
require 



JV(0) = C + C = f N , 

A 



after the burst at t = 0, thus 



N(t g -t e ) = e 



fN P 



+ - e 
A V 







(38) 



(39) 



when the Earth formed, and 



We) 



/# p + f(e^-l) 



(40) 



today, when the galaxy reaches its age f g . Since the production 
rate must be 



P = 

the present abundance is 



(WW P 



N = 7V p e^ 



/ + 



(W) 
A(L - ? e ) 



( e *%-*.) _ !) 



(41) 



(42) 



in terms of the produced abundance 7V p . Supernova theory 
says that th e produced abundance r atio of the isotopes 235 U 



and 238 U is (Cowan et al. 



1987 



1991) 



A^235,p 
7^238,p 



1.4 + 0.2 



(43) 



Taking the ratio of ( 42 1 for the present abundances of 235 U and 
238 U, inserting the decay constants from (30 1, the abundance 



ratios from ( 34 1 and (43 1, and the age of the Earth f e from ( 35 1 



yields an equation which contains only the age of the Galaxy 
f g in terms of the assumed fraction /. This gives 



(6.3±0.2)Gyr /=1 
(8.0±0.6)Gyr / = 0.5 
(12±2)Gyr / = 



(44) 



More detailed models (|Cowan et al.\ |1987| |1991[ |Fowler| 



1989; Meyer and Schramm 1986) yield comparable results. 
Common assumptions and results from galaxy-formation the- 
ory assert that at least 1 Gyr is necessary before galactic 
disks could have been assembled. Therefore, nuclear cos- 
mochronology constrains the age of the Galaxy to fall approx- 
imately within 



7 Gyr <t g < 13 Gyr 



(45) 



Interestingly, nuclear cosmochronology has also been ap- 
plied to metal-poor and therefore presumably very old stars, 
in which heavy elements are nonetheless overabundant. For 
example, (Cayrel et al. 2Q01| > constrain the age of such a star 
to (12.5 + 3) Gyr. 



B. Stellar ages 

Another method for measuring the age of the Universe 
caused much trouble for cosmologists for a long time. It 
is based on stellar evolution and exploits the fact that the 
time spent by stars on the main sequence of the Hertzspmng- 
Russell or colour-magnitude diagram (see Fig. [5} depends 
sensitively on their mass and thus on their colour. 

Stars are described by the stellar-structure equations, which 
relate the mass M, the density p and the pressure p to the 
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radius r and specify the temperature T and the luminosity L. 
They are 



dp 
dr 



GMp 



dM 
~d7 



= Anr p , 



(46) 



which simply state hydrostatic equilibrium and mass conser- 
vation, and 



dr 

"d7 



3Li<p 



Anr 2 a^cT 3 



dL 
dr 



: 4nr pe . 



(47) 



which describe radiative energy transport and energy produc- 
tion, k is the opacity of the stellar material, e is the en- 
ergy production rate per unit mass, and osb is the Stefan- 
Boltzmann constant. Assuming k is independent of temper- 
ature, the energy-transport equation, mass conservation, hy- 
drostatic equilibrium and the equation-of-state for an ideal gas 
yield the scaling relations 



RT 4 , pM 

, M~pR\ p~ H — 

p R 



p~pT 



(48) 



The second pair of equations gives T ~ M/R, which yields 
L ~ M 3 when inserted into the first pair. 

The total lifetime t of a star on the main sequence must 
scale as Lt ~ M because the total energy radiated, Lt, must 
be a fraction of the total rest-mass energy. Together with the 
earlier result, we find 



M 
T 



M 



-2 



'-2/3 



(49) 



According to the Stefan-Boltzmann law, the star's luminosity 
must be 



R 2 T 4 



R 1 



M 3 

y4 



but we also know from above that T ~ M/R. Thus 



R- 



M i R 4 
M 4 



R 4 
M 



R ~ M 



1/2 



M 



1/2 



(50) 



(51) 



and the lifetime t on the main sequence turns out to scale as 
t ~ T~ l , by Eq. (49 1. Afterwards, stars move away from the 
main sequence towards the giant branch. Thus, as a coeval 
stellar population ages, the point in its Hertzspmng-Russell 
diagram up to which the main sequence remains populated 
moves towards lower luminosities and temperatures along 
(L, T) ~ (t _3 ^ 2 ,t _1 ). Old, coeval stellar populations exist: 
they are the globular clusters which surround the centre of 
the Galaxy in an approximately spherical halo. Therefore, the 
main-sequence turn-off points can be used to derive lower lim- 
its to the age of the Galaxy and the Universe. 

In practice, such age determinations proceed by adapting 
simulated stellar-evolution tracks to the Hertzspmng-Russell 
diagrams of globular clusters and assigning the age of the best- 



fitting stellar-evolution model to the cluster (|Renzini and Fusi 



Peccil [1988} |Straniero and Chieffi| [1991] |VandenBerg| |2000) . 
The simulated stellar-evolution tracks depend on the assumed 
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FIG. 5 Hertzspmng-Russell or colour-magnitude diagram of the 
globular cluster M 68. The magnitude is a logarithmic measure of the 
luminosity, colours are defined as magnitude differences. Stars popu- 
late the main sequence while burning hydrogen, then turn away from 
it towards the giant branch. The main-sequence turn-off is clearly 
visible (courtesy of Achim WeiB). 



metalicity of the stellar material, which changes the opacity 
and thus the energy transport through the stars. 

Since observations cannot tell the luminosity of the turn- 
off point on the main sequence, but only its apparent bright- 
ness, age determinations from globular clusters require obser- 
vations measuring both the apparent brightness of the main se- 
quence turn-off, and the distance to the globular cluster. There 
are several ways for estimating cluster distances. One uses 
the period-luminosity relation of a class of variable stars, the 
RR Lyrae stars, which are similar to the Cepheids that will 
play a central role in the next section. Another method ex- 
ploits that stars on the horizontal branch have a typical lumi- 
nosity and can thus be used to calibrate the cluster distance. 

Therefore, uncertainties in the distance determinations di- 
rectly translate to uncertainties in the age determinations. If 
the distance is overestimated, so is the luminosity, which im- 
plies that the age is underestimated, and vice versa. Globular 
clusters typically appeared to have ages well above estimates 



based on the cosmological parameters assumed (|Sandage and 



Cacciari 1990 Vandenberg 1983 1. In the past decade or so, 
this has changed because improvements in stellar-evolution 
theory and direct distance measurements with the Hipparcos 
satellite have lowered the globular-cluster ages, while recently 
determined cosmological parameters yield a higher age for the 



Universe as assumed before (Chaboyer et al. 1998 Salaris 



et a/.||1997a[ ). Now, globular-cluster ages imply 



t > (12.5 ± 1.3) Gyr 



(52) 



for the age of the Universe, allowing for a Gyr between the Big 
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Bang and the formation of the globular clusters (Chaboyer 
1998 ). The uncertainty in this result may still be considerable, 



though (see, e.g. |Grundahl et al. | [2000) 



C. Cooling of white dwarfs 

A third method for cosmic age determinations is based on 
the cooling of white dwarfs (see |Hansen and Liebert|20 03 for 
a review). They are the end products in the evolution of low- 
mass stars and form when the nuclear fuel in the stellar cores 
is exhausted. This happens typically when a core of carbon 
and oxygen has formed and temperature necessary for further 
fusion cannot be reached. The C-O core then shrinks until it is 
stabilised by electron degeneracy pressure. At this point, the 
mean free path of the electrons becomes effectively infinite 
and the core turns isothermal. 

The pressure of a degenerate electron gas is independent 
of temperature and depends only on the density. In the non- 
relativistic limit, p oc p 5 / 3 or 



M5/3 

R 5 



and £ 
R 



GMp GM 2 



R 1 - 



R 5 



(53) 



where the latter two equations express hydrostatic equilib- 
rium. Together, these relations imply R oc M~'/ 3 , showing 
that more massive white dwarfs have smaller radii. Their sur- 
face gravity, 



9' 



M 



M 5/3 



(54) 



which can be determined from their spectra, is thus a direct 
measure for their mass. It turns out that the vast majority of 
white dwarfs are born with similar masses, M w d w (0.55 + 



0.05) M Q (|Cool et al. 19961, which is a conspiracy of stellar 



evolution: more massive stars lose more mass through stellar 
winds before they turn into white dwarfs. 

Except for the latent heat of crystallisation, white dwarfs 
cool passively by radiating the thermal energy stored by their 
mass. Being born with essentially the same mass, they start 
with very similar amounts of thermal energy. Their tem- 
peratures and luminosities are determined by the radiation 
processes and the opacities of their atmospheres, which in 
turn depend on their metalicities. Age inferences from white 
dwarfs thus require that their metalicities be known and the 
energy transport through their atmospheres be modelled. 

After an initial short phase of rapid neutrino cooling, the 
energy loss slows down. This change in cooling time causes 
the white dwarfs to pile up at the temperature and the luminos- 
ity where the energy loss due to neutrinos falls below that due 
to radiative cooling. Thus, as time proceeds, their luminos- 
ity distribution develops a peak which is later slowly depopu- 
lated. Models for the white-dwarf cooling sequence allow the 
construction of time-dependent, theoretical luminosity distri- 
butions from which, by comparison with the observed lumi- 
nosity distribution, the age of a white-dwarf population can be 
determined. This methods yields 



for the age of white dwarfs in the Galactic disk (Chaboyer 
[T9^|Oswalt"eFar] [19961 |Salaris et q/.||1997b) . If we assume 
that massive spiral disks form at redshifts below z < 3, the 
implied age of the Universe is approximately 



f ~ (11 ± 1.4) Gyr. 



(56) 



The age of the white-dwarf population in the globular cluster 
M 4 has been measured to be ~ 12.1 Gyr (Ha nsen et a/.||2004] l. 



D. Summary 

Combining results, we see that the age of the Universe, as 
measured by its supposedly oldest parts, is at least > 1 1 Gyr, 
and this places serious cosmological constraints. In the frame- 
work of the Friedmann-Lemaitre models, they can be inter- 
preted as limits on the cosmological parameters. Suppose 
we lived in an Einstein-de Sitter universe with Q m o = 1 and 
£2ao = 0. Then, we know from ( fT6] > that 



fo= — > 11 Gyr 
ili 



H <2x 10~ 18 s 



(57) 



or Hq < 61 kms 1 Mpc" 1 in conventional units. The Hubble 
constant is measured to be larger than this, which can imme- 
diately be interpreted as an indication that we are not living in 
an Einstein-de Sitter universe. 



III. THE HUBBLE CONSTANT 

A. Hubble constant from Hubble's law 

1 . Hubble's law 

Slipher discovered preceding 1920 that distant galaxies typ- 



ically move away from us (Slipher 1927 1. Hubble found that 
their recession velocity grows with distance, 



v=H D, 



(58) 



and determined t he constant of proportion ality as Ho » 
559kms~ 1 Mpc 1 ( Hubble and Humason||l93l[ ). We had seen 
in ( |2"T| that all distance measures in a Friedmann-Lemaitre 
universe follow the linear relation 



D 



cz 



(59) 



to first order in z <s 1. Since cz — v is the velocity according 



t wi * (9.5+1) Gyr 



(55) 



to the linearised relation between redshift and velocity, ( 59 1 is 
exactly the relation Hubble found. 

The problem with any measurement of the Hubble constant 
from ( |58"] l is that, while ( f59| holds for an idealised, homoge- 
neous and isotropic universe, real galaxies have peculiar mo- 
tions on top of their Hubble velocity which are caused by the 
attraction from local density inhomogeneities. For instance, 
galaxies in our neighbourhood feel the gravitational pull of a 
cosmologically nearby supercluster called the Great Attractor 
and accelerate towards it. The galaxy M 3 1 in Andromeda is 



10 



moving within the Local Group of galaxies and approaches 
the Milky Way. 

Thus, the peculiar velocities of the galaxies must either be 
known well enough, for which a model for the velocity field is 
necessary, or they must be observed at so large distances that 
any peculiar motion is unimportant compared to their Hub- 
ble velocity. Requiring that the typical peculiar velocities, of 
order (300 . . . 600) km s _1 , be less than 10% of the Hubble ve- 
locity, galaxies with redshifts 



z > 10 x 



(300... 600) km s" 



0.01 ...0.02 



(60) 



need to be observed. This is already so distant that it is hard or 
impossible to apply simple geometric distant estimators. This 
illustrates why accurate measurements of the Hubble constant 
are so difficult: nearby galaxies, whose distances are more ac- 
curately measurable, do not follow the Hubble expansion well, 
but the distances to galaxies far enough to follow the Hubble 
law are very hard to measure (see (Sandage and Tammann 
|1974| > for an interesting account of the difficulties). 



2. The distance ladder 

Measurements of the Hubble constant from Hubble's law 
thus require accurate distance measurements out to cosmolog- 
ically relevant distance scales. Since this is impossible in one 
step, the so-called distance ladder must be applied, in which 
each step calibrates the next. 

The most fundamental distance measurement in the astro- 
physical toolbox is the trigonometric parallax caused by the 
annual motion of the Earth around the Sun. By definition, a 
star at a distance of a parsec perpendicular to the Earth's or- 
bital plane has a parallax of one arc second. Astrometric mea- 
surement accuracies of order 10~ 5 " are thus necessary to mea- 
sure distances of order lOkpc. Although the astrometric satel- 
lite Hipparcos could measure distances only to < 1 kpc, its 
results were used to calibrate the distance ladder with nearby 
stars. In this way, it was possible to determine the distance 
to the Large Magellanic Cloud as D LM c = (54.0 ± 2.5) kpc 
( [Feast ancT Catchpole 1997 1. An independent distance esti- 
mate based on the elliptical ring that appeared around the su- 
pernova r987Ajaye_DLMc = (46.77 + 0.04) kpc ( |Gould|1995 
see also Panagia et al. | 1991 1. 



The next important step in the distance ladder is formed 
by the Cepheids. These are stars in late evolutionary stages 
which undergo periodic variability. The underlying instability 
is driven by the temperature dependence of the atmospheric 
opacity in these stars, which is caused by the transition be- 
tween singly and doubly ionised Helium. The cosmologically 
important aspect of the Cepheids is that their variability period 
t and their luminosity L are related, 



Lcct 1 



(61) 



in the V band (|Feast and Walker) |1987[ |Madore and Freed-| 
1991 Udalski et al. 1999| l, hence their luminosity can 



be inferred from their variability period, and their distance 



from the ratio of their luminosity to the flux S observed from 
them, 



D 



L 

4^S 



(62) 



At the relevant distances, any distinction between differently 
defined distance measures is irrelevant. 

It is crucially important here that the calibration of the 
period-luminosity relation depends on the metalicity of the 
pulsating variables used, and thus on the stellar population 
they belong to. Hubble's originally much too high result for 
Ho was corrected when Baade realised that he had confused a 
metal-poor class of variable stars with the metal-rich Cepheids 
flBaade|[T956l >. 

Measuring the periods of Cepheids and comparing their ap- 
parent brightnesses in different galaxies thus allows the de- 
termination of the relative distances to the galaxies. For ex- 
ample, comparisons between Cepheids in the LMC and the 
Andromeda galaxy M 3 1 show ( Feast and Catchpole 1997] > 



D K 



Dlmc 



= 16.4 ± 1.1 , 



(63) 



while Cepheids in the member galaxies of the Virgo cluster 
yield (Freedman et al. 200 1| | 



D 



Virgo 



^LMC 



305 ± 16 . 



(64) 



Of course, for the Cepheid method to be applicable, it must 
be possible to resolve at least the outer parts of distant galaxies 
into individual stars and to reliably identify Cepheids among 
them. This was one reason why the Hubble Space Telescope 
was proposed (Freedman et al. 1994), to apply the superb 



resolution of an orbiting telescope to the measurement of Hq. 
Cepheid distance measurements are possible to distances < 
30Mpc. 

Scaling relations within classes of galaxies provide addi- 
tional distance indicators. In the three-dimensional parameter 
space spanned by the velocity dispersion <r D of their stars, the 
radius R e enclosing half the luminosity, and the surface bright- 
ness h within R e , elliptical galaxies populate the tight fun- 



damental plane (Djorgovski and Davis 1987] Dressier et al. 



1987i defined by 



R R oco-l A r * 5 



Since the luminosity must be 



L oc I e R* 



the fundamental-plane relation implies 



L <* o?V 7 



(65) 



(66) 



(67) 



Such a relation follows directly from the virial theorem if the 
mass-to-light ratio in elliptical galaxies increases gently with 



mass, 



M 



oc M 



0.2 



(68) 
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Thus, if it is possible to measure the surface brightness I e 
and the velocity dispersion cr v of an elliptical galaxy, the 
fundamental-plane relation gives the luminosity, which can be 
compared to the flux to find the distance. 

An interesting way for determining distances to galaxies 
uses the fluctuations in their surface brightness (Ton ry and| 



Schneider 1988 Tonry et al. 1997 2000 1. The idea behind 



this method is that the fluctuations in the surface brightness 
will be dominated by the rare brightest stars, and that the opti- 
cal luminosity of the entire galaxy will be proportional to the 
number N of such stars. Assuming Poisson statistics, the fluc- 
tuation level will be proportional to yfN, from which N and 
Loc N can be determined once the method has been calibrated 
with galaxies whose distance is known otherwise. Again, the 
distance is then found by comparing the flux to the luminosity. 

Planetary nebulae, which are late stages in the evolution 
of massive stars, have a luminosity function with a steep up- 
per cut-off ( |Hui et al.\ |1993| |Jacoby| |1980| l. Moreover, their 
spectra are dominated by sharp nebular emission lines which 
facilitate their detection even at large distances because they 
appear as bright objects in narrow-band filters tuned to the 
emission lines. Since the cut-off luminosity is known, it can 
be converted to a distance. 

An important class of distance indicators are supernovae of 



type la. As will be described in § IX. B their luminosities can 
be standardised, allowing their distances to be inferred from 
their fluxes. 

Although they are not standard (or standardisable) candles, 
core-collapse supernovae of type II can also be used as dis- 



tance indicators through the Baade-Wesselink method (Baade 



1926 Wesselink 1946[ ). Suppose the spectrum of the su 



pernova photosphere can be approximated by a Planck curve 
whose temperature can be determined from the spectral lines. 
Then, the Stefan-Boltzmann law fixes the total luminosity, 



L = cr sbR T 



(69) 



The photospheric radius, however, can be inferred from the 
expansion velocity of the photosphere, which is measurable 
by the Doppler shift in the emission lines, times the time af- 
ter the explosion. When applied to the supernova SN 1987A 
in the Large Magellanic Cloud, the Baade-Wesselink method 
yields a distance of 



D LM c = (49 + 3)kpc 



(70) 



based on infrared photometry (Schmidt et al. 



1992 1 which 



agrees with the direct distance measurements (Gould |1995| l. 
Note that the Bade-Wesselink method gives a one-step geo- 
metrical distance indicator bypassing the distance ladder. 



3. The HST Key Project 

All these distance indicators were used in the HST Key 
Project (Freed man et al.\ |2001| l to determine accurate dis- 
tances to 26 galaxies between (3...25)Mpc, and five very 
nearby galaxies for testing and calibration. 

Double-blind photometry was applied to the identified dis- 
tance indicators. Since Cepheids tend to lie in star-forming 



regions and are thus attenuated by dust, and since their period- 
luminosity relation depends on metalicity, respective correc- 
tions had to be carefully applied. Then, the measured veloci- 
ties had to be corrected by the peculiar velocities, which were 
estimated by a model for the flow field. The estimated lu- 
minosities of the distance indicators could then be compared 
with the extinction-corrected fluxes to determine distances, 
whose proportionality with the velocities corrected by the pe- 
culiar motions finally gave the Hubble constant. A weighed 
average over all distance indicators is (Freed man et a/.|[2001[ ) 



H = (72 + 8)kms" 1 Mpc" 



(71) 



where the error is the square root of the systematic and statisti- 
cal errors summed in quadrature. An alternative interpretation 
of the data is summarised in (Sandage et al. , 2006), who find 
a 14 % smaller value for Ho because of a different calibration. 



B. Gravitational Lensing 

A completely different method for determining the Hubble 
constant uses gravitational lensing. Masses bend passing light 
paths towards themselves and therefore act similarly to convex 
glass lenses. As in ordinary geometrical optics, this effect can 
be described applying Fermat's principle to a medium with an 
index of refraction ( |Schneider et al.\ |1992| l 



1 



2<D 

72" 



(72) 



where <b is the Newtonian gravitational potential. 

If it is strong enough, the bending of the light paths causes 
multiple images to appear from single sources. Compared to 
the straight light paths in absence of the deflecting mass dis- 
tribution, the curved paths are geometrically longer, and they 
have to additionally propagate through a medium whose in- 
dex of refraction is n > 1. This gives rise to a time delay 
( Blandford and Narayan 1986]> which has a geometrical and 



a gravitational component, 
1 + z d D d D A 



(73) 



where 9 are angular coordinates on the sky and /3 is the angu- 
lar position of the source, if/ is the appropriately scaled New- 
tonian gravitational potential of the deflector, projected along 
the line-of-sight. According to Fermat's principle, images oc- 
cur where r is extremal, i.e. V#t = 0. The prefactor contains 
the distances D^ s .ds from the observer to the deflector, to the 
source, and from the deflector to the source, respectively, and 
the redshift Zd of the lens. 

The projected lensing potential iff is related to the surface- 
mass density E of the deflector by 



VV = 2 — = 2k , 
y 

where the critical surface-mass density is 



AnG D d Dd s 



(74) 



(75) 
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The dimension-less time delay t from (73 1 is related to the 
true physical time delay t by 



t oc 



Ho 



(76) 



where the proportionality constant is a dimension-less com- 



with the Hubble radius cH. 



o 



bination of the distances 

and the deflector redshift 1 + z&- Equation f76) shows that the 
true time delay is proportional to the Hubble time, as it can 



intuitively be expected (Refsdal 1964 1. 



Time delays are measurable in multiple images of a vari- 
able source. The variable signal arrives after different times in 
the images seen by the observer. If the deflector is a galaxy, 
time delays are typically of order days to months and therefore 
observable with a reasonable monitoring effort. 

Interestingly, it can be shown in an elegant, but lengthy 



calculation ( Kochanek 2002 1 that measured time delays can 
be inverted to find the Hubble constant from the approximate 
equation 



#o * A(l -<*» + £<*>(!/ -1) 



(77) 



where A and B are constants depending on the measured im- 
age positions and time delays, (k) is the mean scaled surface- 
mass density of the deflector averaged within an annulus 
bounded by the image positions, and rj a 2 is the logarith- 
mic slope of the deflector's density profile. 

Therefore, if a model exists for the gravitationally-lensing 
galaxy, the Hubble constant can be found from the positions 
and time delays of the images. Applying this technique to five 
different lens systems, (Kochanek 2002) found 



H = (73 ±8) km s" 1 



(78) 



assuming that the lensing galaxies have radially constant 
mass-to-light ratios. 

This result is remarkable because it was obtained in one 
step without any reference to the extragalactic distance lad- 
der. On the other hand, it is also problematic because it is ob- 
tained only if significantly more concentrated density profiles 
of the lensing galaxies are assumed than obtained e.g. from 
the kinematics of the stars in galaxies. If less concentrated 
lenses are adopted which agree with the kinematic constraints, 
Hq derived from gravitati onal time delays d rops substantially 



to values near 50km s (Kochanek 20031. This hints at an 



as yet unexplained discrepancy between measurements of Ho 
and the measured time delays within the CDM framework. 



C. Summary 

The results given so far are in excellent agreement with the 
value 



H = (70.1 ± 1.3)kms _1 Mpc" 1 



(79) 



derived from CMB measurements (see § [VI] and Tab. [I| as- 
suming spatial flatness, K — 0. adopting it, we can calibrate 



several important numbers that scale with some power of the 
Hubble constant. First, in cgs units, the Hubble constant is 



Ho = (2.26 + 0.04) x 10~ 1S s , 
which implies the Hubble time 



m l = (14.01 +0.26) Gyr 



and the Hubble radius 

c 



(80) 



(81) 



= (1 .327 ± 0.025) x 10 28 cm = (4.299 + 0.08) Gpc . (82) 



H 



i) 



The critical density of the Universe turns out to be 
3H 2 

Pcro = = (9.15+0.3) x 10- 30 gcnr 3 . (83) 

It corresponds to five protons per m 3 , or approximately one 
galaxy per Mpc 3 . 



IV. BIG-BANG NUCLEOSYNTHESIS 

A. The origin and abundance of Helium-4 

1 . Elementary considerations 

Stellar spectra show that the abundance of 4 He in stellar 
atmospheres ranges between 0.2 < Y < 0.3 by mass (the Sun 
has Y = 0.263), i.e. about a quarter of the baryonic mass in 
the Universe is composed of 4 He. It is produced in stars in the 
course of hydrogen burning. Per 4 He nucleus, the amount of 
energy released corresponds to 0.7% of the masses involved, 
or 

AE = Amc 2 = 0.007 (2m p + 2m n )c 2 * 0.028 m p c 2 

* 26 MeV * 4.2 x 10~ 5 erg . (84) 

Suppose a galaxy such as ours, the Milky Way, shines with 
a luminosity of L w 10 I0 L G ~ 3.8 x 10 43 ergs _1 for a good 
fraction of the age of the Universe, say for r = 5 x 10 9 yr as 
1.5 x 10 17 s. Then, it releases a total energy of 



Lt ^ 5.7 x 10 60 er 



(85) 



The number of 4 He nuclei required to produce this energy is 

AjV = — ~ 5 - 7xl ° 6 ° ^l,4xl0 65 , (86) 
AE 4.2 x 10- 5 

which amounts to a 4 He mass of 



M He ~ Am v AN ~ 9.3 x 10 41 g 



(87) 



Assume further that the galaxy's stars were all composed of 
pure hydrogen initially, and that they are all more or less 
similar to the Sun. Then, the mass in hydrogen was Mu ~ 
10 1() M @ s 2x 10 43 g initially, and the final 4 He abundance by 
mass expected from the energy production is 



9.3 x 10 41 r „ 



(88) 
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much less than the 4 He abundance actually observed. This 
discrepancy is exacerbated by the fact that 4 He is destroyed 
in later stages of the evolution of massive stars, a process af- 
fected also by mixing in stellar interiors. 

We thus see that the amount of 4 He observed in stars was 
highly unlikely produced by these stars themselves under rea- 
sonable assumptions during the lifetime of the galaxies. We 
must therefore consider that most of the 4 He which is now 
observed may have existed already well before the galaxies 
formed. 

Nuclear fusion of 4 He and similar light nuclei in the early 
Universe is possible only if the Universe was hot enough for a 
sufficiently long period during its early evolution. The nuclear 
binding energies of order ~ MeV imply that temperatures had 
to have fallen below T ~ 10 6 MeV x 1.16 x 10 4 K/MeV = 
1.16 x 10 10 K before nucleosynthesis could begin. On the 
other hand, temperatures needed to be still high enough for 
charged nuclei to surmount the Coulomb barrier. Typical 
temperatures during primordial nucleosynthesis were of order 
kT < 100 keV. Since the temperature of the (photon back- 
ground in the) Universe is now 7o ~ 3 K as we shall see later, 
this corresponds to times when the scale factor of the Universe 



1.16 x 10 9 



2.59 x 10" 



(89) 



At times so early, the actual mass density and a possible 
cosmological constant are entirely irrelevant for the expansion 
of the Universe, which is only driven by the radiation density. 
Thus, the expansion function can be simplified to read E(a ) = 



Q}^a 2 , and we find for the cosmic time according to 



15 i 



t(a) 



~TJ7 — f a ' da ' 
Clll 2 H Jo 



2Q r 1 ( ( 2 // 



2.40 x 10 1 V s 



(90) 

where we have inserted the Hubble constant from ( f79l and the 
radiation-density parameter today, £2 r o ~ 8.51 x 10~ 5 , which 
will be justified later [see ( |157| >]. It contains the contributions 
from photons and three neutrino species, which are the parti- 
cles relativistic at the relevant time. 

Inserting a nuc from ( [89] > into ( |90| yields a time scale for nu- 
cleosynthesis of order a few minutes. It is instructive for later 
purposes to establish a relation between time and temperature 
based on ((90]). Since T = T /a, 

t = 2.40 x 10" (ffs. 0.89 yy~ 2 S - (91) 



efficiently enough for the subsequent fusion of 4 He, but not 
too efficiently because otherwise too much deuterium would 
be left over after the end of nucleosynthesis. Realising this, 
(Gamow 1948) suggested that the amount of deuterium pro- 
duced had to be "just right", which he translated into the intu- 
itive criterion 



n^{crv)t « 1 



(93) 



where «b is the baryon number density, (cry) is the velocity- 



averaged cross section for the reaction ( 92 1, and t is the time 
available for the fusion, which we have seen in ( |9T| ) to be 
set by the present temperature of the cosmic radiation back- 
ground, 7q, and the temperature T required for deuterium fu- 
sion. 

Thus, from an estimate of the baryon density «b in the Uni- 
verse, from the known velocity-averaged cross section (cry), 
and from the known temperature required for deuterium fu- 
sion, Gamow's criterion enables an estimate of the present 
temperature 7q of the cosmic radiation background. Based 



on this, (Alpher and Herman 1949[ ) were able to predict 
To » (1...5)K! After these remarkably simple and far- 
reaching conclusions, we shall now study primordial nucle- 
osynthesis and consequences thereof in more detail. 



3. Elements produced 



The fusion of deuterium (92 1 is the crucial first step. Since 
the photodissociation cross section of d is large, destruction of 
d is very likely because of the intense photon background until 
the temperature has dropped way below the binding energy of 
d, which is only 2.2 MeV, corresponding to 2.6 x 10 10 K. In 
fact, substantial d fusion is delayed until the temperature falls 
to T = 9xl0 8 Kor/tr « 78keV! As |9l) shows, this happens 
t « 150 s after the Big Bang. 

From deuterium, 3 He and tritium t can be built, which can 
both be processed to 4 He. These reactions are now fast, im- 
mediately converting the newly formed d. In detail, these re- 
actions are 



d + p 
d + d 
d + d 
3 He + « 



3 He + y , 
3 He + n , 
t + p , and 
t + p , 



(94) 



followed by 



2. The Gamow criterion 

A crucially important step in the fusion of 4 He is the fusion 
of deuterium d, 



p + n — > d + y 



(92) 



because the direct fusion of 4 He from two neutrons and two 
protons is extremely unlikely. The conditions in the early uni- 
verse must have been such that deuterium could be formed 



3 He + J 
t + d 



+ He + p 
*He + n 



and 



(95) 



Fusion reactions with neutrons are irrelevant because free 
neutrons are immediately locked up in deuterons once deu- 
terium fusion begins, and passed on to t, 3 He and 4 He in the 
further fusion steps. 

Since there are no stable elements with atomic weight A = 
5, addition of protons to 4 He is unimportant. Fusion with d 
is unimportant because its abundance is very low due to the 



14 



efficient follow-up reactions. We can therefore proceed only 
by fusing 4 He with t and 3 He to build up elements with A = 7, 

t + 4 He -» 7 Li + y , 
3 He + 4 He — > 7 Be + y , followed by 
7 Be + e~ 7 Li + v e . (96) 



Some 7 Li is destroyed by 



7 Li + /?^2 4 He 



(97) 



The fusion of two 4 He nuclei leads to 8 Be, which is unstable. 
Further fusion of 8 Be in the reaction 



8 Be + 4 He 



12. 



C + y 



(98) 



is virtually impossible because the low density of the reaction 
partners essentially excludes that a 8 Be nucleus meets a 4 He 



nucleus during its lifetime. While the reaction (98 i is possi- 
ble and extremely important in stars, it is suppressed below 
any importance in the early Universe. The absence of stable 
elements with A = 8 thus prohibits any primordial element 
fusion beyond 7 Li. 



4. Helium abundance 

Once stable hadrons can form from the quark-gluon plasma 
in the very early universe, neutrons and protons are kept in 
thermal equilibrium by the weak interactions 



<-> n + v e 



n + e <-» p + v e 



(99) 



until the interaction rate falls below the expansion rate of the 
Universe. While equilibrium is maintained, the abundances 
n„ and n p are controlled by the Boltzmann factor 

HSf • <100) 

where Q = 1 .3 MeV is the energy equivalent of the mass dif- 
ference between the neutron and the proton. 

The weak interaction freezes out when T « 10 10 K or kT =s 
0.87 MeV, which is reached f * 2s after the Big Bang. At this 
time, the n abundance by mass is 



X„(0) = 



n„ + n, : 



l + exp^jj ,0.17 
(101) 

Afterwards, the free neutrons undergo /3 decay with a half life 1 
ofT„ = (885.7 + 0.8) s, thus 



X„ = X„(0)e 



-t/T n 



(102) 



When d fusion finally sets in at tj « 150 s after the Big Bang, 
the neutron abundance has dropped to 



X n (t d ) « X„(0)e- 



-tdhn 



0.14 



(103) 



1 Particle Data Group, http://pdg.lbl.gov/ 



Now, essentially all these neutrons are collected into 4 He be- 
cause the abundances of the other elements can be neglected 
to first order. This simple estimate yields a 4 He abundance by 
mass of 



Y ~ 2X„{t d ) * 0.28 



(104) 



because the neutrons are locked up in pairs to form 4 He nu- 
clei. The Big-Bang model thus allows the prediction that 4 He 
must have been produced such that its abundance is approxi- 
mately 28% by mass, which is in remarkable agreement with 
the observed abundance and thus a strong confirmation of the 
Big-Bang model. 



5. Expected abundances and abundance trends 

Precise abundances of the light elements as produced by the 
primordial fusion must be calculated solving rate equations 
based on the respective fusion cross sections. Uncertainties 
involved concern the exact values of the cross sections and 
their energy dependence, and the precise life time of the free 
neutrons. Since primordial nucleosynthesis happens during 
the radiation-dominated era, the expansion rate is exclusively 
set by the radiation density. Then, the only other parameter 
controlling the primordial fusion processes is the baryon den- 
sity, if the neutron lifetime is taken for granted. 



An excellent recent review of these matters is (Steigman 
2007), to which we refer for further detail. 

Ignoring possible modifications by non-standard tempera- 
ture evolution, the only relevant parameter defining the pri- 
mordial abundances is the ratio between the number densities 
of baryons and photons. Since both densities scale like cT 3 or, 
equivalently, like r 3 , their ratio r\ is constant. Anticipating the 
photon number density to be determined from the temperature 
oftheCMB, 



«B 



io-'V, 



i/io = 273fi B /* 2 



(105) 



Thus, once we know the photon number density, and once 
we can determine the parameter rj from the primodial element 
abundances, we can infer the baryon number density. Typi- 
cal 2-cr uncertainties at a fiducial 77 parameter of 7710 = 5 are 
0.4% for 4 He, 15% for d and 3 He, and 42% for 7 Li. The 4 He 
abundance depends only very weakly on 77 because the largest 
fraction of free neutrons is swept up into 4 He without strong 
sensitivity to the detailed conditions. 

The principal effects determining the abundances of d, 3 He 
and 7 Li are the following: with increasing 77, they can more 
easily be burned to 4 He, and so their abundances drop as 77 in- 
creases. At low 77, an increase in the proton density causes 7 Li 



destruction in the reaction ( 97 1, while the precursor nucleus 
7 Be is more easily produced if the baryon density increases 
further. This creates a characteristic "valley" of the predicted 
7 Li abundance near 77 m (2 . . . 3) x 10" 10 . 
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B. Observed element abundances 



1. Principles 



ties, 3 He is commonly excluded from primordial abundance 
measurements. 



Of course, the main problem with any comparison between 
light-element abundances predicted by primordial nucleosyn- 
thesis and their observed values is that much time has passed 
since the primordial fusion ceased, and further fusion pro- 
cesses have happened since. Seeking to determine the primor- 
dial abundances, observers must therefore either select objects 
in which little or no contamination by later nucleosynthesis 
can reasonably be expected, in which the primordial element 
abundance may have been locked up and separated from the 
surroundings, or whose observed element abundances can be 
corrected for their enrichment during cosmic history in some 
way. 

Deuterium can be observed in cool, neutral hydrogen gas 
(HI regions) via resonant UV absorption from the ground 
state, or in radio wavebands via the hyperfine spin-flip tran- 
sition, or in the sub-millimetre regime via dH molecule lines. 
These methods all employ the fact that the heavier d nucleus 
causes small changes in the energy levels of electrons bound 
to it. 3 He is observed through the hyperfine transition in its 
ion 3 He + in radio wavebands, or through its emission and ab- 
sorption lines in HII regions. 4 He is of course most abundant 
in stars, but the fusion of 4 He in stars is virtually impossible 
to be corrected precisely. Rather, 4 He is probed via the emis- 
sion from optical recombination lines in HII regions ( |Izotov| 



and Thuan 



2004). Measurements of Li must be performed 



in old, local stellar populations. This restricts observations 
to cool, low-mass stars because of their long lifetime, and to 
stars in the Galactic halo to allow precise spectra to be taken 
despite the low 7 Li abundance. 



2. Evolutionary corrections 



Stars brooded heavy elements as early as z ~ 6 or even 



higher (Ciardi and Ferrara 2005 ). Any attempts at measuring 
primordial element abundances must therefore concentrate on 
gas with a metal abundance as low as possible. The depen- 
dence of the element abundances on metalicity allows extrap- 
olations to zero enrichment. 

Such evolutionary corrections are low for d because it is ob- 
served in the Lyman-or forest lines, which arise from absorp- 
tion in low-density, cool gas clouds at high redshift. Likewise, 
they are low for the measurements of 4 He because it is ob- 
served in low-metalicity, extragalactic HII regions. Probably, 
little or no correction is required for the 7 Li abundances deter- 



mined from the spectra of very metal -poor halo stars (Char- 
bonnel and Primas 2005 ), but a lively discussion is going on 
( [Asplund et aZ.||2006||Korn et a/.||2006 1. 

Inferences from 3 He are different because 3 He is produced 
from d in stars during the pre-main sequence evolution. It 
is burnt to 4 He during the later phases of stellar evolution in 
stellar cores, but conserved in stellar exteriors. Observations 
indicate that a net destruction of 3 He must happen, possibly 
due to extra mixing in stellar interiors. For these uncertain- 



3. Specific results 

Due to the absence of strong evolutionary effects and its 
steep monotonic abundance decrease with increasing jj, d is 
perhaps the most trustworthy baryometer. Since it is produced 
in the early Universe and destroyed by later fusion in stars, all 
d abundance determinations are lower bounds to its primordial 
abundance. Measurements of the d abundance at high redshift 
are possible through absorption lines in QSO spectra, which 
are likely to probe gas with primordial element composition 
or close to it. Such measurements are challenging in detail 
because the tiny isotope shift in the d lines needs to be distin- 
guished from velocity-shifted hydrogen lines, H abundances 
from saturated H lines need to be corrected by comparison 
with higher-order lines, and high-resolution spectroscopy is 
required for accurate continuum subtraction. 

A deuterium abundance of 



(2-68!°;g) x 10- 



(106) 



in high-redshift QSOs relative to hydrogen appears consistent 
with most relevant QSO spectra, although substantially lower 
values have been derived (e.g. |Pettini and Bowen 2001 1. A 
substantial depletion from the primordial value is unlikely be- 
cause any depletion should be caused by d fusion and thus 
be accompanied by an increase in metal abundances, which 
should be measurable. 

Some spectra which were interpreted as having < 10 times 
the d abundance from ( |107[ > may be due to lack of spectral 
resolution. The d abundance in the local interstellar medium 
is typically lower, 



i_hi_ 



1.5 x 10~ 5 



(107) 



which is consistent with d consumption due to fusion pro- 
cesses. Conversely, the d abundance in the Solar System is 
higher because d is locked up in the ice on the giant planets. 

In low-metalicity systems, 4 He should be near its primor- 
dial abundance, and a metalicity correction can be applied. 
Possible systematic uncertainties are due to modifications by 
underlying stellar absorption, collisional excitation of ob- 
served recombination lines, and the exact regression towards 
zero metalicity. Values between Y p = 0.240 + 0.006 and 



T p = 0.2477 + 0.0029 are considered realistic and likely (Pe 
|imbertef a/!] [2007) . 

Observations of the 7 Li abundance aim at the oldest stars 
in the Galaxy, which are halo (Pop-II) stars with very low 
metalicity. They should have locked up very nearly primor- 



dial gas, but may have processed it (Pinsonneault et al. 2002 



Salaris and Weiss 2.00 1 1 > . Cool stellar atmospheres are diffi- 



cult to model, and Li may have been produced by cosmic -ray 



spallation on the interstellar medium (Steigman and Walker 
\T992) . 

In the limit of low stellar metalicity, the observed 7 Li abun- 
dances group around an approximately flat line (the so-called 
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Spite plateau ( |Spite and Spite 1982 )), which is asymptotically 
independent of metalicity, 



A( 7 Li) = 12 + log(n Li /n H ) = 2.2 ± 0.1 



(108) 



and shows very little dispersion. Stellar rotation is important 



here because it enhances mixing in stellar interiors (Pinson- 
|neaultef al.\\2M>) . 

The Spite plateau is unlikely to reflect the primordial 7 Li 
abundance, but corrections are probably moderate. A possible 
increase of 7 Li with the iron abundance indicates low produc- 
tion of 7 Li, but the probable net effect is a depletion with re- 
spect to the primordial abundance by no more than ~ 0.2 dex. 
A conservative estimate yields 



2.1 < A( 7 Li) < 2.5 . 



(109) 



In absence of depletion, this value falls into the valley ex- 
pected in the primordial 7 Li at the boundary between destruc- 
tion by protons and production from 8 Be. However, if 7 Li was 
in fact depleted, its primordial abundance was higher than the 
value i 



109 1, and then two values for 77^ are possible. 



4. Summary of results 



Through the relation ( 105 I, the density of visible baryons 



alone implies 7710 > 1.5. The deuterium abundance derived 
from absorption systems in the spectra of high-redshift QSOs 
indicates 7710 = 5.65... 6. 38. The 7 Li abundance predicted 
from this value of 77 is then A( 7 Li) p = 3.81 . . .4.86, which is 
at least in mild conflict with the observed value A( 7 Li) = 2.1- 
2.3, even if a depletion by 0.2 dex due to stellar destruction is 
allowed. 

The predicted primordial abundance of 4 He is then Y p = 
0.2477 . . . 0.2489, which overlaps with the measured value 
Yp = 0.228 ... 0.248. Thus, the light-element abundances 
draw a consistent picture starting from the observed deuterium 
abundance. 

We thus find that Big-Bang nucleosynthesis alone implies 

0.0207 < Q B /7 2 < 0.0234 or 0.0399 < O b < 0.045 

(HO) 

if conclusions are predominantly based on the deuterium 
abundance in high-redshift absorption systems. We shall later 
see that this result is in fantastic agreement with independent 
estimates of the baryon density obtained from the analysis of 
structures in the CMB. 

At this point, it is important to note that primordial nu- 
cleosynthesis depends on the number of relativistic particle 
species, and thus on the number of light neutrino families. The 
theoretical abundances given here assume three lepton fami- 
lies in the standard model of the electroweak interaction. In 
fact, light-element abundances had been used to constrain the 
number of lepton families before it was measured in particle 




FIG. 6 Theoretically expected abundances of the light ele ments deu- 
terium, 3 He, 4 He and 7 Li as functions of 77,0 (following jSteigman 



2007 )). The observed deuterium abundance is marked by the two 
horizontal solid lines. 



V. THE MATTER DENSITY IN THE UNIVERSE 
A. Mass in galaxies 

1. Stars 

Given the luminosity of a stellar population, what is its 
mass? If all stars were like the Sun, the answer was trivial, 
but this is not the case. We shall focus the discussion here on 
stars which fall on the main sequence of the colour-magnitude 
diagram. Stars are formed with an initial mass distribution, the 
"initial mass function", which is for our purposes reasonably 



approximated by the Salpeter form (Salpeter 1955 1 



dlnM 



oc M 



■1.35 



(HI) 



Expressing the mass M in solar units, m = M/M e , and nor- 
malising the mass distribution to unity in the mass range 
mo < m < 00 yields 



dN /m () \ 235 dm 

= 1.35_£ _ 

dm \ m I mo 



(112) 



detectors (Steigman et al. 1977 Yang et al. 1984 1. 



The accepted lower mass limit for a star is mo = 0.08 be- 
cause nuclear fusion cannot set in in objects of lower mass. 
However, we are interested in stars that are able to produce 
visible or near-infrared light so that we can translate the re- 
spective measured luminosities into mass. For a simple esti- 
mate, assume that stars have approximate Planck spectra, for 
which Wien's law holds, relating the wavelength A max of the 
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peak in the Planck curve to the temperature T, 



At the same time, (II 151) needs to be satisfied, hence 



A max = 0.29 I ^ I cm . 



(113) 



The effective Solar temperature T = 5780 K implies A max = 
5.0 x 10~ 5 cm. Stars releasing the majority of their energy 
in the optical and near-infrared regime should have /l max :$ 
lum = 10 4 cm and thus T > 2900 K * 0.5 T @ . We saw in 
( plj l that the temperature scales as T ~ M 1 ^ 2 ; thus, T > 0.5 T Q 
implies ra m 0.25. 



We also saw following (48 1 that the luminosity scales as 



L ~ M . With / = L/Lq, the mean mass-to-light ratio of the 
visible stellar population is thus expected to be 



dN ( m )dm- r 
dm\l ) J 



dN dm 

dm m 2 



1.35 
3.35 ml 



6.4 . 
(114) 

This shows that an average stellar population visible in the op- 
tical and near-infrared spectral ranges is expected to require 
« 6.4 solar masses for one solar luminosity. In order to pro- 
duce, say, 10 10 L G , a galaxy thus needs to have a mass of at 
least as 6.4 x 10 10 M . This estimate assumes a homogeneous 
stellar population at the beginning of its evolution, neglects 
corrections from the exact spectral-energy distribution of the 
stars and the emission of giant stars, but (m/l) > 6.4 still 
seems appropriate for evolved populations. 



2. Galaxies 

The rotation velocities of stars orbiting in spiral galaxies 
are observed to rise quickly with radius and then to remain 
roughly constant ( Sofue and Rubinl |2001[). If measurements 



are continued with neutral hydrogen beyond the radii out to 
which stars can be seen, these rotation curves are observed to 
continue at an approximately constant level (Kre gel and van| 
|derKruitl|2004) . 

In a spherically-symmetric mass distribution, test particles 
on circular orbits have orbital velocities of 



2 , , GM{r) 



Flat rotation curves thus imply that M(r) oc r, requiring the 
density to drop like p(r) oc r~ 2 . This is much flatter than the 
light distribution, which shows that spiral galaxies are charac- 
terised by an increasing amount of dark matter as the radius 
increases. 

A mass distribution with p oc r~ 2 has formally infinite mass, 
which is physically implausible. However, at finite radius, the 
density of the galaxy falls below the mean density of the sur- 
rounding universe. The spherical collapse model (Gunn and 



|Gott||1972| often invoked in cosmology shows that a spherical 
mass distribution can be considered in dynamical equilibrium 
if its mean overdensity is approximately 200 times higher than 
the mean density p. Let R be the radius enclosing this over- 
density, and M the mass enclosed, then 



M _ 3M 



200p 



M _ SOOnpR 2 
~R ~ 3 



(116) 



SOOnpR 2 



"rot 

G 



R 



I 3i4 \ 
\%QQjtGp) 



1/2 



Inserting typical numbers yields 
R = 290 kpc ( 



200 km s" 



With ( 1 15 i, this implies 



Rip" 

M= — ™ =2.7x 10 12 M G 



\200kms- 1 ; 



(117) 



(118) 



(119) 



G ""\200kms- 1 
Typical luminosities of spiral galaxies are given by the 



Tully-Fisher relation (Tully and Fisher 1977 1, 

L - L { "*» f " 4 
*\220kms- 1 / 



(120) 



with L„ * 2.4 x 10 10 L Q . Thus, the mass-to-light ratio of a 
massive spiral galaxy is found to be 

150 (121) 



m 

7 



in solar units, where it is assumed that the galaxy extends out 
to the virial radius of ss 290 kpc with the same density profile 
r~ 2 . Evidently, this exceeds the stellar mass-to-light ratio by 
far. Of course, the mass-to-light ratio of galaxies depends on 
the limiting radius assumed. Values of m/l m 30 are often 
quoted, which are typically based on the largest radii out to 
which rotation curves can be measured. 



3. The galaxy population 

Galaxy luminosities are observed to be distributed accord- 



ing to the Schechter function (Schechter 1976) 



dJV 
dL 



U \U 



exp 



L 



(122) 



(115) where the normalising factor is O t « 3.7 x 10 3 Mpc 3 , the 



scale luminosity is L t given under ( 120 ) above, and the power- 
law exponent is a « 1 ( jMadgwick et q/.||2002| ). Irrespective 
of which physical processes this distribution originates from, 
its functional form turns out to characterise mixed galaxy pop- 
ulations very well, even in galaxy clusters. 

The luminosity density in galaxies is easily found to be 



-°e-'dl 



r°° dN f °° , 

L— dL = OhZ, l 1 -" 
Jo dL Jo 

= T(2 - qt)0»L» ~ <D«Z, ~ 7.4 x 10 7 



Mpc 3 



(123) 



The average mass-to-light ratio ( 121 1 then allows converting 
this number into a mass density, 



<Pg> = (7/ A 



1.1 x 10" 



Mpc 3 



7.5 x 10" 3I gcm" 3 



(124) 
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and thus, with p cl0 = 9. 15 x 1(T 30 g cm" 3 from Eq. ([83}, 

Q g0 ~ 0.08 . (125) 

Of course, estimates based on the more conservative mass-to- 
light ratio m/l » 30 find values which are lower by a factor 
of ~ 5. In summary, this shows that the total mass expected 
to be contained in the dark-matter halos hosting galaxies con- 
tributes of order 8% to the critical density in the Universe. 



Hence, the mass required to keep cluster galaxies bound de- 
spite their high velocities exceeds the mass in galaxies by 
about an order of magnitude, even if the entire virial mass of 
the galactic halos is accounted for dZwicky 1933 1937]>. The 



stellar mass apparently contributes just about one per cent to 
the total cluster mass. 



2. Mass in the hot intracluster gas 



B. Mass in galaxy clusters 

1 . Kinematic masses 

The next step upward in the cosmic hierarchy are galaxy 
clusters, which were first identified as significant galaxy over- 



densities in relatively small areas of the sky (Herschel 1786 
|1789[ ). Although the majority of galaxies is not in galaxy 
clusters, rich galaxy clusters contain several hundred galax- 
ies, which by themselves contain a total amount of mass 



M„ < 10 2 L, 



3 x 10 14 M R 



(126) 



The mass in stars is of course considerably lower. With the 



mean stellar mass-to-light ratio of m/l w 6.4 from ( 1 14- 1, the 
same luminosity implies 



M. < 1.3 x 10 13 M G 



(127) 



The stellar mass of the Coma cluster, for instance, is infer red 
to be M t> coma * (1.4 ± 0.3) x 10 13 M (TWhite et aZ.||l993i. 



The galaxies in rich galaxy clusters move with typical ve- 
locities of order < 10 3 kms~', measured from the redshifts in 
their spectra. Therefore, only one component of the galaxy 
velocity is observed. Its distribution is characterised by the 
velocity dispersion tr„. 

If these galaxies were not gravitationally bound to the clus- 
ters, the clusters would disperse within < 1 Gyr. Since they 
exist over cosmological time scales, clusters need to be (at 
least in some sense) gravitationally stable. Assuming virial 
equilibrium, the potential energy of the cluster galaxies should 
be minus two times the kinetic energy, or 



GM 
~R~ 



3crf 



(128) 



where M and R are the total mass and the virial radius of the 
cluster, and the factor three arises because the velocity disper- 
sion represents only one of three velocity components. We 
combine ( |128| l with ( |116| l to find 



9cr 2 



1/2 



\m0nGp) 



>.5Mpc ( 



1000 km s" 



(129) 



and, with ( fT28| , 



Galaxy clusters are diffuse sources of thermal X-ray 
emission. Their X-ray continuum is caused by thermal 
bremsstrahlung, whose bolometric volume emissivity is 



j x = Z 2 g«C x n 2 Vr 



(131) 



in cgs units, where Z is the ion charge, gg is the Gaunt factor, 
n is the ion number density, T is the gas temperature, and 



C x = 2.68 x 10' 24 



(132) 



in cgs units, if T is measured in keV. 

A common simple, axisymmetric model for the gas-density 
distribution in clusters, adapted to X-ray observations, is 



n(x) — 



n 



(1 + x 2 ) 3 ^ 2 



(133) 



where r c and 9 C are the physical and angular core radii. The 
line-of-sight projection of the X-ray emissivity yields the X- 
ray surface brightness as a function of the angular radius 9, 



Sx(fi) 



-f 



jxdz = 



V^r(3/3- 1/2) Z 2 g s C x r c n\ VF 



TQyS) 



[1 +x 2 ] 



2-13/5-1/2 



(134) 

where we have combined < | 1 3 1 J > and < | 1 3 3 1 > and assumed for 
simplicity that the cluster is isothermal, so T does not change 
with radius. The latter equation shows that two parameters 



of the density profile ( 133 1, namely the slope /3 and the (an- 



gular) core radius 6> c , can be read off the observable surface- 
brightness profile. 

The missing normalisation constant can be obtained from 
the X-ray luminosity, 



Lx = 4tt I j x r 2 dr = 4nr 3 c Z 2 g s C x V^ng 
Jo 



V^r(3/S-3/2) 
41X3/3) 



(135) 

and a spectral determination of the temperature T. Finally, the 
total mass of the X-ray gas enclosed in spheres of radius R is 



: f n{ 
Jo 



M X (R) = An n{r)r A dr 
Jo 



(136) 



which is a complicated function for general ji. For /3 = 2/3, 
which is a commonly measured value, 



R R 

M X (R) = 4nr^n \ arctan — 

1 r c Tr 



(137) 



which is of course formally divergent for R — > oo because the 



/ q-^ y wnicn is or course rormany divergent rc 

1.7x10 M b ( 1Q0Qkms -i J • (130) density falls off oc r~ 2 for /? = 2/3 and r 
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Inserting typical numbers, we set Z = 1 = gg , f3 = 2/3 and 
assume a hypothetic cluster with Lx = 10 45 erg s _1 , a temper- 
ature of kT = lOkeV and a core radius of r c = 250 kpc = 



7.75 x 10 23 cm. Then, (I135J yields the central ion density 



no = 5 x 10 3 cm 3 



and thus the central gas mass density 



Po = »ip«o = 8.5 x 10 27 gem 



(138) 



(139) 



Based on the virial radius \\29\ and on the mass ( |137[ ), we 
find the total gas mass 



1.0 x 10 14 M P 



(140) 



This is of the same order as the cluster mass in galaxies, and 
approximately one order of magnitude less than the total clus- 
ter mass. 

It is reasonable to believe that clusters are closed systems 
in the sense that there cannot have been much material ex- 
change between their interior and their surroundings ( |White 
\et al.\ |1993[ ). If this is indeed the case, and the mixture be- 
tween dark matter and baryons in clusters is typical for the 
entire universe, the matter-density parameter should be 



a 



m() 



M 



Ml 



M. + M x 



10C2; 



bO 



0.4 , 



(141) 



given the determination ( 1 10 1 of Qbo from primordial nucle- 
osynthesis. More precise estimates based on detailed investi- 
gations of individual clusters yield 



mO 



0.3 . 



(142) 



3. Alternative cluster mass estimates 



Cluster masses can be estimated in several other ways. One 
of them is directly related to the X-ray emission discussed 
above. The hydrostatic Euler equation for a gas sphere in equi- 
librium with the spherically-symmetric gravitational potential 
of a mass M(r) requires 



1 dp 
p dr 



GM(r) 



(143) 



where p and p are the gas density and pressure, respectively. 
Assuming an ideal gas, the equation of state is p — nkT, where 
n = p/(pm p ) is the particle density and T the temperature, p. 
is the mean particle mass in units of the proton mass m p . A 
mixture of hydrogen and helium with helium fraction Y has 
p = 4/(8 - 5Y), or p * 0.59 for Y = 0.25. If we further 
simplify the problem assuming an isothermal gas distribution, 
we can write 



kT dp GM(r) 



pm p p dr 
or, solving for the mass 

M(r) = - 



rkT d In p 
Gpm p d In r 



(144) 



(145) 



For the /3 model introduced in ( 133 I, the logarithmic density 
slope is 



d In p d In n 
d In r d In r 



= -3/3 



1 + x 2 ' 



(146) 



thus the cluster mass is determined from the slope of the X-ray 
surface brightness and the cluster temperature, 



M(r) 



3/3rkT x 2 
Gpm p 1 + x 2 



(147) 



With the typical numbers used before, i.e. R m 2.5 Mpc, p 1 
2/3 and kT = lOkeV, the X-ray mass estimate gives 



M(R) w 1.8 X 1O 15 M , 



(148) 



in reassuring agreement with the mass estimate ( |130| l from 
galaxy kinematics, and again an order of magnitude more than 
either the gas mass or the stellar mass. 

A third, completely independent way of measuring cluster 
masses is provided by gravitational lensing. Without going 
into any detail, we mention here that it can generate image 
distortions from which the projected lensing mass distribution 
can be reconstructed. Mass estimates obtained in this way 
by and large confirm those from X-ray emission and galaxy 
kinematics, although interesting discrepancies exist in detail 
( Cypriano etai\ |2"004) poekstra| |2"0"07] l. 

None of the cluster mass estimates is particularly reliable 
because they are all to some degree based on stability and 
symmetry assumptions. For mass estimates based on galaxy 
kinematics, for instance, assumptions have to be made on 
the shape of the galaxy orbits, the symmetry of the gravita- 
tional potential and the mechanical equilibrium between or- 
biting galaxies and the body of the cluster. Numerous as- 
sumptions also enter X-ray based mass determinations, such 
as hydrostatic equilibrium, spherical symmetry and, in some 
cases, isothermality of the intracluster gas. Gravitational lens- 
ing does not need any equilibrium assumption, but inferences 
from strongly distorted images depend very sensitively on the 
assumed symmetry of the mass distribution, and measure only 
the mass enclosed by the lensed images. 



C. Mass density from cluster evolution 

A conceptually very interesting constraint on the cosmic 
mass density is based on the evolution of cosmic structures. 
Abell's cluster catalogue ( Ab ell et al.\ |1989~] l covers the red- 
shift range 0.02 < z < 0.2, which encloses a volume of 
=s 9x 10 s Mpc 3 . There are 1894 clusters in that volume, which 
yields an estimate for the spatial cluster density of 



« c ^2x 10~ 6 Mpc~ 3 



(149) 



It is a central assumption in cosmology that structures 
formed from Gaussian random density fluctuations. The 
spherical collapse model then asserts that gravitationally 
bound objects of mass M form where the linear density con- 
trast, smoothed within spheres of radius R enclosing the aver- 
age mass M, exceeds a critical threshold of i5 c ~ 1 .686, quite 
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independent of cosmology. The probability for this to happen 
in a Gaussian random field with a standard deviation cr R (z) is 



Pc(z) 



^ erfc ( r _^ c 

2 \ V2o*(z) 



where 



o-r(z) = o- m) D + (z) 



(150) 



(151) 



because the linear growth of the density contrast is determined 
by the growth factor, for which a fitting formula was given in 

Now, the present-day standard deviation ctrq must be cho- 
sen such as to reproduce the observed number density of clus- 



ters given in ( 149 1. The fraction of cosmic matter contained 
in clusters is approximated by 



Mn c 
p cr ofi, 



m0 



3 x 10- 3 O-{, . 



(152) 



The standard deviation <x in ( |150| > must now be chosen such 
that this number is reproduced, which yields 



O-R0 



0.61 
0.72 



^m0 



1.0 

0.3 



(153) 



Equations ( 150 1 and ( |151| l can now be used to estimate how 
the cluster abundance should change with redshift. Simple 
evaluation reveals that the comoving cluster abundance is ex- 
pected to drop very rapidly with increasing redshift if f2 m o is 
high, and much more slowly if Q m o is low. Qualitatively, this 
behaviour is easily understood. If, in a low-density universe, 
clusters do not form early, they cannot form at all because the 
rapid expansion due to the low matter density prevents them 
from growing late in the cosmic evolution. From the observed 
slow evolution of the cluster population as a whole, it can be 
concluded that the matter density must be low. Estimates ar- 
rive at 

£} ra0 * 0.3 , (154) 
in good agreement with the preceding determinations. 

D. Conclusions 



at, say, up to one minute after the Big Bang. According to 
(90 1, the scale factor at this time was a « 10~ 10 , and thus the 
matter density was of order p m as 10 3 VcrO ~ 10 gem" 3 . The 
horizon size is ru « ct « 1.8 x 10 12 cm, thus the mass en- 
closed by the horizon is as 3 x 10 4 M G , which limits possible 
black-hole masses from above. 

Gravitational microlensing was used to constrain the 
amount of dark, compact objects of subsolar mass in the halo 
of the Milky Way. Although they were found to contribute 
part of the mass, they can certainly not account for all of 
it ( |Alcock et ai\ [2000) |Lasserre et~aT\ [2000]!. Black holes 
with masses of the order 10 b M G should have been found 
through their microlensing effects (Wambsganss and Paczyn 
|ski|[T99 2). The abundance of lower-mass black holes is lim 
ited by stellar-dynamics arguments, in particular by the pres- 
ence of cold disks (Kless eiTand Burke"rt]|1996[|Rix and Lake| 
[T993] l. 

We are thus guided to the conclusion that the dark matter 
is most probably not baryonic and not composed of compact 
dark objects. We shall see in § VLB. 2 that and why the most 
favoured hypothesis now holds that it is composed of weakly 
interacting massive particles. 

Neutrinos, however, are ruled out because their total mass 
has been measured to be way too low (< 2eV 2 ). Their corre- 
spondingly high velocities require objects more massive than 
galaxies to keep them bound. The formation of galaxies would 
be much delayed in this case until larger-scale objects could 
fragment into smaller pieces. This is the opposite of the ob- 
served hierarchy of cosmic structure formation, which clearly 
shows that galaxies appear substantially earlier than galaxy 
clusters. Therefore, the conclusion seems inevitable that the 
dark matter must be cold, i.e. consisting of particles moving 
non-relativistically (Peebles, 1982). A sequence of numerical 
simulations has shown that structure formation in a universe 
filled with cold dark matter can be brought into agreement 



with the observed large-scale cosmic structures (Davis et al. 
[T9851 |White eTaL\\TWl) . 



VI. THE COSMIC MICROWAVE BACKGROUND 



A. The isotropic CMB 

1 . Thermal history of the Universe 



What does it all mean? The preceding discussion should 
have demonstrated that the matter density in the Universe is 
considerably less than its critical value, approximately one 
third of it. Since only a small fraction of this matter is vis- 
ible, we call the invisible majority dark matter. 

What is this dark matter composed of? Can it be baryons? 
Tight limits are set by primordial nucleosynthesis, which pre- 
dicts that the density in baryonic matter should be Q B x 0.04, 
cf. ( |1 10[ ). In the framework of the Friedmann-Lemaitre mod- 
els, the baryon density in the Universe can be higher than this 
only if baryons are locked up in some way before nucleosyn- 
thesis commences. They could form black holes before, but 
their mass is limited by the mass enclosed within the horizon 



How does the Universe evolve thermally? We have seen 
that the abundance of 4 He shows that the Universe must have 
gone through an early phase which was hot enough for the 
nuclear fusion of light elements. But was there thermal equi- 
librium? Thus, can we speak of the "temperature of the Uni- 
verse"? 

From isotropy, we must conclude that the Universe ex- 
panded adiathermally: no net heat can have flowed between 
any two volume elements in the Universe because any flow 



Particle Data Group, http://pdg.lanl.gov/ 
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would have singled out a direction, which is forbidden by 
isotropy. An adiathermal process is adiabatic if it proceeds 
slow enough for equilibrium to be maintained. Then, it is 
also reversible and isentropic. Of course, irreversible pro- 
cesses must have occurred during the evolution of the Uni- 
verse. However, as we shall see later, the entropy in the Uni- 
verse is so overwhelmingly dominated by the photons of the 
microwave background radiation that no entropy production 
by irreversible processes can have added a significant amount 
of entropy. Thus, we assume that the Universe has in fact ex- 
panded adiabatically. 

As the Universe expands and cools, particles are diluted 
and interaction rates drop, so thermal equilibrium must break 
down at some point for any particle species because collisions 
become too rare. Very early in the Universe, however, the 
expansion rate was very high, and it is important to check 
whether thermal equilibrium can have been maintained de- 
spite the rapid expansion. The collision probability between 
any two particle species will be proportional to their number 
densities squared, oc « 2 , because collisions are dominated by 
two-body encounters. The collision rate, i.e. the number of 
collisions experienced by an individual particle with others 
will be oc «, which is oc gt 3 for non-relativistic particles. Thus, 
the collision time scale was r co u oc a 3 . 

According to Friedmann's equation, the expansion rate in 
the very early Universe was determined by the radiation den- 
sity, and thus proportional to oc a/a oc or 2 , and the expansion 
time scale was r exp oc a 2 . Equilibrium could be maintained as 
long as the collision time scale was sufficiently shorter than 
the expansion time scale, 



Trail < T e 



(155) 



which turns out to be well satisfied in the early Universe when 
a <K 1 . Thus, even though the expansion rate was very high 
in the early Universe, the collision rates were even higher, and 
thermal equilibrium could indeed be maintained. 

The final assumption is that the components of the cos- 
mic fluid behave as ideal gases. By definition, this requires 
that their particles interact only via short-ranged forces, which 
implies that partition sums can be written as powers of one- 
particle partition sums and that the internal energy of the fluids 
does not depend on the volume occupied. This is a natural as- 
sumption which holds even for charged particles because they 
shield opposite charges on length scales small compared to 
the size of the observable universe. 

It is thus well justified to assume that there was thermal 
equilibrium between all particle species in the early universe, 
that the constituents of the cosmic "fluid" can be described as 
ideal gases, and that the expansion of the universe can be seen 
as an adiabatic process. In later stages of the cosmic evolu- 
tion, particle species will drop out of equilibrium when their 
interaction rates fall below the expansion rate of the Universe. 
As long as all species in the Universe are kept in thermody- 
namic equilibrium, there is a single temperature characteris- 
ing the cosmic fluid. Once particle species drop out of ther- 
mal equilibrium because their interaction rates decrease, their 
temperatures, if defined, may begin deviating. Even then, we 



characterise the thermal evolution of the Universe by the tem- 
perature of the photon background. 



2. Mean properties of the CMB 

As discussed before, the CMB had been predicted in order 
to expla in the abundance of the light elements, in particular 



of 4 He ( |Alpher and Herman 



1949 



Gamow 



1948). It was 



serendipitously discovered by Penzias and Wilson in 1965 
(Penzias and Wilson 1965 i. Measurements of the energy 



density in this radiation were mostly undertaken at long (ra- 
dio) wavelengths, i.e. in the Rayleigh-Jeans part of the CMB 
spectrum. To firmly establish that the spectrum is indeed the 
Planck spectrum expected for thermal black-body radiation, 
the FIRAS experiment was placed on-board the COBE satel- 
lite, where it measured the best realisation of a Planck spec- 
trum ever observed ( |Fixsen and M ather 2002 ; Mather ~f al.\ 
[T994] l. 

The temperature of the Planck curve best fitting the latest 
measurement of the CMB spectrum by the COBE satellite is 



T = 2.726 K . 



(156) 



which implies the mean number density «cmb = 405 cm 3 
of CMB photons and the energy density mcmb — 4.17 x 
10~ 13 erg cirT 3 in the CMB, equivalent to a mass density of 



Pr,CMB = 4.63 x 10 
eter in radiation is 



■34 , 



gem . Using ( 14 1, the density param- 



t0 



8.51 x 10" 



(157) 



which shows that the scale factor at matter-radiation equality 

was 

Q^n 1 a 

fl - = ^3^r 3 - 0xl ° ■ (158) 

This calculation includes three neutrino species besides the 
photons, which were relativistic at the time of matter-radiation 
equality. 

The number density of baryons in the Universe is approxi- 
mately 



^BOPcr a ~ -, „_7 _3 

«b ~ — ~ 2.3 x 10 cm 



(159) 



confirming that the photon-to-baryon ratio is extremely high, 



405 



2.3 x 10- 7 



1.8 x 10 y . 



(160) 



3. Decoupling of the CMB 

When and how was the CMB set free? While the Universe 
was sufficiently hot to keep electrons and protons separated 
(we neglect heavier elements here for simplicity), the photons 
scattered off the charged particles, their mean free path was 
short, and the photons could not propagate. When the Uni- 
verse cooled below a certain temperature, electrons and pro- 
tons combined to form hydrogen, free charges disappeared, 
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the mean free path became virtually infinite and photons could 
freely propagate. 

The recombination reaction 



p + e <-> H + y 



(161) 



can thermodynamically be described by Saha's equation for 
the ionisation fraction x, 



1 -x 



4 V2^(3)77 \ kT 



2\3/2 



0.26 



kT 



(162) 

where x is the ionisation energy of hydrogen, x - 13.6 eV, 
and £ is the Riemann Zeta function. 

Notice that Saha's equation contains the inverse of the r\ pa- 
rameter ( 160 1, which is a huge number due to the high photon- 
to-baryon ratio in the Universe. This counteracts the expo- 
nential which would otherwise guarantee that recombination 
happens when kT * x> i- e - at T ~ 1.6 x 10 5 K. Recom- 
bination is thus delayed by the high photon number, which 
illustrates that newly formed hydrogen atoms are effectively 
reionised by sufficiently energetic photons until the temper- 
ature has dropped well below the ionsation energy. Putting 
.v 



0.5 in ( 162 1 yields a recombination temperature of 



kT n 



0.3 eV 



3000 K 



(163) 



and thus a recombination redshift of z rec ~ 1100. Since this is 
in the early matter-dominated era, the age of the Universe was 
then 



-r 

Jo 



da 



2a 



3/2 



aH(a) 
374 kyr, 



[VT+a(l -2a) + 2a V2 ] 
(164) 



cf. Eq. ( 15 1, where a := a eq /a n 



0.33. 



Recombination does not proceed instantaneously. The ion- 
isation fraction x drops from 0.9 to 0.1 within a temperature 
range of approximately 200 K, corresponding to a redshift 
range of 



Az 



dz 
dr 



AT 



or a time interval of 



At 



Aa 
~aH 



d / T 



Az 



AT 

AT* — * 75 (165) 



//oVCWl+z) 5 / 2 



50 kyr . 



(166) 



We are thus led to conclude that the CMB was released when 
the Universe was approximately 374,000 years old, during a 
phase that lasted approximately 50,000 years. We have de- 
rived this result merely using the present temperature of the 
CMB, the photon-to-baryon ratio, the Hubble constant and the 
matter density parameter O m o- The cosmological constant or 
a possible curvature of the Universe do not matter here. 



B. Structures in the CMB 



1. The dipole 



The Earth is moving around the Sun, the Sun is orbiting 
around the Galactic centre, the Galaxy is moving within the 
Local Group, which is falling towards the Virgo cluster of 
galaxies. We can thus not expect the Earth to be at rest with 
respect to the CMB. We denote the net velocity of the Earth 
with respect to the CMB rest frame by u e . Lorentz transfor- 
mation shows that, to lowest order in v e /c, the Earth's motion 
imprints a dipolar intensity pattern on the CMB with an am- 
plitude of 



AT 



(167) 



The dipole's amplitude has been measured to be m 1.24raK, 



from which the Earth's velocity is inferred to be (Fixsen et al. 

[T996) 



(371 + 1.5)kms" 



(168) 



This is the highest-order deviation from isotropy in the CMB, 
but it is irrelevant for our purposes since it does not allow any 
conclusions on the Universe at large. 



2. Expected amplitude of CMB fluctuations 

It is reasonable to expect that density fluctuations in the 
CMB should match density fluctuations in the matter because 
photons were tightly coupled to baryons by Thomson scat- 
tering before recombination. For adiabatic fluctuations, the 
density contrast 6 in matter is 3/4 that in radiation. Since the 
radiation density is oc T 4 , a density contrast 6 is thus expected 
to produce relative temperature fluctuations of order 



3 4T 3 6T 

4 T 4 



6T 

Y 



(169) 



Obviously, there are large-scale structures in the Universe to- 
day whose density contrast reaches or even substantially ex- 
ceeds unity. Assuming linear structure growth on scales where 
#0 ~ 1, and knowing the scale factor at recombination, we can 
thus infer that relative temperature fluctuations of order 



6T 
~T 



So 



1 



1 



3 \D + {a KC )J 3a 



10- 



(170) 



should be seen in the CMB, i.e. fluctuations of order mK, 
similar to the dipole. Such fluctuations, however, were not 
observed, although cosmologists kept searching increasingly 



desperately for decades after 1965 (see (Uson and Wilkinson 
1984) for an example). 

Why do they not exist? The estimate above is valid only 
under the assumption that matter and radiation were tightly 
coupled. Should this not have been the case, density fluctua- 
tions did not need to leave such a pronounced imprint on the 
CMB. In order to avoid the tight coupling, the majority of 
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matter must not interact electromagnetically. Thus, we con- 
clude from the absence of mK fluctuations in the CMB that 
matter in the Universe must be dominated by something that 



does not interact with light (IPeebles 1982 1. This is perhaps 



the strongest argument in favour of not-electromagnetically 
interacting dark matter. 



3. Expected CMB fluctuations 

Before we come to the results of CMB observations and 
their significance for cosmology, let us summarise which 
physical effects we expect to imprint structures on the CMB 



(Bond and 


Efstathiou 


1984 


Hu and Sugiyama 


1996 


Pee- 


bles and Yu| 


1970 


Seljak anc 


Zaldarriaga 


1996 


Sunyaev and 



Zeldovich |1970 i. The basic hypothesis is that the cosmic 
structures that we see today formed via gravitational instabil- 
ity from seed fluctuations in the early Universe, whose infla- 
tionary origin is likely, but yet unclear. This implies that there 
had to be density fluctuations at the time when the CMB was 
released. Via Poisson's equation, these density fluctuations 
were related to fluctuations in the Newtonian potential. Pho- 
tons released in a potential fluctuation (50 lost energy if the 
fluctuation was negative, and gained energy when the fluctua- 
tion was positive. This energy change can be translated to the 
temperature change 



ST 

Y 



1 54> 



(171) 



which is called the Sachs-Wolfe effect ( |Sachs and Wolfe] 



1967). The factor 1/3 in front is caused by the gravitational 
redshift being offset by time retardation in the gravitational 
field of the perturbation. 

Let us briefly look into the expected statistics of the Sachs- 
Wolfe effect. We introduced the power spectrum of the density 



fluctuations in ( 27 1 as the variance of the density contrast in 
Fourier space. Poisson's equation implies 



5<b cc - 



k 2 



(172) 



for the Fourier modes of the gravitational potential. Thus 
the power spectrum of the temperature fluctuations due to the 
Sachs- Wolfe effect is 



P T oc Pa, cc 




k <^k 



eq 



(173) 



eq 



according to ( |28] l with n s = 1. This shows that the Sachs- 
Wolfe effect can only be important at small k, i.e. on large 
scales, and dies off quickly towards smaller scales. 

The main constituents of the cosmic fluid were dark mat- 
ter, baryons, electrons and photons. Overdensities in the dark 
matter compressed the fluid due to their gravity until the rising 
pressure in the tightly coupled baryon-electron-photon fluid 
was able to counteract gravity and drive the fluctuations apart. 
In due course, the pressure sank and gravity won again, and 
so forth: the baryon-electron-photon fluid underwent acous- 
tic oscillations, but not the dark matter, which was decoupled. 



Since the pressure was dominated by the photons, whose pres- 
sure is a third of their energy density, a good approximation 
to the sound speed of the tightly coupled photon-baryon fluid 

was 




0.58 c . 



(174) 



Only such density fluctuations could undergo acoustic oscilla- 
tions which were small enough to be crossed by sound waves 
in the available time. The largest comoving length that could 
be travelled by sound waves was the comoving sound horizon 



f 

Jo 



c,dt 




y/a) = 163.3 Mpc 



a Ho V^mO 

(175) 

Larger-scale density fluctuations could not oscillate. We saw 
in ( 18 1 and ( 19 1 that the comoving angular-diameter distance 
from today to scale factor a rec is fK[w( a iec)] = w(a lec ) if we 
assume spatial flatness, K — 0. In excellent approximation, 



f K [w(a lec )] = w(a rec ) = 3.195 



-0.4 



H \ 0.3 

Thus, the sound horizon sets an angular scale of 



w(a rec ) 



= 0.66° . 



(176) 



(177) 



This angular scale can be read off the first acoustic peak in 
the CMB power spectrum (see § VLB. 5 below). Its relation 



to the physical sound horizon ( 175 1 depends almost precisely 
on the sum of Q. m o and Qao if all else remains fixed. Thus, 
the location of the first acoustic peak determines the spatial 
curvature of the cosmological model. When combined with 
measurements of Ho, latest data (cf. Tab. show that K — 
to high precision. 

A third effect influencing structures in the CMB is caused 
by the fact that, as recombination proceeds, the mean-free 
path of the photons increases. If n e = xn B is the electron 
number density and <x T is the Thomson cross section, the co- 
moving mean-free path is 



A « (xn B cr T ) 



(178) 



As the ionisation fraction x drops towards zero, the mean-free 
path approaches infinity. After N scatterings, the photons will 
have diffused by 

A D ~ VAU (179) 
on average. The number of scatterings per unit time is 

dN ~ xn B o- T cdt , (180) 
and thus the diffusion scale is given by 

cdt 



4 * J A 2 dN » J 



X«bO"T 



(181) 
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Inserting x » 1 /2 as a typical value into the latter integrand, 
we can approximate 



2cAt 
n B o"T 



(182) 



Numerical evaluation returns a comoving damping length of 
order 50Mpc, corresponding to angular scales of 6 D * 10' 
on the sky. The effect of this diffusion process is called Silk 
damping (Silk 1968). We thus expect three principal mech- 
anisms to shape the appearance of the microwave sky: the 
Sachs- Wolfe effect on the largest scales, the acoustic oscilla- 
tions on scales smaller than the sound horizon, and Silk damp- 
ing on scales smaller than a few arc minutes. 



4. CMB polarisation 

If the CMB does indeed arise from Thomson scattering, 
interesting effects must occur because the Thomson scatter- 
ing cross section is polarisation sensitive, and can thus con- 
vert unpolarised into linearly polarised radiation. Suppose an 
electron is illuminated by unpolarised radiation from the left, 
then the radiation scattered towards the observer will be lin- 
early polarised in the perpendicular direction. Likewise, un- 
polarised radiation incoming from the top will be linearly po- 
larised horizontally after being scattered towards the observer. 
Thus, if the electron is irradiated by a quadrupolar intensity 
pattern, the scattered radiation will be partially linearly po- 
larised. The polarised intensity is expected to be of order 10% 
of the total intensity. The polarised radiation must reflect the 
same physical effects as the unpolarised radiation, and the two 
must be cross-correlated. Much additional information on the 
physical state of the early Universe is thus contained in the 
polarised component of the CMB, apart from the fact that the 
mere detection of the polarisation supports the physical pic- 
ture of the CMB's origin. 



The shape of the CMB power spectrum reflects the three 
physical mechanisms identified above: at small I (on large 
scales), the Sachs-Wolfe effect causes a feature-less plateau, 
followed by pronounced maxima and minima due to the 
acoustic oscillations, damped on the smallest scales (largest 
I) by Silk damping. 
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FIG. 7 Schematic appearance of the CMB power spectrum due to 
the three dominating physical effects defining its shape: the Sachs- 
Wolfe effect on largest scales, the Silk damping on smallest scales, 
and the acoustic oscillations in between. 

The detailed shape of the CMB power spectrum depends 
sensitively on the cosmological parameters, which can in turn 
be determined by fitting the theoretically expected to the mea- 
sured C/. This is the main reason for the detailed and sensi- 
tive CMB measurements pioneered by COBE, continued by 
ground-based and balloon experiments, and culminating re- 
cently in the spectacular results obtained by the WMAP satel- 
lite. 



5. The CMB power spectrum 



6. Secondary anisotropies and foregrounds 



Fourier transformation is impossible on the sphere, but the 
analysis of the CMB proceeds in a completely analogous 
way by decomposing the relative temperature fluctuations into 
spherical harmonics, finding the spherical-harmonic coeffi- 
cients 



ST 

ai m = I d z — Y, m (6) , 



and from them the power spectrum 



m=—l 



(183) 



(184) 



which is the equivalent on the sphere to the three-dimensional 



power spectra defined in ( 27 1. The average over m expresses 
the expectation of statistical isotropy. If parts of the observed 
sky need to be masked, care has to be taken to orthonormalise 
the spherical harmonics on the remaining domain. 



By definition, the CMB is the oldest visible source of 
photons because all possible earlier sources could not shine 
through the hot cosmic plasma. Therefore, every source that 
produced microwave photons since, or that produced photons 
which became redshifted into the microwave regime by now, 
must appear superposed on the CMB. The CMB is thus cov- 
ered by layers of foreground emission that have to be unveiled 
before the CMB can be observed. 

Broadly, the CMB foregrounds can be grouped into point 
sources and diffuse sources. The most important among the 
point sources are infrared galaxies at high redshift and bodies 
in the Solar System such as the major planets, but even some 
of the asteroids. 

The population of infrared sources at high redshift is poorly 
known, but the angular resolution of CMB measurements has 
so far been too low to be significantly contaminated by them. 
Future CMB observations will have to remove them carefully 
( |Toffolatti etai] [1998) . 
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Microwave radiation from bodies in the Solar System has 
so far been used to calibrate microwave detectors. CMB ob- 
servations at an angular resolution below ~ 10' are expected 
to detect hundreds of minor planets. 

Diffuse CMB foregrounds are mainly caused by our Galaxy 
itself. There are three main components: synchrotron emis- 
sion, emission from warm dust, and bremsstrahlung (Dame 



et al. [2001] |Finkbeiner 2003 Haslam et al. 1982| |Schlegel 



eFoT] [1998) . 

Synchrotron radiation is emitted by relativistic electrons in 
the Galaxy's magnetic field. It is highly linearly polarised and 
has a power-law spectrum falling steeply from radio towards 
microwave frequencies. It is centred on the Galactic plane, 
but shows filamentary extensions from the Galactic centre to- 
wards the Galactic poles. 

The dust in the Milky Way is also concentrated in the 
Galactic plane. It is between (10 . . . 20) K warm and thus sub- 
stantially warmer than the CMB itself. It has a thermal spec- 
trum modified by frequency-dependent self absorption. Due 
to its higher temperature, the warm dust has a spectrum rising 
with increasing frequency in the frequency window where the 
CMB is usually observed. 

Bremsstrahlung radiation is emitted by ionised hydrogen 
clouds (HII regions) in the Galactic plane. Its spectrum has 
the shape typical for thermal free-free radiation, falling expo- 
nentially at photon energies near and above the gas tempera- 
ture, but is flat at CMB frequencies. Further sources of mi- 
crowave radiation in the Galaxy are less prominent. Among 
them are line emission from CO molecules embedded in cool 
gas clouds. 

The falling synchrotron spectrum, the flat spectrum of 
the free-free radiation, and the rising spectrum of the warm 
dust create a window for CMB observations between ~ 
(100 . . . 200) GHz. The different spectra of the foregrounds, 
and their non-Planckian character, are crucial for their proper 
removal from the CMB data. Therefore, CMB measurements 
have to be carried out in multiple frequency bands. 

Secondary anisotropics are caused by propagation effects 
rather than photon emission. The most important are the 
integrated Sachs-Wolfe effect and the (thermal and kinetic) 
Sunyaev-Zel'dovich effects in galaxy clusters. 

The integrated Sachs-Wolfe effect is caused by gravitational 
potential wells deepening while crossed by CMB photons. It 
is determined by the line-of-sight integral of the time deriva- 
tive of the potential fluctuations caused by the density fluctu- 
ations between us and the CMB. By cross-correlating CMB 
temperature fluctuations with structures in the galaxy distri- 
bution, the integrated Sachs-Wolfe effect has indeed been de- 



tected (Gaztanaga et al. 2006 Giannantonio et al. 2006) . It 



constrains the ratio D + (a)/a. 



The Sunyaev-Zel'dovich effect (Sunyaev and Zeldovich 



1972) is due to inverse Compton scattering of CMB photons 
off hot electrons in the intracluster plasma. On average, the 
photons gain energy and are thus moved from the low- to the 
high-frequency part of the spectrum. When observed through 
a galaxy cluster, the CMB therefore appears fainter at low and 
brighter at high frequencies, with the transition at 217 GHz. 
This gives galaxy clusters a peculiar spectral signature in the 



microwave regime as they cast shadows below, emit above, 
and vanish at 217 GHz. Once the angular resolution of CMB 
detectors will drop towards a few arc minutes, a large number 
of galaxy clusters are expected to show up in this way. Besides 
this thermal Sunyaev-Zel'dovich effect, there is the kinetic ef- 
fect caused by the bulk motion of the cluster as a whole, which 
causes CMB radiation to be scattered by the electrons moving 
with the cluster. Very few clusters have so far been detected 
in CMB data, but thousands are expected to be found in future 
missions. 



7. Measurements of the CMB 



Wien's law ( 113 i shows that the CMB spectrum peaks at 
^max * 0.11 cm, or at a frequency of v max m 282 GHz. As we 
saw, Silk damping sets in below ~ 10' arc minutes, thus most 
of the structures in the CMB are resolvable for rather small 
telescopes. According to the formula 



A8x 1.44- 



D 



(185) 



relating the diffraction-limited angular resolution A9 to the ra- 
tio between wavelength and telescope diameter D, we find that 
mirrors with 



D < 1.44: 



55 cm 



(186) 



are sufficient to achieve sufficient angular resolution up to the 
Silk-damping scale 6> D . Detectors are needed which are sensi- 
tive to millimetre and sub-mm radiation and reach jjK sensi- 
tivity, while the telescope optics can be kept rather small and 
simple. 

Two types of detector are commonly used. The first are 
bolometers, which measure the energy of the absorbed ra- 
diation by the temperature increase it causes. They have 
to be cooled to very low temperatures, typically in the ;«K 
regime. The second are so-called high electron mobility de- 
tectors (HEMTs), in which the currents caused by the incom- 
ing electromagnetic field are measured directly. They measure 
amplitude and phase of the waves and are thus polarisation- 
sensitive by construction, which bolometers are not. Polarisa- 
tion measurements with bolometers are possible with suitably 
shaped so-called/eet/ horns guiding the radiation into the de- 
tectors. 

Since water vapour in the atmosphere both absorbs and 
emits microwave radiation through molecular lines, CMB ob- 
servations need to be carried out either at high, dry and cold 
sites on the ground (e.g. in the Chilean Andes or at the South 
Pole), or from balloons rising above the troposphere, or from 
satellites in space. 

After the breakthrough achieved with COBE (Boggess 
et al. 1992), progress was made with balloon experiments 



and with ground-based interferometers. The balloons covered 
a small fraction of the sky (typically ~ 1%) at frequencies 
between 90 and 400 GHz, while the interferometers observe 
even smaller fields at somewhat lower frequencies (typically 
around 30 GHz). 
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The first discovery of the CMB polarisation and its cross- 
correlation with the CMB temperature was achieved with the 
DASI interferometer ( |Kovac et a/.||2002] l. The existence, lo- 
cation and height of the first acoustic peak had been firmly 



established (Wang et al. 20031 before the NASA satellite 



Wilkinson Microwave Anisotropy Probe (WMAP for short) 
was launched, but the increased sensitivity and the full-sky 



coverage of WMAP produced breath-taking results ( Spergel 



et al 



2007 ). WMAP is still operating, measuring the CMB 



temperature at frequencies between 23 and 94 GHz with an 
angular resolution of > 15'. The sensitivity of WMAP is 
barely high enough for polarisation measurements. 

By now, data from the first five years of operation have been 
published, and cosmological parameters have been obtained 
fitting theoretically expected to the measured temperature- 
fluctuation power spectrum and the temperature-polarisation 
power spectrum. Results were given in Tab. [I] adapted from 
dKomatsu etaT\ [2008] l. 
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FIG. 8 Power spectrum of CMB temperature fluctuations as mea- 
sured from the 5-year data of WMAP and several additional ground- 
based experiments. (Courtesy of the WMAP team) 

Although the CMB provides insight into the cosmological 
parameters on its own, it is most powerful when combined 
with other data sets. In particular, the Hubble constant is not 
an independent measurement from the CMB alone. Only by 
assuming a flat ACDM universe, it can be inferred from the 
location of the first acoustic peak in the CMB power spectrum 
to be Ho - 70.1 ± 1.3kms _1 Mpc~', which agrees perfectly 
with the results of the Hubble Key Project and at least one 
interpretation of the gravitational-lens time delays. 

A European CMB satellite mission is under way: ESA's 
Planck satellite has been launched in May 2009. It will ob- 
serve the microwave sky in nine frequency bands between 
30 and 857 GHz with about ten times higher sensitivity than 
WMAP, and an angular resolution of > 5' ( |Tauber|[20oi) . Its 
wide frequency coverage will be very important for substan- 
tially improved foreground subtraction. It will also have suf- 
ficient sensitivity to precisely measure the CMB polarisation 



in some of its frequency bands ( Lamarre et al. 2003 Valen 
2007]!. Moreover, it is expected that Planck will 



ziano et al. 



Sunyaev-Zel'dovich effect. 

VII. COSMIC STRUCTURES 
A. Quantifying structures 

1. Introduction 

We have seen before that there is a very specific prediction 
for the power spectrum of cold dark matter density fluctua- 
tions in the Universe, characterised by ( 28 1. Recall that its 



shape was inferred from the simple assumption that the rms 
mass of density fluctuations entering the horizon should be in- 
dependent of the time when they enter the horizon, and from 
the fact that perturbation modes entering during the radiation 
era are suppressed until matter begins dominating. Given the 
simplicity of the argument, and the corresponding strength of 
the prediction, it is very important for cosmology to find out 
whether the actual power spectrum of matter density fluctu- 
ations does in fact have the expected shape, and if it has, to 
determine the only remaining parameter, i.e. the normalisa- 
tion of the power spectrum. 

Since the comoving wave number k eq of the peak in the 
power spectrum is set by the comoving horizon radius at 
matter-radiation equality d26l, 



k e q 



2n 



i-'eq 



'eq 



1/2 



On 



(187) 



it is proportional to the matter-density parameter Q mU . A mea- 
surement of k eq would thus provide an independent and very 
elegant determination of Q m rj. 

Since the power spectrum is defined in Fourier space, it is 
not obvious how it should be measured. The correlation func- 
tion is related to the power spectrum by 



1 f°° sinfcx , 
f(x) = —= P{k)——tfdk, 
2n l Jo kx 



whose inverse transform is 



P(k)=An f av) 

Jo kx 



sinkx 



x Ax . 



(188) 



(189) 



The power spectrum can thus in principle be determined from 
a measured correlation function ij(x). 



2. Measuring the correlation function 

The correlation function of the three-dimensional density 
field cannot be measured directly. Gravitational lensing by 
large-scale structures comes close by measuring the weighed 
line-of-sight projection of the matter power spectrum; see 



detect a large number of galaxy clusters through their thermal 



§ VIII below. Starting from the assumption that galaxies may 
be tracers of the underlying density field, we can use their cor- 
relation function as an estimate for that of the matter. 
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Suppose we divide space into cells of volume dV small 
enough to contain at most a single galaxy. Then, the prob- 
ability of finding one galaxy in <W\ and another in dV 2 is 



dP = <«(x 1 )«(f2)>dyidy 2 



(190) 



where n is the number density of the galaxies as a function of 
position. If we introduce a density contrast for the galaxies in 
analogy to the density contrast for the matter, 



6n = - - 1 , 
n 



and assume for now that 6n - 5, we find from \\9Q\ with 
n = n(l + 6) 

dP = n 2 <(l+(Ji)(l+5 2 )>dy I dy 2 = n 2 [l+£(x)]dVidV 2 , (192) 

where x is the comoving distance between the two volume 
elements. This shows that the correlation function quantifies 
the excess probability above random for finding galaxy pairs 
at a given distance. 

Thus, the correlation function can be measured by count- 
ing galaxy pairs and comparing the result to the Poisson ex- 
pectation, i.e. to the pair counts expected in a random point 
distribution. Symbolically, 



1+f 



(DP) 
(RR) 



(193) 



where D and R represent the data and the random point set, re- 
spectively. Several other estimators for ^ have been proposed 
which are all equivalent in the ideal situation of an infinitely 
extended point distribution. For finite point sets, their noise 
properties differ (Hamilton, 1993b, Landy and Szalay 1993] >. 
The recipe for measuring f (x) is thus to count pairs separated 
by x in the data D and in the random point set R, or between 
the data and the random point set, and to use one of the esti- 
mators proposed. 

The obvious question is then how accurately £ can be de- 
termined. The simple expectation in the absence of clustering 
is 



= . 



if) 



1 



(194) 



where N p is the number of pairs found. Thus, the Poisson 
error on the correlation function is 

A£ 1 



(195) 



This is a lower limit to the actual error, however, because 
the galaxies are in fact correlated. It turns out that the re- 
sult ( |195[ ) should be multiplied with 1 + A-nhJo,, where J$ is 
the volume integral over £ within the galaxy-survey volume 
( |Peebles||1973| l. The true error bars on £ are therefore hard to 
estimate. 

Having measured the correlation function, it would in prin- 



ciple suffice to carry out the Fourier transform (189 1 to find 
P(k), but this is difficult in reality because of the inevitable 



sample limitations. Consider ( 188 i and an underlying power 
spectrum of CDM shape, falling off oc k~ 3 for large k, i.e. on 



small scales. For fixed x, the integrand in ( 188 i falls off very 
slowly, which means that a considerable amount of small- 
scale power is mixed into the correlation function. Since £ 
at large x is small and most affected by measurement errors, 
any uncertainty in the large-scale correlation function is prop- 
agated to the power spectrum even on small scales. 

A further problem is the uncertainty in the me an ga laxy 



number density n. Since 1 + £ oc n according to ( 192 1, the 



(191) uncertainty in £ due to an uncertainty in n is 



A? 



An 



(196) 



showing that £ cannot be measured with an accuracy better 
than the relative accuracy of the mean galaxy density. 



3. Measuring the power spectrum 

Given these problems with real data, it seems appropri- 
ate to estimate the power spectrum directly. The function to 
be transformed is the density field sampled by the galaxies, 
which can be represented by a sum of Dirac delta functions 
centred on the locations of the N galaxies, 



n(x) = ^5 D (x- Xi) 



(197) 



The Fourier transform of the galaxy density contrast is then 



SJk) 



1 N 

/=1 



(198) 



In the absence of correlations, the Fourier phases of the indi- 
vidual terms are independent, and the variance of the Fourier 
amplitude for a single mode becomes 



-it.?, 



i 

N 



(199) 



This is the so-called shot noise present in the power spec- 
trum due to the discrete sampling of the density field. The 
shot-noise contribution needs to be subtracted from the power 
spectrum of the real, correlated galaxy distribution, 



P(k) 



- y. \w\ 2 



1 

N 



(200) 



where the sum extends over all m Fourier modes with wave 
number k contained in the survey. 

This is not the final result yet, because any real survey typ- 
ically covers an irregularly shaped volume from which parts 
need to be excised because they are overshone by stars or un- 
usable for any other reason. The combined effect of mask and 
irregular survey volume is described by a window function 
f(x) which multiplies the galaxy density, 

n(x) ^ f(x)n(x) , (1 + 6 g ) ^ f(x)(l + 6 g ) , (201) 
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implying that the Fourier transform of the mask needs to be 
subtracted. 

Moreover, the Fourier convolution theorem says that the 
Fourier transform of the product f(x)6 g (x) is the convolution 
of the Fourier transforms f(k) and 5 g (k), 



fS = f* ff, 



f(k')6 g (k' - k)d 3 k' 



(202) 



If the Fourier phases of / and 6 g are uncorrected, which is 
the case if the dimensions of the survey are large enough com- 
pared to the size 2n/k of the density mode, this translates to a 
convolution of the power spectrum, 



^obs — Ml 



\f(k)\ 2 



(203) 



This convolution smoothes the observed compared to the true 
power spectrum and changes its amplitude. 

If the Poisson error dominates in the survey, the different 
modes 6 g (k) can be shown to be uncorrected, and the standard 
deviation after summing over the m modes with wave number 
k is ^jTm/N, which yields the minimal error bar to be attached 
to the power spectrum. 

Thus, the shot noise contribution and the Fourier transform 
of the window function need to be subtracted, and the window 
function needs to be deconvolved, and the amplitude needs to 
be corrected for the effective volume covered by the window 
function before the measured power spectrum can be com- 
pared to the theoretical expectation. 



4. Biasing 

What we have determined so far is the power spectrum of 
the galaxy number- density contrast 6n rather than that of the 
matter density contrast 6. Both are related by a possibly scale- 
dependent bias factor b(k), such that 



5n(k) = b(k)5Jk) . 



(204) 



Clearly, different types of objects sample the underlying mat- 
ter density field in different ways. Galaxy clusters, for in- 
stance, are much more rare than galaxies and are thus expected 
to have a substantially higher bias factor than galaxies. Obvi- 
ously, the bias factor enters the power spectrum squared, e.g. 



P gal = b 2 eal (k)P(k). 



(205) 



It constitutes a major and possibly systematic uncertainty in 
the determination of the matter power spectrum from the 
galaxy power spectrum. 



5. Redshift-space distortions 

For any quantification of three-dimensional structures in 
galaxy surveys, the three-dimensional positions x, of the 
galaxies in the survey need to be known. Distances can be 
inferred only from the galaxy redshifts and thus from galaxy 



velocities. These, however, are composed of the Hubble ve- 
locities, from which the distances can be determined, and the 
peculiar velocities, 



V = ^Hubble + fpec 



(206) 



which are caused by local density perturbations and are unre- 
lated to the galaxy densities. Since observations of individual 
galaxies do not allow any separation between the two velocity 
components, distances are inferred from the total velocity u 
rather than the Hubble velocity as they should be, 



D = 



^Hubble + fpec 



= Aiue + AD , 



(207) 



giving rise to a distance error 6D = v vsc /Hq, the so-called 
redshift-space distortion. 

The redshift-space distortions create a peculiar pattern 



through which they can be corrected (Hamilton 1993a 



|Kaiser|[l987[ ). Consider a matter overdensity such as a galaxy 
cluster, containing galaxies moving with random virial veloc- 
ities in it. The virial velocities of order 1000km s _1 scatter 
around the systemic cluster velocity and thus broaden the red- 
shift distribution of the cluster galaxies. In redshift space, 
therefore, the cluster appears stretched along the line-of-sight, 
which is called the finger-of-god effect. 

In addition, the cluster is surrounded by an infall region 
where the galaxies are not virialised yet, but move in an or- 
dered, radial pattern towards the cluster. Galaxies in front of 
the cluster thus have higher, and galaxies behind the cluster 
have lower recession velocities compared to the Hubble ve- 
locity, leading to a flattening of the infall region in redshift 
space. 

A detailed analysis shows that the redshift-space power 
spectrum P z is related to the real-space power spectrum P by 



P z (k) = P(k)(l+fr 2 f 



(208) 



where fi is the direction cosine between the line-of-sight and 
the wave vector k, and [5 is related to the bias parameter b 
through 



mm) 



(209) 



and /(O m ) is the logarithmic derivative of the growth factor 
D + (a), 



dlnD + (a) 
din a 



0.6 



(210) 



Thus, the characteristic pattern of the redshift-space distor- 
tions around overdensities allows a measurement of the bias 
factor (Hawkins et al. 2003) >. Another way of measuring b is 
based upon gravitational lensing ( Hoekstra et al. 2002a Si- 
|mon et aT\ |2007| ). Measurements of b show that it is in fact 
only weakly scale-dependent, near unity for "ordinary" galax- 



ies, but depends mildly on galaxy luminosity and type ( |Cole 
eFaL\ |2005| |Lahav et al. \ \HW\ |Norberg et al. | [200T 
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6. Baryonic acoustic oscillations 

As we have seen in the discussion of the CMB, the cos- 
mic fluid underwent acoustic oscillations on comoving scales 
smaller than the sound horizon ( |175| l w s - 163.3 Mpc, cor- 
responding to a comoving wave number k s = 2n/w s = 
0.038 Mpc 1 . When the CMB decoupled, the oscillations 
ceased, leaving structures in the cosmic matter distribution 
with a fundamental wavelength of w s and its overtones. The 
other scale characterising the cosmic structures is the horizon 



radius at matter-radiation equality (26 1, which was responsi- 
ble to set the peak location of the matter-fluctuation power 
spectrum at k eq = 0.01 Mpc -1 . 

Thus, at wave numbers as 3.8 times the peak scale and 
above, we expect the wave-like imprint of these baryonic 
acoustic oscillations (BAOs) on top of the otherwise smooth 
dark-matter power spectrum. Mode coupling due to non- 
linear evolution of cosmic structures must have stretched and 
distorted this pattern to some degree. As we shall see below, 
the BAOs have indeed been discovered in the largest galaxy 
surveys. They play an important role in attempts to recover the 
history of the cosmic expansion rate (see (Eisenstein, 2005 \ 
for a recent review). 



B. Measurements and results 

1 . The power spectrum 

Spectacularly successful measurements of the galaxy 
power spectrum became recently possible with the two largest 
galaxy surveys to date, the Two-Degree Field Galaxy Red- 
shift Survey (2dFG RS, (|Colless| [T999) ) and the Sloan Digital 
Sky Survey (SDSS, ( |York et a/.]|2000) ). As anticipated in the 
preceding discussion, an enormous effort had to be made to 
identify galaxies, measure their redshifts, select homogeneous 
galaxy subsamples by luminosity and colour as a function of 
redshift so as not to compare and correlate apples with or- 
anges, estimate the window function of the survey, determine 
the average galaxy number density, correct for the convolution 
with the window function and for the bias, and so forth. 

Moreover, numerical calibration experiments were carried 
out in which all measurement and correction techniques were 
applied to simulated data in the same way as to the real data 
to determine reliable error estimates and to test whether the 
full sequence of analysis steps ultimately yields an unbiased 
result. 

Based on 221,414 galaxies, the 2dFGRS consortium de- 

2005) . 



rived a power spectrum of superb quality ( Cole et al. 



First and foremost, it agrees well with the power-spectrum 
shape expected for cold dark matter on the small-scale side of 
the peak. This is a highly remarkable result on its own. The 
2dFGRS power spectrum also clearly shows a tum-over to- 
wards larger scales, signalling the presence of the peak. The 
survey is still not quite large enough to clearly show the peak, 
but the peak location can be estimated from the flattening of 
the power spectrum. Its proportionality to £2 m0 allows an in- 
dependent determination of the matter density parameter. Fi- 



nally, and most spectacularly, the power spectrum shows the 
baryonic acoustic oscillations, whose amplitude allows an in- 
dependent determination of the ratio between the density pa- 
rameters of baryons and dark matter. 

Apart from the fact that the galaxy power spectrum supports 
the CDM shape on small scales, the results obtained from the 
2dFGRS can be summarised as follows: 



Q m0 = 0.233 ±0.022 
O b o/O m o = 0.185 ±0.046 . 



(211) 



A Hubble constant of h — 0.72 is assumed here. Indirectly, 
the baryon density is constrained to be Qbo ~ 0.04, which is 
in perfect agreement with the value derived from primordial 
nucleosynthesis and the measured abundances of the light el- 
ements. 

Based on 522,280 galaxies, the power spectrum inferred 
from the SDSS (Percival et al. 2Q07\ is also well compatible 
with the CDM shape. The estimate for the matter density pa- 
rameter overlaps with that from the 2dFGRS on large scales, 
^mO = 0.22 ± 0.04, but increases when small scales are in- 
cluded. 
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FIG. 9 Galaxy power spectra obtained from the 2dFGRS ( Cole et al. 
2005) and two releases of the SDSS (|Percival et a/.| |2007[|Tegmark| 



et al. 



|2004| >. The expected shape of the CDM power spectrum is 



well reproduced, and the difference in amplitudes may be attributed 
to scale-dependent galaxy biasing, (from ( |Percival et a/.[|2007] >) 

It should be kept in mind, however, that the assumption of 
linear biasing may turn out to be too simplistic for precise 
cosmological inferences to be drawn from galaxy power spec- 
tra. Nonetheless, power spectra for different types of galaxy 
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will provide invaluable information on the physics of galaxy 
formation. 



viously, the complete weight function 

(w — w')w' 
W(w',w) = - — 



(214) 



VIII. COSMOLOGICAL WEAK LENSING 



A. Cosmological light deflection 

1 . Deflection angle, convergence and shear 

Gravitational lensing was mentioned two times before: first 



in § III.B as a means for measuring the Hubble constant 
through the time delay caused by gravitational light deflec- 
tion, and second as a means for measuring cluster masses 
in § V.B.3 For cosmology as a whole, gravitational lensing 



has also developed into an increasingly important tool (see 



dBartelmann and Schneider] |200T] |Mellier| [19991 jRefregier 
|2003D for reviews). 

Matter inhomogeneities deflect light. Working out this ef- 
fect in the limit of a small Newtonian gravitational potential, 
<1> <k c 2 , leads to the deflection angle 

S{9) = \\ dw' /g( ^~ W) V ± 0[/^(»/)e] ■ (212) 
c z Jo Jk(w) 

It is determined by a weighted integral over the gradient of 
the Newtonian gravitational potential <S) perpendicular to the 
line-of-sight into direction 9 on the observer's sky, and the 
weight is given by the comoving angular-diameter distance 
/k(iv) defined in Q. The integral extends along the comoving 
radial distance w' along the line-of-sight to the distance w of 
the source. 

Equation (212 1 can be intuitively understood. Light is de- 



flected due to the pull of the dimension-less Newtonian grav- 
itational field V ± 0/c 2 perpendicular to the otherwise unper- 
turbed line-of-sight, and the effect is weighed by the ratio be- 
tween the angular-diameter distances from the deflecting po- 
tential to the source, fx(u) - w'), and from the observer to the 
source, /k(w). Thus, a lensing mass distribution very close 
to the observer gives rise to a large deflection, while a lens 
near the source, w' m w, has very little effect. The factor of 
two is a relic from general relativity and is due to space-time 
curvature, which is missing in Newtonian gravity. Assuming 
K = 0, we can replace fx(w) by w. 

It is important to realise that the deflection itself is not ob- 
servable. If all light rays emerging from a source were de- 
flected by the same angle on their way to the observer, no no- 
ticeable effect remained. What is important, therefore, is not 
the deflection angle itself, but its change from one light ray to 
the next. This is quantified by the derivative of the deflection 
angle with respect to the direction 9, 



dag 

89: 



, (w - w')w' c? 2 <t> ,^ 

dw' - '— irir (u/e) . 

w oxjOXj 



(213) 



The additional factor w' in the weight function arises because 
the derivative of the potential is taken with respect to comov- 
ing coordinates x, rather than the angular components Ob- 



vanishes at the observer, w' = 0, and at the source, w' = w, 
and peaks approximately half-way in between. 

For applications of gravitational lensing, it is important to 



distinguish the trace-free part of the matrix (213 i from its 
trace, 



tr y^j„ dwmw ' 



w) 



( «9 2 <D 3 2 <D 
8x 2 <9x? 



(w'9) . (215) 



"1 " A 2 > 

The derivatives of <S> can be combined to the two-dimensional 
Laplacian, which can be replaced by the three-dimensional 
Laplacian because the derivatives parallel to the line-of-sight 
do not contribute to the integral ( |215| l. Thus, we find 



tr 



da/ 
89, 



I dw 

Jo 



W(w',w) AO . 



(216) 



Next, we can use Poisson's equation to replace the Lapla- 
cian of O by the density. In fact, we have to take into account 
that light deflection is caused by density perturbations, and 
that we need the Laplacian in terms of comoving rather than 
physical coordinates. Thus, 

a- 2 A<£> = 4nGp6 , (217) 

where 5 is the density contrast and 

3i/ 2 

P = Poa~ 3 = p cl -Q m0 fir 3 = —^Q. m0 a' 2 (218) 

OKLr 

is the mean matter density. Thus, Poisson's equation reads 

AO = \Hl^l m o S - , (219) 

and ( |216| l becomes 

3// ( 2 O m o 



\39j 



Jr>W 
dw' 
o 



W(w',w)-=2k, (220) 
a 



where we have introduced the (effective) convergence k. 
The trace-free part of the matrix ( |213| l is 



da-, 1 / da>j 



- 6ijK = 



n 72 

72 -71 



(221) 



89 j 2 1 \d9j) 89 j 
which defines the so-called shear components y,. Specifically, 

'3 2 <D 8 2 $>\ 



72 



i r 

= — dw'W(w',w) 
c " Jo 

dw W(w', w) 



<9x 2 8x\ t 
<9 2 <t> 
8x\8x2 



(222) 



Combining results, we can write the matrix of deflection- 
angle derivatives as 



daii 



k + 7\ 72 

72 K-J\ 



(223) 
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This matrix contains the important information on how an im- 
age is magnified and distorted. In the limit of weak gravita- 
tional lensing, the size of a lensed image is changed by the 
relative magnification 



dpi — 2k . 



(224) 



while the image distortion is given by the shear components. 
In fact, an originally circular source with radius r will appear 
as an ellipse with major and minor axes 

r r 

(225) 



k — y 



1 — K + y 



where y = (y\ + yi) ' 2 . The ellipticity of the observed image 
of a circular source thus provides an estimate for the shear, 



e = 



a — b 



a + b 1 — k 



y 



(226) 



2. Power spectra 

Of course, the exact light deflection expected along a par- 
ticular line-of-sight cannot be predicted because the mass dis- 
tribution along that light path is unknown. However, we can 
predict the statistical properties of weak lensing from those of 
the density-perturbation field. We are thus led to the following 
problem: Suppose the power spectrum P(k) of a Gaussian ran- 
dom density-perturbation field 5 is known, what is the power 
spectrum of any weighed projection of 5 along the line-of- 
sight? The answer is given by Limber's equation. Suppose 
the weight function is q(w) and the projection is 

g0) = I dw q(w')6(w'&) . (227) 
Jo 

If q(w) is smooth compared to 5, i.e. if the weight function 
changes on scales much larger than typical scales in the den- 
sity contrast, then the power spectrum of g is 



-f 

Jo 



w' z \ w' 



(228) 



where I is a two-dimensional wave vector which is the Fourier 
conjugate variable to the two-dimensional position 3 on the 
sky. 

Strictly speaking, Fourier transforms are inappropriate be- 
cause the sky is not an infinite, two-dimensional plane. In- 
stead, the appropriate set of orthonormal base functions are 
the spherical harmonics. However, lensing effects are usually 
observed in areas whose solid angle is very small compared 
to the full sky. If this is so, the survey area can be approx- 
imated by a section of the plane locally tangent to the sky, 
and Fourier transforms can be used in this so-called flat-sky 
approximation. 



Equation (220 1 is clearly of the form (227 1 with the weight 
function 



3 H% W(w',w) 

q(w ) = -O m0 ^ 

2 c L a 



(229) 



thus the power spectrum of the convergence is, according to 
Limber's equation, 



PM = 



9QL Hi 



m() 



Jo 



4 c 

with a new weight function 

W(w',w) = 



W 2 (w', w) P 



I 



Mw')j ' 



WW, w) 
aw' 



(230) 



(231) 



While it requires huge data sets and extremely careful data 
analysis to observe the differential magnification bp or the 



convergence k ( |Scranton et al. 2005| l, image distortions can 
in principle be measured in a more straightforward way. With 
a brief excursion through Fourier space, it can easily be shown 
that the power spectrum of the shear is exactly identical to that 
of the convergence, 



PyiD = PM 



(232) 



Thus, the statistics of the image distortions caused by cos- 
mological weak lensing contains integral information on the 
power spectrum of the matter fluctuations. 

Since the shear is defined on the two-sphere (the observer's 
sky), its power spectrum is related to its correlation function 
% y through the two-dimensional Fourier transform 



r a z , 

where J v is the ordinary Bessel function of order v. 



„ 



Idl 

Ik 



P 7 {T)Ml<t>) , (233) 



3. Correlation functions 

In principle, shear correlation functions are measured by 
comparing the ellipticity of one galaxy with the ellipticity of 
other galaxies at an angular distance <p from the first. Ellip- 
ticities are oriented, of course, and one has to specify against 
what other direction the direction of, say, the major axis of a 
given ellipse is to be compared with. Since correlation func- 
tions are measured comparing the shear of galaxy pairs, a 
preferred direction is defined by the line connecting the two 
galaxies of the pair under consideration. 

Let a be the angle between this direction and the major axis 
of the ellipse, then the tangential and cross components of the 
shear are defined by 



y + = y cos 2a , y x = y sin 2a 



(234) 



The factor two is important because it accounts for the fact 
that an ellipse is mapped onto itself when rotated by an angle 
7r. This illustrates that the shear is a spin-2 field: It returns into 
its original orientation when rotated by n rather than 2n. 

The correlation functions of the tangential and cross com- 
ponents of the shear are 



£ ++ (0) = (y + (S)y + (e + . 



1 r 

2 Jo 



Idl 
27 



(235) 
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and 



fxxW = <rx(0)rx(0- 



1 r 

2 Jo 



id/ 

2^ 



P,(/)[Jo(/0)-J 4 W)] , 
(236) 

while the cross-correlation between the tangential and cross 
components must vanish, £ +x (0) = 0. This suggests to define 
the correlation functions f ± = ± £ xx , which are related to 
the power spectrum through £ + = £ r and 



f 

Jo 



id/, 
2^ 



P K (/)J 4 (/(*) 



(237) 



The principle of all measures for cosmic shear is the same: 
They are integrals of the weak-lensing power spectrum times 
filter functions describing the specific response of the mea- 
surement to the underlying power spectrum of density fluc- 
tuations. The width of the filter functions controls the range 
of density-perturbation modes k contributing to one specific 
mode I of weak-lensing on the sky. 

We can now estimate typical numbers for the cosmological 
weak-lensing effect. The power A* in the weak-lensing quan- 
tities such as the cosmic shear is given by the power spectrum 



P K (l) found in (230 1, times the volume in /-space, 



l 2 P K (l) 



(238) 



Assuming a cosmological model with Q m o = 0.3 and £2 A o 
0.7. the CDM power spectrum and a reasonable source red- 

1 /2 

shift distribution, A K ' (I) is found to peak on scales I cor- 
responding to angular scales 27r/Z of 2' ...3', and the peak 
reaches values of 0.04 . . . 0.05. This shows that cosmological 
weak lensing will typically cause source ellipticities of a few 
per cent, and correlations have a typical angular scale of a few 
arc minutes. Details depend on the measure chosen through 
the filter function. 



B. Cosmic-shear measurements 



1 . Typical scales and requirements 
How can cosmic gravitational lensing effects be measured? 



As shown in ( 226 1, the ellipticity of a hypothetic circular 
source is an unbiased estimator for the shear. However, typi- 
cal sources are not circular, but to first approximation ellipti- 
cal themselves. Thus, measuring their ellipticities yields their 
intrinsic ellipticities in the first place. 

Let e (s) be the intrinsic source ellipticity. It is a two- 
component, spin-2 quantity. The cosmic shear adds to that 
ellipticity, such that the observed ellipticity is 



e (s) + y 



(239) 



in the weak-lensing approximation. What is observed is there- 
fore the sum of the signal, y, and the intrinsic noise component 



On sufficiently deep observations, < 20 galaxies per square 
arc minute are routinely detected. Since the full moon has 
half a degree diameter, it covers a solid angle of I5 2 n = 700 



square arc minutes, or < 14, 000 of such distant, faint galax- 
ies! From this point of view, the sky is covered by densely 
patterned "wall paper" of distant galaxies. Thus, it is possi- 
ble to average observed galaxy ellipticities. Assuming their 
shapes are intrinsically independent, the intrinsic ellipticities 
will average out, and the shear will remain, 



<e> * <6 (s) ) + <r> * <r> ■ 



(240) 



It is a fortunate coincidence that the typical angular scale of 
cosmic lensing, which we found to be of order a few arc min- 
utes, is large compared to the mean distance between back- 
ground galaxies, which is of order Vl/20 « 0.2'. This en- 
ables averages over background galaxies without cancelling 
the cosmic shear signal. If y varied on scales comparable to 
or smaller than the mean galaxy separation, any average over 
galaxies would remove the lensing signal as well. 

The intrinsic ellipticities of the faint background galaxies 
have a distribution with a standard deviation of <x e ss 0.3. 
Averaging over N of them, and assuming Poisson statistics, 
yields expectation values of 



<e (s) ) = , 5e = <(e (s) ) 2 > 1/2 = 



(241) 



for the mean and its intrinsic fluctuation. 

A rough estimate for the signal-to-noise ratio of a cosmic 
shear measurement can proceed as follows. Suppose the cor- 
relation function £ is measured by counting pairs of galaxies 
with a separation within 59 of 9. As long as 6 is small com- 
pared to the side length of the survey area A, the number of 
pairs will be 



1 9 

N p = - xn 2 Ax2n969 , 



(242) 



and thus the Poisson noise due to the intrinsic ellipticities will 
be 



2a f 



noise 



n ^nA959 



(243) 



where the factor of two arises because of the two galaxies in- 
volved in each pair. The signal is the square root of the corre- 
lation function £ which we can approximate as 

f * l 2 P K (l)6 In / » l 2 P K (l) j * l 2 PM j , (244) 

where we have used in the last step that 9 = 2n/l. Thus, the 
signal-to-noise ratio is estimated to be 



S_ 

N 



V? lnd9yf^AP~ K n y/n 3 AP K 69 



noise 



2<x, 



(245) 



Evidently, the signal-to-noise ratio, and thus the significance 
of any cosmic-lensing detection, grows with the survey area 
and decreases with the intrinsic ellipticity of the source galax- 
ies. In evaluating (245 1 numerically, we have to take into ac- 
count that l 2 P K (l) must be a dimension-less number, which 
implies that the power spectrum P K must have the dimension 
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steradian. Therefore, either the survey area A and the number 
density n in \2A5\ must be converted to steradians, or P K must 
be converted to square arc minutes first. 

The signal-to-noise ratio increases approximately linearly 
with scale. Assuming 60/9 = 0.1, it is S /N ~ 5 on a scale 
of 0.5' for a survey of one square degree area. This shows 
that, if the cosmic shear should be measured on such small 
scales with an accuracy of, say, five per cent, a survey area of 
A w (20/5) 2 « 16 square degrees is needed since the signal- 
to-noise ratio scales like the square root of the survey area. 
On such an area, the ellipticities of 16 x 3600 x 30 ~ 2 x 10 6 
background galaxies have to be accurately measured. 

Matters are more complicated in reality, but the orders- 
of-magnitude are well represented by this rough estimate. 
Bearing in mind that typical fields-of-view of telescopes 
large enough for detecting sufficiently many faint background 
galaxies reach one tenth of a square degree up to one square 
degree, and that typical exposure times are of order one hour 
for that purpose, the total amount of telescope time for a weak- 
lensing survey of about 100 square degrees is estimated to be 
several hundred telescope hours. With perhaps five hours of 
telescope time per night, and perhaps half of the nights per 
year usable, it is easy to see that the time needed for such 
surveys is measured in months. 

Since typical sizes of the faint background galaxies mea- 
sure fractions of arc seconds, shape measurements require a 
pixel resolution of, say, 0.1". A total survey area of 100 
square degrees must therefore be resolved into 100 x 3600 x 
3600/0. 1 2 w 1.3 X 10 11 pixels. 

These estimates neglect all sources of noise other than the 
shot noise caused by the finite number of galaxies. On angu- 
lar scales below a few arc minutes, the cosmic variance caused 
by field-to-field variations in the shear signal due to the large- 



scale cosmic structures must be added ( Kilbinger and Schnei- 
|der]|2004) . 



2. Ellipticity measurements 

The determination of image ellipticities is straightforward 
in principle, but difficult in practice ( |Kaiser et al.\ 1 1995| >. Of- 
ten, the surface-brightness quadrupole 



J I(x)x iX jd 2 
j I(x)d 2 x 



(246) 



is measured, from whose principal axes the ellipticity can be 
read off. Real galaxy images, however, are typically far from 
ideally elliptical. They are structured or otherwise irregular. 
In addition, if they are small, they are coarsely resolved into 
just a few pixels, so that only a crude approximation to the 



integral in ( 246 1 can be found. 

How the image of a point-like source, such as a star, ap- 
pears on the detector is described by the point-spread function 
(PSF). The PSF may be anisotropic if the telescope optics is 
slightly astigmatic, and this anisotropy may, and will in gen- 
eral, depend on the location in the focal plane. The image is a 
convolution of the ideal image shape prior to any distortion by 



the atmosphere and the telescope optics. Any accurate mea- 
surement of image ellipticities requires a PSF correction or 
deconvolution, for which the PSF must of course be known. It 
is commonly measured off the images of stars in the field. 

Many other effects may distort the PSF and thus the images 
in systematic ways. For instance, if the CCD chips are not 
exactly perpendicular to the optical axis of the telescope, or if 
the individual chips of a CCD mosaic are not exactly copla- 
nar, or if the telescope is slightly out of focus, systematic im- 
age deformations may result which typically vary across the 
focal plane. They have to be measured and corrected. This 
is commonly achieved by fitting the measured PSF by low- 
order, two-dimensional polynomials on the focal plane. Since 
part of the image distortions may depend on time due to ther- 
mal deformation, changing atmospheric conditions and such, 
PSF corrections will also depend on time and have to be de- 
termined and applied with much care. 

Even if the surface-brightness quadrupole of the image on 
the detector can be accurately determined, the image appears 
affected by imperfections of the telescope optics and by the 
turbulence in the atmosphere, the so-called seeing. Due to 
the wave nature of light and the finite size of the telescope 
mirror, the telescope will h ave fi nite resolution. The angular 
resolution limit is given by ( 185 1. With A « 6 x 10~ 5 cm and 
D = 400 cm, the angular resolution is AO w 0.04", much 
smaller than needed for our purposes. 

The turbulence of the Earth's atmosphere effectively con- 
volves images with a Gaussian whose width depends on the 
site, the weather and other conditions. Typical seeing ranges 
around 1". Under very good conditions, it can shrink to 
~ 0.5" or even less. Clearly, if an image of sub-arc second 
size is convolved with a Gaussian of similar width, any ellip- 
ticity is substantially reduced. 

Systematic effects may remain which need to be detected 
and quantified. Any coherent image distortions caused by 
gravitational lensing must be describable by the tidal gravi- 
tational field, i.e. by second-order derivatives of a scalar po- 
tential. In analogy to the fi-field in electromagnetism, such 
distortion patterns are called E-modes. Similarly, distortion 
patterns which are described as the curl of a vector field are 
called B-modes. They cannot be due to gravitational lensing 
and thus signal systematic effects remaining in the data. Such 
fi-mode contaminations could recently be strongly reduced or 



suppressed by improved algorithms for PSF correction ( |Hoek- 
|stra| |2004| l. Absence of B-mode contamination does not al- 
low the implication that the results are free of systematics, 
though, because optical distortions also tend to create spuri- 
ous E-modes. 



3. Results 

Despite the smallness of the effect and the many difficul 
ties in measuring it, much progress in cosmic-shear observa 
tions has been achieved in the past few years (Heymans et al. 



2006| |Massey et al. 2007a| l. Current and ongoing surveys 
in particular the Canada-France-Hawaii Legacy Survey (Fu 
et al.\ |2008[ |Semboloni et al.\ |2006| l, combined with well 
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developed, largely automatic data-analysis pipelines, have 
succeeded in producing cosmic-shear correlation functions 
with very small error bars covering angular scales from below 
an arc minute to several degrees. The best correlation func- 
tions could be shown to be at most negligibly contaminated 
by B-modes. 

The power spectrum P K {1) depends crucially on the non- 
linear evolution of the dark-matter power spectrum. This, and 
the exact redshift distribution of the background galaxies, are 
the major uncertainties now remaining in the interpretation of 
cosmic-shear surveys. Apart from that, the measured cosmic- 
shear correlation functions agree very well with the theoretical 
expectation from CDM density fluctuations in a spatially-flat, 
low-density universe. 



As (230 1 shows, the weak-lensing power spectrum P K (l) 
depends on the product of a factor Q^ due to the Poisson 
equation, times the amplitude A of the matter power spec- 
trum. An additional weak dependence on cosmological pa- 
rameters is caused by the geometric weight function W(w', w). 
The cosmic-shear correlation function thus measures approx- 
imately the product AQj^L a < 1, which means that the am- 
plitude of the power spectrum is nearly degenerate with the 
matter density parameter. Only if it is possible to constrain 
Q. m o or A in any other way can the degeneracy be broken. 

We shall see later how this may work. The amplitude of the 
power spectrum A is conventionally described by a parame- 
ter 0"^ which will be defined and described in more detail in 
§ |X| Weak lensing thus constrains the product crgQ^, and lat- 
est measure ments find cr 8 (Q m o/0.25) 064 ~ 0.784 + 0.043 |Fu| 
[eFari[2008l l. 

Weak gravitational lensing is a fairly new field of cosmo- 
logical research. Within a few years, it has considerably 
matured and returned cosmologically interesting constraints. 
Considerable potential is attributed to weak lensing in wide- 
area surveys in particular when combined with photometric 
redshift information, because this is expected to allow con- 
straints on the growth of cosmic structures. 



IX. SUPERNOVAE OF TYPE IA 

A. Standard candles and distances 

1 . The principle 

Before discussing supernovae of type la and their cosmo- 
logical relevance, let us set the stage with a few illustrative 
considerations. Suppose we had a standard candle whose lu- 
minosity, L, we knew precisely. Then, according to the defi- 



nition of the luminosity distance in ( 20 1, the distance can be 
inferred from the measured flux, S , through 



D] um 
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(247) 



Besides the redshift z, the luminosity distance will depend on 
the cosmological parameters, 
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FIG. 10 Constraints in the fi m o-cr 8 plane from weak lensing on large 
angular scales in the CFHTLS. The Universe is assumed spatially flat 
here, (from ( |Fu et q/.||2008l >) 

which can in principle be used to determine cosmological pa- 
rameters from a set of distance measurements from a class of 
standard candles. 

For this to work, the standard candles must be at a suit- 
ably high redshift for the luminosity distance to depend on the 
cosmological model. As we have seen in (21 1, all distance 
measures share the low-z limit 
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(249) 



and lose their sensitivity to all cosmological parameters except 

In reality, we rarely know the absolute luminosity L even of 
cosmological standard candles. The problem is that they need 
to be calibrated first, which is only possible from a flux mea- 
surement once the distance is known by other means, such as 
from parallaxes in case of the Cepheids. Supernovae, how- 
ever, which are the subject of this chapter, are typically found 
at distances which are way too large to allow direct distance 
measurements. Therefore, the only way out is to combine dis- 
tant supernovae with local ones, for which the approximate 
distance relation ( |249| l holds. 

Any measurement of flux S ; and redshift z; of the z'-th stan- 
dard candle in a sample then yields an estimate for the lumi- 
nosity L in terms of the squared inverse Hubble constant, 



L = 4nSi 



(250) 



Since all cosmological distance measures are proportional to 
the Hubble length c/Hq, the dependence on Hq on both sides 



Aum - Aumfc £^m0> ^A0, Hq, . . .) , 



(248) 



of ( 247 1 cancels, and the determination of cosmological pa- 
rameters other than the Hubble constant becomes possible. 
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Thus, the first lesson to learn is that cosmology from distant 
supernovae requires a sample of nearby supernovae for cali- 
bration. 

Of course, this nearby sample must satisfy the same crite- 
rion as the distance indicators used for the determination of 
the Hubble constant: their redshifts must be high enough for 
the peculiar velocities to be negligible, thus z > 0.02. If the 



redshifts are low enough for the linear approximation ( 249 1 to 
hold, the interpretation of the nearby sample is independent of 
the cosmological model. 

It is important to note that it is not necessary to know the 
absolute luminosity L even up to the uncertainty in Hq. If 
L is truly independent of redshift, cosmological parameters 



could still be determined through ( 247 1 from the shape of the 
measured relation between flux and redshift even though its 
precise amplitude may be unknown. It is only important that 
the objects used are standard candles, but not how bright they 



2. Requirements and degeneracies 

Let us now collect several facts about cosmological infer- 
ence from standard candles. Since we aim at the determina- 
tion of cosmological parameters, say Q m o, it is important to 
estimate the accuracy that we can achieve from measurements 
of the luminosity distance. Suppose we restrict the attention to 
spatially flat cosmological models, for which D A o = 1 — Q m Q. 
Then, because the dependence on the Hubble constant was 
cancelled, Q m o is the only remaining relevant parameter. We 
estimate the accuracy through first-order Taylor expansion, 
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about a fiducial model, such as a ACDM model with Q m0 
0.3. 

At a fiducial redshift of z ~ 0.8, we find numerically 
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which shows that a relative distance accuracy of 
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is required to achieve an absolute accuracy of AQ m o. For 
AQ m o m 0.02, say, distances thus need to be known to w 1%. 

This accuracy requires sufficiently large supernova sam- 
ples. Assuming Poisson statistics for simplicity and distance 
measurements to N supernovae, the combined accuracy is 
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That is, an accuracy of AQ m o ~ 0.02 can be achieved from 
« 100 supernovae whose individual distances are known to 
* 10%. 

Anticipating physical properties of type-la supernovae, 
their intrinsic peak luminosities in blue light are L m 3.3 x 



10 erg s , with a relative scatter of order 10%. Given un- 
certainties AL in the luminosity L and AS in the flux measure- 
ment S , error propagation on ( |247[ i yields the relative distance 
uncertainty 
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Even if the flux could be measured precisely, the intrinsic lu- 
minosity scatter currently forbids distance determinations to 
better than 10%. 

Fluxes have to be inferred from photon counts. For vari- 
ous reasons to be clarified later, supernova light curves should 
be determined until ~ 35 days after the peak, when the lu- 
minosity has typically dropped to w 2.5 x 10 42 ergs _1 . The 
luminosity distance to z ~ 0.8 is » 5 Gpc, which implies 
fluxes S * 1.1 x 10~ 14 ergs~' cirT 2 at peak and S ~ 8.7 x 
10~ 16 erg s _1 cirT 2 35 days later. 

Dividing by an average photon energy of 5 x 10~ 12 erg, mul- 
tiplying with the area of a typical telescope mirror with 4 m 
diameter, and assuming a total quantum efficiency of 30%, 
we find detected photon fluxes of S Y ~ 85 s~' at peak and 
Sy ~ 7 s _1 35 days afterwards. These fluxes are typically dis- 
tributed over a few CCD pixels. 

Supernovae need to be identified against the background of 
their host galaxies and the sky brightness. For distant super- 
novae, the sky brightness on the area of the supernova image 
is typically 10-1000 times higher than the supernova itself, 
dominating the noise budget. Let N$n an d N s ^ y me number 
of photons received from the supernova and the sky, respec- 
tively, on the supernovae image, and assume for simplicity 
A^ s ky ~ IOOA^sn- Then, an estimate for the signal-to-noise ra- 
tio for the detection is 
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Signal-to-noise ratios of > 10 up to 35 days after the max- 
imum thus require A^sn ~ 10 4 photons from the supernova. 
This implies exposure times of order 10 4 /7 w 1400 s, or 
about 20 minutes. Typical exposure times are thus of or- 
der (15 . . . 30) minutes to capture supernovae out to redshifts 
Z ~ 1. Then, the photometric error around peak luminosity 
is certainly less than the remaining scatter in the intrinsic lu- 
minosity, and relative distance accuracies of order 10% are 
within reach. 

However, a major difficulty is the fact that the identifica- 
tion of type-la supernovae requires spectroscopy. Sufficiently 
accurate spectra typically require exposures lasting 1-3 hours 
per spectrum on the world's largest telescopes, such as ESO's 
Very Large Telescope which consists of four individual mir- 
rors with 8 m diameter each. 



3 As we shall see later, type-la supernovae are standardisable rather than 
standard candles, and the standardising procedure is currently not able to 
reduce the scatter further. 
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In order to see what we can hope to constrain by measuring 
luminosity distances, we form the gradient of Di um in the Q m Q- 
Qao plane, 



g = 
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at a fiducial ACDM model with O m o = 0.3. When normalised 
to unit length, it turns out to point into the direction 
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This vector rotated by 90° then points into the direction in 
the Q m o-QAo plane along which the luminosity distance does 
not change. Thus, near the fiducial ACDM model, the param- 
eter combination 
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is degenerate. The degeneracy direction, characterised by 
the vector %(njT)g - (0.65,0.76) r , encloses an angle of 
arctan(0.76/0.65) = 49.5° with the Q m o axis, almost along 
the diagonal from the lower left to the upper right corner of 
the parameter plane. Thus, it is almost perpendicular to the 
degeneracy direction obtained from the curvature constraint 
due to the CMB. This illustrates how parameter degeneracies 
can very efficiently be broken by combining suitably different 
types of measurement. Moreover, combining supernova data 
from a wide redshift range partially lifts the parameter degen- 
eracy obtained from them already. 



B. Supernovae 

1 . Types and classification 

Supernovae are "eruptively variable" stars. A sudden rise 
in brightness is followed by a gentle decline. They are 
unique events which at peak brightness reach luminosities 
comparable to those of an entire galaxy, or (10 10 . . . 10 11 ) L Q . 
They reach their maxima within days and fade within several 
months. 

Supernovae are traditionally characterised by their early 
spectra ( Filippenko 1997). If hydrogen lines are missing, they 
are of type I, otherwise of type II. Type-la supernovae show 
silicon lines, unlike type-Ib/c supernovae, which are distin- 
guished by the prominence of helium lines. Normal type-II 
supernovae have spectra dominated by hydrogen. They are 
subdivided according to their lightcurve shape into type-IIL 
and type-IIP. Type-lib supernova spectra are dominated by he- 
lium instead. 

Except for type-la, supernovae arise due to the collapse of 
a massive stellar core, followed by a thermonuclear explo- 
sion which disrupts the star by driving a shock wave through 
it. Core-collapse supernovae of type-II arise from stars with 
masses between (8 . . . 30) M , those of type-I (i.e. types Ib/c) 
from more massive stars ( Woosley et al. 2002| l. 



Type-la supernovae, which we are dealing with here, arise 
when a white dwarf is driven towards the Chandrasekhar mass 
limit by mass overflowing from a companion star. In a binary 
system, the more massive star evolves faster and can reach its 
white-dwarf stage before its companion leaves the main se- 
quence and becomes a red giant. When this happens, and the 
stars are close enough, matter will flow from the expanding 
red giant on the white dwarf. 

Electron degeneracy pressure can stabilise white dwarfs up 
to the Chandrasekhar mass limit of ~ 1.4 M Q (Chandrasekhar 



1984 Hillebrandt and Niemeyer 2000). Because the core ma- 
terial is degenerate, its pressure is independent of its temper- 
ature. The mass accreted from the companion star increases 
the pressure until nuclear burning can begin in isolated places 
somewhere in the core. The electron degeneracy is lifted, 
the temperature rises dramatically, and the thermonuclear run- 
away sets in. Neutrinos produced in inverse beta decays carry 
away much of the explosion energy unnoticed because they 
can leave the supernova essentially without further interac- 
tion. 

This thermonuclear runaway destroys the white dwarf. 
Since this type of explosion (more precisely, deflagration) in- 
volves an approximately fixed amount of mass, it is physically 
plausible that the explosion releases a fixed amount of energy. 
Thus, the Chandrasekhar mass limit is the main responsible 
for type-la supernovae to be approximate standard candles. 

The nuclear fusion processes in type-la supernovae con- 
verts the carbon and oxygen in the core of the white dwarf into 
56 Ni, which later decays through 56 Co into stable 56 Fe. Ac- 
cording to detailed numerical explosion models, the nuclear 
fusion is started at random points near the centre of the white 



dwarf (Ro pkeef al.\\2QQ7\ . 

The presence of silicon lines in the type-la spectra indicates 
that not all of the white dwarf's material is converted into 56 Ni. 
This shows that there is no explosion, but a deflagration, in 
which the flame front propagates at velocities below the sound 
speed. The deflagration can burn the material fast enough if 
it is turbulent, because the turbulence dramatically increases 
the surface of the flame front and thus the amount of material 
burnt per unit time. Theoretical models predict that ~ 0.5 M 
of 56 Ni are typically produced. 

The peak brightness is reached when the deflagration front 
reaches the former white dwarf's surface and drives it as a 
rapidly expanding envelope into the surrounding space. The y 
photons released in the nuclear fusion processes are redshifted 
by scattering off the expanding material and finally leave the 
explosion site as X-ray, UV, optical and infrared photons. 

Once the thermonuclear fusion has ended, additional en- 
ergy is released by the /3 decay of 56 Co into 56 Fe with a half 
life of 77.12 days. The exponential nature of the radioactive 
decay causes the typical exponential decline phase in super- 
nova light curves (Colgate and McKee 1969 Truran et al. 



1967 1. Since the supernova light has to propagate through the 
expanding envelope before we can see it, the opacity of the en- 
velope and thus its metalicity are important for the appearance 



of the supernova ( Blinnikov et al. 2006 ). 
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2. Observations 

Since supernovae are transient phenomena, they can only 
be detected by sufficiently frequent monitoring of selected ar- 
eas in the sky. Typically, fields are selected by their accessi- 
bility for the telescope to be used and the least degree of ab- 
sorption by the Galaxy. Since a type-la supernova event lasts 
for about a month, monitoring is required every few days. 

Supernovae are then detected by differential photometry, in 
which the average of all preceding images is subtracted from 
the last image taken. Since the seeing varies, the images ap- 
pear convolved with point-spread functions of variable width 
even if they are taken with identical optics, thus the objects on 
them appear more or less blurred. Before they can be mean- 
ingfully subtracted, they therefore have to be convolved with 
the same effective point-spread function. This causes sev- 
eral complications in the later analysis procedure, in particular 
with the photometry. 

Of course, this detection procedure returns many variable 
stars and supernovae of other types, which are not standard 
candles and have to be removed from the sample. Pre- 
selection of type-la candidates is done by colour and the light- 
curve shape, but the identification of type-la supernovae re- 
quires spectroscopy in order to identify the decisive silicon 
lines at 6347 A and 6371 A. Since these lines move out of the 
optical spectrum at redshifts z > 0.5, near-infrared observa- 
tions are crucially important for the high-redshift supernovae 
relevant for cosmology. 

Nearby supernovae, which are needed for calibration, re- 
veal that type-la supernovae are not standard candles, but 
show a substantial scatter in luminosity. It turned out that 
there is an empirical relation between the duration of the su- 
pernova event and its peak brightness in that brighter super- 



novae last longer (Phillips et al. 1999 1. This relation between 



the light-curve shape and the brightness can be used to stan- 
dardise type-la supernovae. It was seen as a major problem 
for their cosmological interpretation that the origin for this re- 
lation was unknown, and that its application to high-redshift 
supernovae was based on the untested assumption that the re- 
lation found and calibrated with local supernovae would also 
hold there. Recent simulations indicate that the relation is 
an opacity effect ( |Kasen and Woosley 2007 1: brighter su- 
pernovae produce more 3b Ni and thus have a higher metalic- 
ity, which causes the envelope to be more opaque, the energy 
transport through it to be slower, and therefore the supernova 
to last longer. 

Thus, before a type-la supernova can be used as a standard 
candle, its duration must be determined, which requires the 
light-curve to be observed over sufficiently long time. It must 
be taken into account here that the cosmic expansion leads 
to a time dilation, due to which supernovae at redshift z ap- 
pear longer by a factor of (1 + z). We note in passing that 
the confirmation of this time dilation effect indirectly supports 
the cosmic expansion. After the standardisation, the scatter 
in the peak brightnesses of nearby supernovae is substantially 
reduced. This encourages (and justifies) their use as standard- 
isable candles for cosmology. 

The remaining relative uncertainty is now typically between 



(10... 15)% for individual supernovae. Since, as we have 
seen following ( |253| l, we require relative distance uncertain- 
ties at the per cent level, of order a hundred distant supernovae 
are needed before meaningful cosmological constraints can be 



placed, which justifies the remark after (254i 



An example for the several currently ongoing supernova 



surveys is the Supernova Legacy Survey (SNLS, ( |Astier et al. 



2006 )) in the framework of the Canada-France-Hawaii Legacy 



Survey (CFHTLS), which is carried out with the 4-m Canada- 
France-Hawaii telescope on Mauna Kea. It monitors four 
fields of one square degree each five times during the 18 days 
of dark time between two full moons (lunations). 

Differential photometry is performed to find out variables, 
and candidate type-la supernovae are selected by light-curve 
fitting after removing known variable stars. Spectroscopy on 
the largest telescopes (mostly ESO's VLT, but also the Keck 
and Gemini telescopes) is then needed to identify type-la su- 
pernovae. To give a few characteristic numbers, the SNLS has 
taken 142 spectra of type-la candidates during its first year of 
operation, of which 91 were identified as type-la supernovae. 

The light curves of these objects are observed in several 
different filter bands. This is important to correct for interstel- 
lar absorption. Any dimming by intervening material makes 
supernovae appear fainter, and thus more distant, and will 
bias the cosmological results towards faster expansion. Since 
the intrinsic colours of type-la supernovae are characteristic, 
any deviation of the observed from the intrinsic colours sig- 
nals interstellar absorption which is corrected by adapting the 
amount of absorption such that the observed is transformed 
back into the intrinsic colour. 

This correction procedure is expected to work well unless 
there is material on the way which absorbs equally at all wave- 
lengths, so-called "grey dust" ( |Aguirre| |1999| l. This could 
happen if the absorbing dust grains are large compared to the 
wavelength. Currently, it is quite difficult to conclusively rule 
out grey dust, although it is implausible based on the interstel 



lar absorption observed in the Galaxy ( |Ostman and Mortsell 
|2005l|RiesseFar||2004) . 

After applying the corrections for absorption and duration, 
each supernova yields an estimate for the luminosity distance 
to its redshift. Together, the supernovae in the observed sam- 
ple constrain the evolution of the luminosity distance with 
redshift, which is then fit varying the cosmological parame- 
ters except for Hq, i.e. typically Q m o and Qao- This yields an 
"allowed" region in the Q m Q-Q A0 plane compatible with the 
measurements which is degenerate in the direction calculated 
before. 

More information or further assumptions are necessary to 
break the degeneracy. The most common assumption, justi- 
fied by the CMB measurements, is that the Universe is spa- 
tially flat. Based upon it, the SNLS data yield a matter density 
parameter of 



£2 m0 = 0.263 + 0.037 



(260) 



This is a remarkable result. First of all, it confirms the other 
independent measurements we have already discussed, which 
were based on kinematics, cluster evolution and the CMB. 
Second, it shows that, in the assumed spatially flat universe, 
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FIG. 1 1 Constraints in the £2 m0 -fl A o plane obtained from type-la su- 
pernovae in the SNLS. (from ( |Astier et al. 20061). 



the dominant contribution to the total energy density must 
come from something else than matter, possibly the cosmo- 
logical constant. 

It is important for the later discussion to realise in what way 
the parameter constraints from supernovae differ from those 
from the CMB. The fluctuations in the latter show that the 
Universe is at least nearly spatially flat, and the density pa- 
rameters in dark and baryonic matter are near 0.24 and 0.04, 
respectively. The rest must be the cosmological constant, or 
the dark energy. Arising early in the cosmic history, the CMB 
itself is almost insensitive to the cosmological constant, and 
thus it can only constrain it indirectly. 

Type-la supernovae, however, measure the angular- 
diameter distance during the late cosmic evolution, when the 
cosmological constant is much more important. As ( 259 1 
shows, the luminosity distance constrains the difference be- 
tween the two parameters, 



where the degenerate parameter P is determined by the mea- 
surement. Assuming Oao = 1 - Q m o as in a spatially-flat 
universe yields 



p = i -2.no. 
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0.43 



(262) 



from the SNLS first-year result ( 260 1, illustrating that the sur- 



vey has constrained the density parameters to follow the rela- 
tion 



1.17 Q m0 + 0.43 



(263) 



The relative acceleration of the universe, a I a, is given by 
|7J, which can be simplified to read 



- =H \Cl A0 ~ — 



a 



(264) 



if matter is pressure-less. Thus, the expansion of the universe 
accelerates today (a = 1) if a = Hq(Q.ao - ^mo/2) > 0, or 
£2 A o > C2 m o/2. Given the measurement ( |263| l, the conclusion 
seems inevitable that the Universe's expansion does indeed 
accel erate today (jAstier et aT||2006||Leib undgut and Blondin 
[2Ti03i|Perimutter et a/.||1999[|Riess et q/.||1998) . 

If the Universe is indeed spatially flat, then the transition 
between decelerated and accelerated expansion happened at 



1 - 0.263 
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a = 0.56 



(265) 



or at redshift z ~ 0.78. Luminosity distances to supernovae at 
higher redshifts should show this transition, and in fact they 
do ( |Riess ef ari|2004l l. 



3. Potential problems 

The problem with possible grey dust has already been men- 
tioned: While the typical colours of type-la supernovae al- 
low the detection and correction of the reddening coming with 
typical interstellar absorption, grey dust would leave no trace 
in the colours and remain undetectable. However, grey dust 
would re-emit the absorbed radiation in the infrared and add 
to the infrared background, which is quite well constrained. It 
thus seems that grey dust is not an important contaminant, if 
it exists. 

Gravitational lensing is inevitable for distant supernovae 



£2a () = 1.17£2 m0 + P. 



(261) 



( Holz and Wald|[1998[ ). Depending on the line-of-sight, they 
are either magnified or demagnified. Since high magnifica- 
tions due to non-linear structures may occasionally happen, 
the magnification distribution must be skewed towards de- 
magnification to keep the mean of unit magnification. Thus, 
the most probable magnification experienced by supernova is 
below unity. In other words, lensing may lead to a slight de- 
magnification if lines-of-sight towards type-la supernovae are 
random with respect to the matter distribution. In any case, 
the rms cosmic magnification adds to the intrinsic scatter of 
the supernova luminosities. It may become significant for red- 
shifts z > 1. 

It is a difficult and debated question whether supernovae 
at high redshifts are intrinsically the same as at low redshifts 
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where they are calibrated. Should there be undetected sys- 
tematic differences, cosmological inferences could be wrong. 
In particular, it may be natural to assume that metalicities at 
high redshifts are lower than at low redshifts. Since super- 
novae last longer if their atmospheres are more opaque, lower 
metalicity may imply shorter supernova events, leading to un- 
derestimated luminosities and overestimated distances. Sim- 
ulations of type-la supernovae, however, seem to show that 
such an effect is probably not significant. For this and other 
systematic effects, see (fLeibundgut , .200 1) . 

It was also speculated that distant supernovae may be in- 
trinsically bluer than nearby ones due to their possibly lower 
metalicity. Should this be so, the extinction correction, which 
is derived from reddening, would be underestimated, causing 
intrinsic luminosities to be under- and luminosity distances to 
be overestimated. Thus, this effect would lead to an underes- 
timate of the expansion rate and counteract the cosmological 
constant. There is currently no indication of such a colour 
effect. 

Supernovae of types Ib/c may be mistaken for those of type 
la if the identification of the characteristic silicon lines fails 
for some reason. Since they are typically fainter than type- 
la supernovae, they would contaminate the sample and bias 
results towards higher luminosity distances, and thus towards 
a higher cosmological constant. It seems, however, that the 
possible contamination by non-type-la supernovae is so small 
that it has no noticeable effect. 

Several more potential problems exist. It has been argued 
for a while that, if the evidence for a cosmological constant 
was based exclusively on type-la supernovae, it would prob- 
ably not be considered entirely convincing. However, since 
the supernova observations come to conclusions compatible 
with virtually all independent cosmological measurements, 
they add substantially to the persuasiveness of the cosmolog- 
ical standard model. Moreover, recent supernova simulations 
reveal good physical reasons why they should in fact be reli- 
able, standardisable candles. 



X. THE NORMALISATION OF THE POWER SPECTRUM 
A. Introduction 



quite well established, its amplitude still poses a surprisingly 
obstinate problem. We shall see in this section why it is so 
difficult to measure. For this purpose, we shall discuss four 
ways of measuring <x 8 ; the amplitude of large-scale tempera- 
ture fluctuations in the CMB, the cosmic-shear autocorrelation 
function, the abundance and evolution of the galaxy-cluster 
population, and the statistics of Lyman-ff forest lines. 

For historical reasons, the amplitude of the dark-matter 
power spectrum is characterised by the variance of the density 
fluctuations within spheres of 8/r'Mpc radius. More gener- 
ally, one imagines placing spheres of radius R randomly and 
measuring the density-contrast variance within them. Since 
the variance in Fourier space is characterised by the power 
spectrum, it can be written as 



oi= f 
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where Wjt(k) is a window function selecting the k modes con- 
tributing to the variance within the spheres. 

Imagining spheres of radius R in real space, the window 
function should be the Fourier transform of a step function, 
which is inconvenient because it extends to infinite wave 
numbers. It is thus more common to use either Gaussians, 
since they Fourier transform to Gaussians, or step functions 
in Fourier space. For simplicity of the following illustrative 
calculations, we use the latter choice, thus 
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This is a step function dropping to zero at k = 2n/R. Inserting 



this into ( 266 1 gives 
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All modes larger than R contribute to the density fluctuations 
in spheres of radius R because all smaller modes average to 
zero. The normalisation of the power spectrum is usually 
expressed in terms of <x 8 , fixing R to its historical value of 
8 /T 1 Mpc. 



B. Fluctuations in the CMB 



We saw in § VII.B that the measured power spectrum of the 
galaxy distribution follows the CDM expectation in the range 
of wave numbers where current large surveys allow its mea- 
surement. This range can be extended to some degree towards 
smaller scales by measuring the autocorrelation of hydrogen 
absorption lines in the spectra of distant quasars. Such obser- 
vations of the power spectrum of the so-called Lyman-o- forest 
lines show that the power spectr um does indeed turn towards 
the asymptotic behaviour oc £ 3 ([McDonald et al. 2005 Viel 



et al. |2008[ ). In addition, we have seen that the peak location 
agrees with the expectation for a universe with Q. m o as 0.3 and 
h as 0.72. This indicates that the CDM expectation for the 
dark-matter power spectrum is indeed at least very close to its 
real shape, which is a remarkable finding. 

Although the shape of the power spectrum could thus be 



1 . The large-scale fluctuation amplitude 



We saw in § VLB. 3 that the long-wavelength (low-k) tail 



of the CMB power spectrum is caused by the Sachs- Wolfe 
effect, giving rise to relative temperature fluctuations ( |171| i in 
terms of the Newtonian potential fluctuations 50. The three- 
dimensional temperature-fluctuation power spectrum is then 
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The Poisson equation in its form (219 1 implies that the 
power spectra of potential- and density fluctuations are related 
through 



p m 9// o n 2 (D + (a KC )\ 2 P s {k) 



/c 4 



(270) 
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where the linear growth factor D + (a mc ) was introduced to re- 
late the potential-fluctuation power spectrum at the time of 
decoupling to the present density-fluctuation power spectrum 
I'M). 

Now we need to account for projection effects. A three- 
dimensional mode with comoving wave number k and comov- 
ing wavelength A = 2n/k appears under an angle = A/w rec , 
where w KC = w(a Tec ) is the comoving angular-diameter dis- 
tance (176 1 to the CMB. The angular wave number under 



which the mode appears is thus 



2n 



(271) 



Expressing now the power spectrum ( 269 1 in terms of the 



angular wave number I yields, with (270 1 
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where the factor w TI ? c arises because of the transformation 
from spatial to angular wave numbers I [cf. Limber's equation 



(228 1], and the factor w^ ec /l expresses the factor k from the 



squared Laplacian. This shows that the angular power spec- 
trum P T (l) of the large-scale CMB temperature fluctuations 
can only be translated into the amplitude of the dark-matter 
power spectrum A if the cosmological parameters are already 
known well enough. 

A further complication is added by the integrated Sachs- 



Wolfe effect introduced in § VLB. 6 It depends on D+(a)/a 
and adds secondary anisotropies to the CMB unless D + (a) = 
a. The primordial CMB fluctuations are then lower than mea- 
sured and need to be corrected by subtracting the integrated 
Sachs-Wolfe contribution, which adds a further dependence 
on the cosmological parameters. 



2. Translation to tr 8 

Two more complications arise in the translation of the large- 
scale amplitude A to cr 8 . Since structures with wave numbers 
larger than the peak location k eq in the power spectrum con- 
tribute to <t%, the dependence of k eq on the cosmological pa- 
rameters comes in. Finally, the spectral index n s may affect 
the extrapolation from large to small scales substantially be- 
cause of the long lever arm between the scales involved. 

Of course, one could also use the small-scale part of the 
CMB power spectrum for normalising the dark-matter power 
spectrum. Due to the acoustic oscillations, however, this part 
depends in a much more complicated way on additional cos- 
mological parameters, such as the baryon density. Reading 
cr 8 off the low-order multipoles is thus a safer, albeit intricate 
procedure. 

Even if the cosmological parameters are now known well 
enough to translate the low-order CMB multipoles to <x 8 , 
an additional uncertainty remains. We know that, although 
the Universe became neutral ~ 400, 000 years after the Big 
Bang, it must have been reionised after the first stars and 



other sources of UV radiation formed. Since then, CMB pho- 
tons are travelling through ionised material again and experi- 
ence Thomson (or Compton) scattering. The optical depth for 
Thomson scattering is 



o" T cdt 



(273) 



where n e is the number density of free electrons and cr T is the 
Thomson scattering cross section. After propagating through 
the optical depth r, the CMB fluctuation amplitude is reduced 
by exp(-r). 

Of course, the CMB photons cannot disappear through 
Thomson scattering, thus the CMB's overall intensity cannot 
change in this way, but the fluctuation amplitudes are lowered 
in this diffusion process. The optical depth r depends on the 
path length through ionised material. In view of the CMB, 
this means that the degree of fluctuation damping depends on 
the reionisation redshift, i.e. the redshift after which the cos- 
mic baryons were transformed back into a plasma. Unless the 
reionisation redshift is known, we cannot know by how much 
the CMB fluctuations were suppressed. 

So far, the reionisation redshift can be estimated in two 



ways. First, as discussed in § VI.B.4 Thomson scattering 
creates linear polarisation. Of course, the polarisation due 
to reionised material appears superposed on the primordial 
polarisation, but on different angular scales. The character- 
istic scale for secondary polarisation is the horizon size at 
the reionisation redshift, which is much larger than the typ- 
ical scales of the primordial polarisation. Thus, the reionisa- 
tion redshift can be inferred from large-scale features in the 
CMB polarisation, provided the cosmological parameters are 
known well enough to translate angular scales into physical 
scales (ZaldarriagaJ |1997| >. Modern cosmological parameter 
estimates aim at fitting the available data simultaneously for 
all relevant parameters mostly by Monte-Carlo Markov chain 



techniques, as explained in detail e.g. in the textbook by Dur- 

Irerl2008l 

Unfortunately, this is aggravated by the polarised mi- 
crowave radiation from the Milky Way. Synchrotron and dust 
emission can be substantially polarised and mask the CMB 
polarisation, which can only be measured reliably if the fore- 



grounds of Galactic origin can be accurately subtracted (Page 



et al 



2007 ). Thus, the degree to which the foreground polari- 



sation is known directly determines the accuracy of the <x 8 pa- 
rameter derived from the CMB fluctuations. This contributes 
considerably to the remaining uncertainty in the cr 8 derived 
from the 5-year WMAP data given in Tab. [I] 

The other way to constrain the reionisation redshift uses 
the spectra of distant quasars. Light with wavelengths shorter 
than the Lyman-a wavelength cannot propagate through neu- 
tral hydrogen because it is immediately absorbed. There- 
fore, quasar spectra released before the reionisation redshift 
must be completely absorbed blueward of the Lyman-ff emis- 
sion line. The appearance of this so-called Gunn-Peterson ef- 
fect (Gunn and Peterson 1965[ ) at high redshift thus signals 
the transition from ionised into neutral material. Using this 
technique, the reionisation was estimated to end at redshifts 
~ 6 . . . 8 ( Fan et al. 2006 ), while the secondary polarisation of 
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the CMB implies a reionisation redshift of 1 1 .0 ± 1 .4 if instan- 
taneous reionisation is assumed ( |Dunkley et a l. 2008 ). These 
apparently discrepant values do not contradict each other be- 
cause a small admixture of neutral hydrogen is enough to pro- 
duce the Gunn-Peterson effect, which therefore persists until 
reionisation has completed. 



C. Cosmological weak lensing 

Compared to the outlined procedure to obtain <x 8 from the 
CMB, it appears completely straightforward to derive it from 
the cosmic-shear measurements. As we have discussed in 



( VIII. B. 3 i, the cosmic-shear power spectrum is proportional 



to times the amplitude A of the dark-matter power spec- 
trum, which leads to the approximate degeneracy Q.°^ Q cr% » 
const, between cr 8 and the matter-density parameter Q m o- 

A more subtle dependence on Q m o and to some degree also 
on other cosmological parameters is introd uced by the geo- 
metrical weight function W(w', w) shown in (231 1, and by the 



growth of the power spectrum along the line-of-sight. This 
slightly modifies the form of the cr 8 -£2 m o degeneracy, but does 
not lift it. However, knowing £2 m o well enough, we should 
be able to read cr 8 off the cosmic-shear correlation function. 
However, there are three problems associated with that. 

First, the cosmic shear measured on angular scales below 
~ 10' is heavily influenced by the onset of non-linear structure 
growth and the effect this has on the dark-matter power spec- 
trum. While the linear growth factor can be straightforwardly 
calculated analytically, non-linear growth can only be quan- 
tified by means of large numerical simulations and recipes 
derived from them (Pea cock and Dodds] [T99 6; Smith^FaT] 
2003b). Insufficient knowledge of the non-linear dark-matter 
power spectrum is a major uncertainty in the cosmological in- 
terpretation of cosmic shear. 

Second, the amplitude of cosmological weak-lensing ef- 
fects depends on the redshift distribution of the sources used 
for measuring ellipticities. Since these background galaxies 
are typically very faint, it is demanding to measure their red- 
shifts. Two methods have typically been used. One adapts the 
known redshift distribution of sources in narrow, very deep 
observations such as the Hubble Deep Field to the character- 
istics of the observation to be analysed. The other relies on 
photometric redshifts, i.e. redshift estimates based on multi- 
band photometry. Yet, the precise redshift distribution of the 
background sources adds additional uncertainty to estimates 
of cr 8 . 

Third, it is possible that systematic effects remain in weak- 
lensing measurements because the effect is so small, and many 
corrections have to be applied to measured ellipticities before 
the cosmic shear can be extracted. Advanced correction meth- 
ods have been developed which made the B-mode contamina- 
tion almost or completely disappear. This is good news, but it 
does not yet guarantee the absence of other systematic effects 
in the data. 

Nonetheless, cosmic lensing, combined with estimates of 
the matter-density parameter, is perhaps the most promising 
method for precisely determining <x 8 . Table \H\ lists values of 



<x 8 derived from some cosmic-shear measurements under the 
assumption of D. m o = 0.30 in a spatially-flat universe. 



0.86*® 
0.72 ± 0.09 



0.97: 
1.02: 
0.83: 



0.13 
0.16 
0.07 



0.68 ±0.13 
0.85 ± 0.06 
0.80 ±0.1 
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0.74 ± 0.04 
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CFHTLS-wide 



reference 
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IJarvis et al.\ 2003 
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(Bacon et al. 2003 



(Rhodes et al.\ 2004 1 



I Van Waerbeke et al.\ 
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jHeymans et al. 2005 



iHoekstra era/. ,2006 
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2007) 



]Massey et al. 2007b) 



Benjamin et al. 2007| 
JFuer a/. | [2008) 



TABLE II Values for tr 8 derived from cosmic-shear measurements 
under the assumption of a spatially-flat universe with Q. mQ = 0.3. 

The systematic deviations between them are probably due 
to three different uncertainties. First, the redshift distribution 
of the background sources needs to be known; second, non- 
linear evolution of the matter power spectrum need to be ac- 
curately modelled; and third, all kinds of PSF distortions need 
to be corrected. 



D. Galaxy clusters 

1 . The mass function 

Based on the assumption that the density contrast is a Gaus- 



sian random field and the spherical -collapse model, ( Press and 



|Schechter] |1974) derived a mass function for dark-matter ha 
los. It compares the standard deviation cr s of the density- 
fluctuation field to the linear density-contrast threshold 6 C « 
1 .686 for collapse in the spherical-collapse model. The mean 
mass contained in spheres of radius R sets the halo mass, 
which brings the mean (dark-) matter density p into the game. 

The standard deviation cr R is related to the power spectrum. 
For convenience, we introduce an effective slope 

dlnP(fc) 

n = (274) 
din A: 

for the power spectrum, which will of course be scale- 
dependent. On large scales, n * 1, while n — > -3 on small 
scales, i.e. for small halo masses. For galaxy clusters, n « — 1. 

We introduce the non-linear mass scale M, as the mass con- 
tained in spheres of radius R chosen such that cr R = 1 . Since 
cr R grows with the linear growth factor D + {a), the non-linear 
mass grows with time. It is convenient here to express the am- 
plitude of the power spectrum, and thus <x 8 , in terms of M*. In 
sufficient approximation, 

'M, s 



o-R 
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In terms of the dimensionless mass m = M/M*, the Press- 
Schechter mass function can be cast into the form 



N(m, a)dm = \ — 



P<5 C 



nMiD + (a) 



am" 2 exp I - - 



2» 



dm 



2D\(d) 

(276) 

The Press-Schechter mass function, and some improved 
variants of it (TJenkins et oT] |2001| |Sheth and Tormeri] |2002| l, 
have been spectacularly confirmed by numerical simulations 
(Springel et al. 2005 i. It shows that the mass function is 



a power law with an exponential cut-off near the non-linear 
mass scale M,. For galaxy clusters, n at -1, thus a as 1/3, and 



N(m, a)dm oc m 5 ^ 3 exp 



2D 2 + (a) 



,2/3 



dm . 



(277) 



with an amplitude characterised by M„ the mean dark-matter 
density p, and the growth factor D+{a). 

This opens a way to constrain cosmological parameters as 
well as <t 8 with galaxy clusters: if the abundance and evolu- 
tion of the cluster mass function can be measured, they can 
be determined from the mass scale of the exponential cut-off 
and the amplitude of the power-law end. Today, the non-linear 
mass scale is a few times 10 13 M Q . Therefore, the exponential 
cut-off in the halo mass will not be seen in the galaxy mass 
function. Clusters, however, show a pronounced exponential 
cut-off ( |Rosati et al.\ 2002} , and thus their population is very 
sensitive to changes in cr g . In principle, therefore, crg should 
be very well constrained by the cluster population. 



2. What is a cluster's mass? 

The main problem here is how observable cluster proper- 
ties should be related to quantities used in theory. Cluster 
masses, as used in the theoretical mass function ( |277| i, are not 
observables. Cluster observables are the X-ray temperature 
and flux, the optical luminosity and the velocity distribution 
of their galaxies, and their gravitational-lensing effects. Be- 
fore we discuss their relation to mass, let us first see what the 
"mass of a galaxy cluster" could be. 

It is easy to define masses of gravitationally bound, well 
localised objects, such as planets or stars. They have a well- 
defined boundary, e.g. the planetary surfaces or the stellar 
photospheres. This is markedly different for objects like 
galaxies and galaxy clusters. As far as we know, their den- 
sities drop smoothly towards zero like power laws, oc r -( 2 - 3 ) 
Thus, although they are gravitationally bound, it is less ob- 
vious what should be seen as their outer boundary. Strictly 
speaking, there is none. 

The only way out is then to define an outer boundary in 
such a way that it is well-defined in theory and identifiable 
in observational data. A common choice was introduced in 
it defines the boundary by the mean overdensity it 
Although this is problematic as well, it may be as 



§ V.A.2 



encloses. 

good as it gets. Three immediately obvious problems created 
by this definition are that objects like galaxy clusters are often 
irregularly shaped rather than spherical, that the overdensity 
of 200 is as arbitrary as any other, even if it is inspired by 



virial equilibrium in the spherical-collapse model, and that its 
measurement requires a sufficiently accurate density profile to 
be known or assumed. 

How could standardised radii such as /?20o be measured? 
This could for instance be achieved applying equations such 



as ( 147 1 after measuring the slope ft and the core radius of 
the X-ray surface brightness profile together with the X-ray 
temperature, by calibrating an assumed density profile with 
galaxy kinematics based on the virial theorem, or by con- 
straining the cluster mass profile with gravitational lensing. 

Obviously, all these measurements have their own prob- 
lems. Being sensitive to all mass along the line-of-sight, grav- 
itational lensing cannot distinguish between mass bound to a 
cluster or just projected onto it. Any measurement based on 
the virial theorem must of course rely on virial equilibrium, 
which takes time to be established and is often perturbed in 
real clusters because of merging and accretion. The common 
interpretation of X-ray measurements requires the assumption 
that the X-ray gas be in hydrostatic equilibrium with the host 
cluster's gravitational potential. 

This illustrates that it may be fair to say that there is no such 
thing as the mass of a galaxy cluster. Even if measurements of 
cluster "radii" were less dubious, it remained unclear whether 
they mean the same as those assumed in theory, which are re- 
lated to the spherical-collapse model. Interestingly, but not 
surprisingly, cluster masses obtained from numerical simula- 
tions suffer from the same poor definition of the concept of a 
"cluster radius". 

How can we make progress then? Observables such as the 
cluster temperature Tx or its X-ray luminosity Lx should be 
related to the depth of the gravitational potential they are em- 
bedded in, which should in turn be related to some measure of 
the total mass. If we can calibrate such expected temperature- 
mass or luminosity-mass relations, e.g. using numerical simu- 
lations of galaxy clusters, a direct comparison between theory 
and observations seems possible. This is sometimes called an 
external calibration of the required relations. 

Internal or self-calibrations, i.e. calibrations based on clus- 
ter data alone, have become increasingly fashionable over the 
past years. Here, empirical temperature-mass and luminosity- 
mass relations are obtained based on one or more of the mass 
estimates sketched above. 
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(Smith ef al. I 2003a I 
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IPierpaoli et al. 


2003 



Dahle 2006 



TABLE III Values of erg derived from the galaxy-cluster population 
based on different observational data. 
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The result of both procedures is qualitatively the same. It 
allows the conversion of observables to mass, and thus of the 
observed cluster temperature or luminosity functions to mass 
functions, which can then be compared to theory. The shape 
and amplitude of the power spectrum and the growth factor 
can then be adapted to optimise the agreement between ob- 
served and expected mass functions. Clusters at moderate or 
high redshift constrain the evolution of the mass function and 
allow an independent estimate of the matter-density parameter 
Sl m o, as sketched in § |V.C| before. 

In view of the many difficulties listed, it is an astonishing 
fact that, when applied to cluster samples rather than individ- 
ual clusters, the determination of the cluster mass function 
and its evolution seems to work quite well. Values for cr 8 



derived therefrom are given in Tab. Ill Systematic differences 
between these values are most likely due to the uncertainties 
of the calibration procedures applied to the relations between 
cluster masses and observables. 



E. The Lyman-** forest 

The Lyman-a forest lines mentioned before arise from ab- 
sorption in neutral-hydrogen clouds. It is reasonable to as- 
sume that they are located where the dark matter is overdense, 
such that fluctuations in the neutral-hydrogen density are pro- 
portional to the dark-matter density contrast 5. Then, obser- 
vations of the Lyman-of forest lines, in particular their number 
per unit redshift, equivalent-width distribution and correlation 
properties, must contain information on the shape and the nor- 
malisation of the dark-matter power spectrum. 

Retrieving this information would be straightforward if 
the biasing relation between the neutral-hydrogen and dark- 
matter densities was known and simple, the gas temperature 
was known, redshifts were unaffected by peculiar motions, 
and non-linear structure evolution could be neglected. While 
this is not the case in reality, these perturbing effects can be 
calibrated with numerical simulations and corrected. Two dif- 
ferent methods have been proposed; one inverts the Lyman-Q' 
forest directly ( |Nusser and Haehnelt 1999 ), while the other 
adapts power spectra of numerical simulations to those recov- 



ered from Lyman-o- forest observations (Croft et al. 1998 1. 
Both methods have allowed interesting constraints which are 
summarised in Tab. HVl 
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data 
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0.73 ± 0.04 
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Keck spectra 
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l| Viel and Haehnelt| 
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TABLE IV Selection of erg values obtained from Lyman-o- forest 
data alone. 



The method is young, but promising. Values of erg obtained 
so far seem to be typically systematically higher than those de- 
rived from other observables. Possible sources of systematic 



error are the onset of nonlinear structure growth and, perhaps 
most importantly, the temperature of the hydrogen gas and its 
relation to the gas density (often called its equation of state). 



XI. INFLATION AND DARK ENERGY 
A. Cosmological inflation 

1 . Motivation 

In the preceding chapters, we have seen the remarkable 
success of the cosmological standard model, which is built 
upon the two symmetry assumptions underlying the class of 
Friedmann-Lemaitre-Robertson- Walker models which expe- 
rienced a Big Bang a finite time ago. We shall now discuss a 
fundamental problem of these models, and a possible way out. 
Historically, the problem was raised in a different way, but it 
intrudes with the very straightforward realisation that it is by 
no means obvious why the CMB should appear as isotropic as 
it is, and why there should be large coherent structures in it. 

Let us begin with the so-called comoving particle horizon, 
which is the comoving distance that light can travel between 
the Big Bang and time t. It is given by \\15\ with the sound 
speed c s replaced by the light speed c, 



cdt 

a 



■ ( VT+a - yfa) = 282.8 Mpc 



(278) 

with a m 0.33 as defined below ( |164) , On the other hand, we 
have seen in (176i that the comoving angular-diameter dis- 
tance to the CMB is 



c I Q. m (\ 
w(a Iec ) = 3.195— 

Hq \ <j.j 



-0.4 



Gpc , 



(279) 



which implies that the angular size of the particle horizon is 

,-0.1 



w(a KC ) 



1.14° 



~03~ 



(280) 



The physical meaning of the particle horizon is that no 
event between the Big Bang and recombination can exert 
any influence on a given particle if it is more than the hori- 
zon length away. Our simple calculation thus shows that we 
can understand how causal processes could establish identi- 
cal physical conditions in patches of the sky with about one 
degree radius. Points on the CMB separated by larger angles 
were never causally connected before the CMB was released. 
It is therefore not at all plausible how the CMB could have at- 
tained almost the same temperature across the entire sky! The 
simple fact that the CMB is almost entirely isotropic across 
the sky thus poses a problem which the standard cosmological 
model is apparently unable to solve. Moreover, the formation 
of coherent structures larger than the particle horizon remains 
mysterious. This is one way to state the horizon problem. It 
is sometimes called the causality problem: How can coherent 
structures in the CMB be larger than the particle horizon was 
at recombination? 
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Another uncomfortable problem of the standard cosmologi- 
cal model is the flatness, or at least the near-flatness, of spatial 
hypersurfaces of our Universe. To see this, we write Fried- 
mann's equation ([6]) in the form 



H 2 (a) = H 2 (a) 



^total( fl ) ' 



Kc 2 



2 H 2 



which is equivalent to 



^total(fl) - 1 = 



Kc 1 
a 2 H 2 



(281) 



(282) 



The right-hand side typically grows as some power of the time 
t, 



{a 2 H 2 y l = a 2 oc t" ; 



(283) 



e.g. j3 = 2/3 in a matter-dominated universe without cosmo- 
logical constant. Then, 



^total(fl) 1 OC £ 



(284) 



This shows that any deviation of the total density parameter 
f^totai from unity tends to grow with time. Thus, (spatial) flat- 
ness is an unstable property. If it is not very precisely flat in 
the beginning, the Universe will develop away from flatness. 
Since we know that spatial hypersurfaces are now almost flat, 
l^totai( a ) _ 1| S5 1% (cf. Tab. [I]), the deviation from flatness 
must have been of order 



|£W(<»rec) - 1] < 1%X 



10" 



(285) 



or ten parts per million at the time of recombination. Clearly, 
this requires enormous fine-tuning. This is called the flatness 
problem: How can we understand flatness in the late universe 
without assuming an extreme degree of fine-tuning at early 
times? 



2. The idea of inflation 



over, it starts expanding again, but if the reduction was suf- 
ficiently large, it could remain within the causally-connected 
region at least until the present. 

How could such a shrinking comoving Hubble radius be 
arranged? Obviously, we require 



d c 
At aH 



(aH) 



-{aH + aH) 



(aH) 2 



<0 



(287) 



which is possible if and only if a > 0, in other words, if the 
expansion of the Universe accelerates. This appears counter- 
intuitive because the cosmic expansion is dominated by grav- 
ity, which should be attractive and thus necessarily decelerate 
the expansion. Friedmann's equation (j7]i implies the matter 
condition 



pc 2 + 3p<0 



pc^ 
P< ~- 



(288) 



In other words, cosmic acceleration is possible if and only if 
the dominant ingredient of the cosmic fluid has sufficiently 
negative pressure. Equation ([8]) implies the density evolution 



p a p \ a 

3- 1 + -^ <-2- 

p a\ pc z j a 



(289) 



for a cosmic fluid satisfying (288i, showing that its density 
would fall off flatter than a~ 2 , and thus less steeply than the 
matter or radiation densities. Thus, once a component of 
the cosmic fluid with sufficiently negative pressure reaches a 
density comparable to the densities of matter or radiation, it 
quickly starts dominating the cosmic expansion. In the limit- 
ing case p — or p — -pc 2 , ^ reduces to 



a 2,7tG 

- - —r—p - const. > , 
a 3 



which implies the exponential expansion or inflation 

v l/2 

a oc exp 



(SttG 



(290) 



(291) 



Since c/H is the Hubble radius, the quantity r H = c/(aH) is 3 - Slow ro11 ' structure formation, and observational constraints 



the comoving Hubble radius. According to (283 1, it ty pically 
grows with time as r H oc f^/ 2 . Since we can write (282 1 as 



^totai(a) - 1 = Kr H 



(286) 



this is equivalent to the flatness problem. 

This motivates the idea that at least the flatness problem 
may be solved if the comoving Hubble radius could, at least 
for some sufficiently long period, shrink with time. If that 
could be arranged, any deviation of Q tota i(a) from unity would 
be driven towards zero. 

Conveniently, such an arrangement would also remove or 
at least alleviate the causality problem. Possibly, the Hubble 
radius, characterising the radius of the observable universe, 
could be driven inside the horizon and thus move the entire ob- 
servable universe into a causally-connected region. When the 
hypothesised epoch of a shrinking comoving Hubble radius is 



We have seen that we need inflation to solve the flatness and 
causality problems, and inflation needs a form of matter with 
negative pressure. What could that be? Consider a scalar field 
<p with a self-interaction potential V((p). Then, field theory 
shows that pressure and density of the scalar field are related 
by the equation of state 



p,f, - wp^c 2 with w = 



V 



V 



(292) 



Evidently, negative pressure is possible if the kinetic energy of 
the scalar field is sufficiently smaller than its potential energy. 
For the cosmological-constant case, (f> — 0, we have w = — 1 
or p = -pc 2 , in agreement with the conclusion from ( 289 ». 
In other words, a suitably strongly self-interacting scalar field 
has exactly the properties we need. 
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Inflation, i.e. accelerated expansion, broadly requires 2 to 
be sufficiently smaller than V. Moreover, we need inflation 
to operate long enough to drive the total matter density pa- 
rameter sufficiently close to unity for it to remain there to the 
present day. These two conditions are conventionally cast into 
the form 



1 



24nG \ V 



V 



1 V" 

«: 1 and n = «: 1 

' 8nG V 



(293) 



(see the standard textbook on inflation, (Liddl e~and Lyth| 
2000)). They are called the slow-roll conditions. The first 



assures that inflation can set in, because if it is satisfied, the 
potential has a small gradient and cannot drive rapid changes 
of the scalar field. The second restricts the curvature of the 
potential and thus assures that the inflationary condition is sat- 
isfied long enough. 

Estimates show that inflation needs to expand the Universe 
by ~ 50 ... 60 e-foldings ( |Alabidi and Lyfh[ |2006) (i.e. by a 
factor of e 50 60 ) for it to solve the causality and flatness prob- 
lems. Inflation ends once the slow-roll conditions are vio- 
lated. By then, the Universe will have become extremely cold. 
While the density of the inflaton field <f> will be approximately 
the same as at the onset of inflation (as for the cosmological 
constant, this is a consequence of the negative pressure), all 
other matter and radiation fields will have their energy densi- 
ties lowered by factors of a" 3 - 4 , i.e. by < 100 orders of mag- 
nitude. 

Once e approaches unity, the kinetic term cf 1 will dominate 
the potential, and the scalar field will start oscillating rapidly. 
It is assumed that the scalar field then decays into ordinary 
matter which fills or reheats the Universe after inflation is 
over. 

It is an extremely interesting aspect of inflation that it 
also provides a mechanism for seeding structure formation 



(Mukhanov and Chibisov 198 1| >. As any other quantum field, 
the inflaton field <f> must have undergone vacuum oscillations 
because the zero-point energy of a quantum harmonic oscil- 
lator cannot vanish due to Heisenberg's uncertainty principle. 
These vacuum oscillations cause the spontaneous creation and 
annihilation of particle-antiparticle pairs. Once inflation sets 
in, vacuum fluctuation modes are quickly driven out of the 
horizon and lose causal connection. Then, they cannot de- 
cay any more and "freeze in". Thus, inflation introduces the 
breath-taking notion that density fluctuations in our Universe 
today may have been seeded by vacuum fluctuations of the 
inflaton field before inflation set in and enlarged them to cos- 
mological scales. 



This idea has precisely quantifiable consequences (Bran- 


denberger| 
1985pigf 


1984(|GufhandPi| 1982 
;ls 1983). First, linear de 

ii r 


Halliwell and Hawking 


nsity fluctuations created 



by the inflaton field must be adiabatic and, by the central limit 
theorem, they should form a Gaussian random field. This is 
because they arise from incoherent superposition of extremely 
many independent fluctuation modes whose amplitude and 
wave number are all drawn from the same probability distri- 
bution. Under these circumstances, the central limit theorem 
shows that the result, i.e. the superposition of all these modes, 
must be a Gaussian random field. 



Second, it implies that the statistics of density fluctuations 
in the Universe today must be explicable by the statistics of 
vacuum fluctuations in a scalar quantum field. This is in- 
deed the case. The power spectrum resulting from this con- 
sideration is very close to the scale-free Harrison-Zel'dovich- 
Peebles shape introduced in § |I.C.2| 

P s (k) cc k" s , (294) 

with « s w 1 . The spectral index n s would be precisely unity 
if inflation lasted forever. Since this was obviously not so, 
n s must deviate slightly from unity, and detailed calculations 



show that it must be slightly smaller ( Alabidi and Lyth 2006 



Liddle and Lyth 2000 1, 



n s — 1 + 2/7 - 6e . 



(295) 



The latest WMAP measurements (Komatsu et al. , 2008) do in 
fact show that 



tu = 0.960 



+0.014 
-0.013 



(296) 



cf. Tab. [j] The completely scale-invariant spectrum, n s = 1, is 
thus excluded at more than 3cr. The measured deviation of n s 
from unity also restricts the number of e-foldings completed 
by inflation. Under fairly general assumptions, 



N = 54 + 7 



(297) 



( Alabidi and Lyth , 2006 ) based on the three-year WMAP data. 

Another prediction of inflation is that it may excite not only 
scalar, but also tensor perturbations ( |Liddle and Lyth||2000| l. 
Scalar perturbations lead to the density fluctuations, tensor 
perturbations correspond to gravitational waves. Vector per- 
turbations do not play any role because they decay quickly as 
the universe expands. Simple models of inflation predict that 
the ratio r between the amplitudes of tensor and scalar pertur- 
bations, taken in the limit of small wave numbers, is 



16e . 



(298) 



An inflationary background of gravitational waves is in prin- 



ciple detectable through the polarisation of the CMB (Seljak 



and Zaldarriaga 1997| l. Limits of order r < 0.05 are expected 



from the upcoming Planck satellite ( Kamionkowski and Jaffe 
2001 Kinney 1998 ). Together with the result n s + 1 from 



WMAP, it will then be possible to constrain viable inflation 
models, i.e. to constrain the shape of the inflaton potential. 



B. Dark energy 

1 . Motivation 

The CMB shows us that the Universe is at least nearly spa- 
tially flat. Constraints from kinematics, from cluster evolution 
and from the CMB show that the matter density alone cannot 
be responsible for flattening space, and primordial nucleosyn- 
thesis and the CMB show that baryons contribute at a very 
low level only. Something is missing, and it even dominates 
today's cosmic fluid. 
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From structure formation, we know that this remaining con- 
stituent cannot clump on the scales covered by the galaxy sur- 
veys and below. It is thus different from dark matter. The 
terms dark energy, kosmon or quintessence have been coined 
for it ( [Peebles and Ratra| [2003] [19881 |Wetterich| [19881 The 
type-la supernovae tell us that it behaves at least very similar 
to a cosmological constant in the recent cosmic past. Maybe 
the dark energy is a cosmological constant? Nothing cur- 
rently indicates any deviation from this "simplest" assump- 
tion (e.g. Kom atsu et al. |2 008 ). So far, the cosmological con- 
stant is a perfectly viable description for all observational ev- 
idence we have. However, this is deeply dissatisfactory from 
the point of view of theoretical physics. The problem is the 
value of Qaq. As we have seen above, a self-interacting scalar 
field with negligible kinetic energy behaves like a cosmolog- 
ical constant. Then, its density should simply be given by its 
potential V. Simple arguments suggest that V should be the 
Planck mass divided by the third power of the Planck length, 
which turns out to be 120 orders of magnitude larger than the 
cosmological constant derived from observations. Since this 
fails, it seems natural to expect that the cosmological constant 
should vanish, but it does not. The main problem with the 
cosmological constant is therefore, why is it not zero if it is so 
small? 

Anthropic arguments suggest that we observe the Universe 
as it is because it needs to be that way for us to exist. The 
Universe needs to be old enough for stars to have produced 
carbon, on which life must be based to the best of our knowl- 
edge. Yet, it must not be too old for the majority of stars 
to have evolved past their main-sequence life for the stabil- 
ity of benign conditions for life. If the cosmological constant 
was much larger, inflationary expansion would have prohib- 
ited structure formation and thus also the formation of life. It 
may thus be that the cosmological parameters have the values 
we measure because otherwise we would not exist to measure 
them. 

The explanation of inflation by means of an inflaton field 
suggests another way out. As we have seen there, accelerated 
expansion can be driven by a self-interacting scalar field while 
its potential energy dominates. Moreover, it can be shown that 
if the potential V has an appropriate shape, the dark energy has 
attractor properties in the sense that a vast range of initial den- 



sity values can evolve towards the same value today (Zlatev 



\et al.\ |1999| >. At this point at the latest, we leave the realm 
of what may be considered the established cosmological stan- 
dard model, and enter into the very controversial discussion 
of dark energy. 



2. Observational constraints? 

If the dark energy is indeed dynamical and provided by a 
self-interacting scalar field, how can we find out more about 
it? Reviewing the cosmological measurements we have dis- 
cussed so far, it becomes evident that they are all derived from 
constraints on 

• cosmic time, as in the age of the Galaxy or of globular 
clusters, or in primordial nucleosynthesis; 



• distances, as in the spatial flatness derived from the 
CMB, the type-la supernovae or the geometry of cos- 
mological weak lensing; or 

• the growth of cosmic structures, as in the acoustic os- 
cillations in the CMB, the evolution of the cluster pop- 
ulation, the structures in the galaxy distribution or the 
source of cosmological weak-lensing effects. 

We must therefore seek to constrain the dark energy by 
measurements of distances, times, and structure growth. Since 
they can all be traced back to the expansion behaviour of the 
universe as described by Friedmann's equation, we must see 
how the dark energy enters there, and what effects it can seed 
through it. 

Let us therefore assume that the dark energy is a suitably 
self-interacting, homogeneous scalar field. Then, its pressure 
can be described by 



p = w(a)pc 



(299) 



where the equation-of-state parameter w is some function of a. 



According to ( |288| l, accelerated expansion needs w < -1/3, 
and the cosmological constant corresponds to w = — 1, Since 
all cosmological measurements to date are in agreement with 
the assumption of a cosmological constant, we need to ar- 
range things such that w — > -1 today. Suppose we had some 
function w(a), which could either be obtained from a phe- 
nomenological choice, a model for the self-interaction poten- 



sation. Then, ( [289 



tial V((/>) through ( 292 1 or from a suitable ad-hoc parameteri- 



implies 



p a 
- = -3(1 +w)- 
p a 



(300) 



or 



p(a) = po exp j-3 J [1 + J = Po/(a) ■ (301) 

For constant w, this simplifies to 

p(a) = p ( ,exp [-3(1 + ic)lna] = p a' 3(l+w) ■ (302) 

If w = -1, we recover the cosmological-constant case p = 
po = const., for pressure-less material, w = and p cc cr 3 , 
and for radiation, w = 1/3 and p cc a~ 4 . 

Therefore, we can take account of the dynamical dark en- 



ergy by replacing the term Qao in the Friedmann equation ( 10 1 
by Q.£,Eof(a), and the expansion function E(a) turns into 

E(a) = [fi lfl fl- 4 + fi m0 <r 3 + Odbo/Ca) + "ko<T 2 ] 1/2 , (303) 

where Qko = 1 - fi r o - ^ m o - ^deo is a density parameter 
assigned to the spatial curvature. 

We thus see that the equation-of-state parameter enters the 
expansion function in integrated form. Since all cosmological 
observables are integrals over the expansion function, includ- 
ing the growth factor D + (a), this implies that cosmological 
observables measure integrals over the integrated equation-of- 
state function w(a). Note that in cosmological models with 
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dynamical dark energy, the growth factor is more compli- 
cated than described by Eq. (|23]l because dark-energy clus- 



tering needs to be taken into account (e.g. Kunz and Sapone 



2007 1. Needless to say, the dependence of cosmological mea- 



surements on the exact form of w(a) will be extremely weak, 
which in turn implies that extremely accurate measurements 
will be necessary for constraining the nature of the dark en- 
ergy. 



of sources in different redshift shells, the growth factor can be 
probed differentially. First examples for this technique have 
been published. They give rise to the expectation that clarify- 
ing the nature of the dark energy may indeed be feasible in the 
near future. 
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FIG. 12 Relative changes of the growth factor D + and the luminosity 
distance Di ura with w as functions of redshift, assuming constant w. 

In order to illustrate the required accuracies, let us consider 
by how much the angular-diameter distance and the growth 
factor change compared to ACDM upon changes in w away 
from -1, 

dlnDang(z) dlnD + (z) 
dw dw 

as a function of redshift z. Assuming Q. m Q = 0.3 and Qao 
0.7. we find typical values between -0.1 and -0.2 at most. 
Since we currently expect deviations of w from -1 at most 
at the ~ 10% level, accurate constraints on the dark energy 
require relative accuracies of distances and the growth factor 
at the per-cent level. 

Clearly, all models for dynamical dark energy proposed so 
far are purely phenomenological. The most straightforward 
question to be addressed in this situation may be whether a 
cosmological constant can be ruled out in some way. All suit- 
able cosmological information will need to be combined in 
order to progress. 

Currently, the largest hope is put on the BAOs (see 
§ VII.A.6| l and on so-called tomographic measurements. 
BAOs have a characteristic physical scale whose angular scale 
depends on the angular-diameter distance, which is in turn 
sensitive to the dark energy. Tomography attempts to trace 
the evolution of structures throughout cosmic history ( |Hu| 
2002a|b| l. An example is given by weak gravitational lensing: 
Since its geometrical sensitivity peaks approximately half- 
way between the sources and the observer, sources at higher 
redshift also probe more distant, and thus less evolved, cosmic 
structures. If lensing effects can be measured for sub-samples 



I am most grateful to many colleagues for inspiring and 
clarifying discussions, in particular to Peter Schneider, Achim 
WeiB and Simon White, whose detailed comments helped im- 
proving this review substantially. Three anonymous refer- 
ees contributed numerous thoughtful and constructive com- 
ments. This work was supported in part by the German Sci- 
ence Foundation (DFG) through the Collaborative Research 
Centres SFB 439 and TRR 33. 



References 

Abell, G. O., J. H. G. Corwin, and R. P. Olowin, 1989, Astrophys. J. 

Supp. 70, 1. 
Aguirre, A., 1999, Astrophys. J. 525, 583. 

Alabidi, L., and D. H. Lyth, 2006, Journal of Cosmology and Astro- 
Particle Physics 8, 13. 
Alcock, C, R. A. Allsman, D. R. Alves, T. S. Axelrod, A. C. Becker, 

D. P. Bennett, K. H. Cook, N. Dalai, A. J. Drake, K. C. Freeman, 
M. Geha, K. Griest, et al, 2000, Astrophys. J. 542, 281. 

Allen, S. W„ R. W. Schmidt, A. C. Fabian, and H. Ebeling, 2003, 
Mon. Not. Royal Astron. Soc. 342, 287. 

Alpher, R. A., and R. C. Herman, 1949, Physical Review 75, 1089. 

Asplund, M., D. L. Lambert, P. E. Nissen, F. Primas, and V. V. Smith, 
2006, Astrophys. J. 644, 229. 

Astier, P., J. Guy, N. Regnault, R. Pain, E. Aubourg, D. Balam, 
S. Basa, R. G. Carlberg, S. Fabbro, D. Fouchez, I. M. Hook, D. A. 
Howell, etai, 2006, Astron. Astrophys. 447, 31. 

Baade, W., 1926, Astronomische Nachrichten 228, 359. 

Baade, W., 1956, Pub. Astron. Soc. Pacific 68, 5. 

Bacon, D. J., R. J. Massey, A. R. Refregier, and R. S. Ellis, 2003, 
Mon. Not. Royal Astron. Soc. 344, 673. 

Barnes, J., and P. Hut, 1986, Nature 324, 446. 

Bartelmann, M., and P. Schneider, 2001, Phys. Rep. 340, 291. 

Benjamin, J., C. Heymans, E. Semboloni, L. van Waerbeke, H. Hoek- 
stra, T. Erben, M. D. Gladders, M. Hetterscheidt, Y. Mellier, and 
H. K. C. Yee, 2007, Mon. Not. Royal Astron. Soc. 381, 702. 

Blandford, R., and R. Narayan, 1986, Astrophys. J. 310, 568. 

Blinnikov, S. I., F. K. Ropke, E. I. Sorokina, M. Gieseler, M. Rei- 
necke, C. Travaglio, W. Hillebrandt, and M. Stritzinger, 2006, As- 
tron. Astrophys. 453, 229. 

Boggess, N. W., J. C. Mather, R. Weiss, C. L. Bennett, E. S. Cheng, 

E. Dwek, S. Gulkis, M. G. Hauser, M. A. Janssen, T. Kelsall, S. S. 
Meyer, S. H. Moseley, etai, 1992, Astrophys. J. 397, 420. 

Bond, J. R., and G. Efstathiou, 1984, Astrophys. J. Lett. 285, L45. 

Brandenberger, R. H, 1984, Nuclear Physics B 245, 328. 

Brown, M. L., A. N. Taylor, D. J. Bacon, M. E. Gray, S. Dye, 

K. Meisenheimer, and C. Wolf, 2003, Mon. Not. Royal Astron. 

Soc. 341, 100. 

Cayrel, R., V. Hill, T. C. Beers, B. Barbuy, M. Spite, F. Spite, B. Plez, 
J. Andersen, P. Bonifacio, P. Francois, P. Molaro, B. Nordstrom, 
et al, 2001, Nature 409, 691. 

Chaboyer, B., 1998, Phys. Rep. 307, 23. 



48 



Chaboyer, B., P. Demarque, P. J. Kernan, and L. M. Krauss, 1998, 

Astrophys. J. 494, 96. 
Chandrasekhar, S., 1984, Reviews of Modern Physics 56, 137. 
Charbonnel, C, and F. Primas, 2005, Astron. Astrophys. 442, 961. 
Ciardi, B., and A. Ferrara, 2005, Space Science Reviews 116, 625. 
Cole, S., W. J. Percival, J. A. Peacock, P. Norberg, C. M. Baugh, 

C. S. Frenk, I. Baldry, J. Bland-Hawthorn, T. Bridges, R. Cannon, 

M. Colless, C. Collins, etal, 2005, Mon. Not. Royal Astron. Soc. 

362, 505. 

Coles, P., and F. Lucchin, 2002, Cosmology: The Origin and Evolu- 
tion of Cosmic Structure, Second Edition (Cosmology: The Ori- 
gin and Evolution of Cosmic Structure, Second Edition, by Peter 
Coles, Francesco Lucchin, pp. 512. ISBN 0-471-48909-3. Wiley- 
VCH , July 2002.). 

Colgate, S. A., and C. McKee, 1969, Astrophys. J. 157, 623. 

Colless, M., 1999, Royal Society of London Philosophical Transac- 
tions Series A 357, 105. 

Cool, A. M., G. Piotto, and I. R. King, 1996, Astrophys. J. 468, 655. 

Cowan, J. J., F.-K. Thielemann, and J. W. Truran, 1987, Astrophys. 
J. 323, 543. 

Cowan, J. J., F.-K. Thielemann, and J. W. Truran, 1991, Phys. Rep. 
208, 267. 

Croft, R. A. C, D. H. Weinberg, M. Bolte, S. Buries, L. Hernquist, 
N. Katz, D. Kirkman, and D. Tytler, 2002, Astrophys. J. 581, 20. 

Croft, R. A. C, D. H. Weinberg, N. Katz, and L. Hernquist, 1998, 
Astrophys. J. 495, 44. 

Cypriano, E. S., L. J. Sodre, J. -P. Kneib, and L. E. Campusano, 2004, 
Astrophys. J. 613, 95. 

Dahle, H, 2006, Astrophys. J. 653, 954. 

Dame, T. M., D. Hartmann, and P. Thaddeus, 2001, Astrophys. J. 
547, 792. 

Davis, M., G. Efstathiou, C. S. Frenk, and S. D. M. White, 1985, 

Astrophys. J. 292, 371. 
Djorgovski, S., and M. Davis, 1987, Astrophys. J. 313, 59. 
Dodelson, S., 2003, Modern cosmology (Modern cosmology / Scott 

Dodelson. Amsterdam (Netherlands): Academic Press. ISBN 0- 

12-219141-2, 2003, XIII + 440 p.). 
Doroshkevich, A. G., 1970, Astrofizika 6, 581. 
Dressier, A., D. Lynden-Bell, D. Burstein, R. L. Davies, S. M. Faber, 

R. Terlevich, and G. Wegner, 1987, Astrophys. J. 313, 42. 
Dunkley, L, E. Komatsu, M. R. Nolta, D. N. Spergel, D. Larson, 

G. Hinshaw, L. Page, C. L. Bennett, B. Gold, N. Jarosik, J. L. 

Weiland, M. Halpern, et al, 2008, ArXiv e-prints 803. 
Durrer, R., 2008, The cosmic microwave background (Cambridge 

University Press). 
Efstathiou, G., M. Davis, S. D. M. White, and C. S. Frenk, 1985, 

Astrophys. J. Supp. 57, 241. 
Eisenstein, D. J., 2005, New Astronomy Review 49, 360. 
Eke, V. R., C. M. Baugh, S. Cole, C. S. Frenk, and J. F. Navarro, 

2006, Mon. Not. Royal Astron. Soc. 370, 1 147. 
Fan, X., C. L. Carilli, and B. Keating, 2006, Ann. Rev. Astron. As- 
trophys. 44, 415. 

Feast, M. W., and R. M. Catchpole, 1997, Mon. Not. Royal Astron. 
Soc. 286, LI. 

Feast, M. W., and A. R. Walker, 1987, Ann. Rev. Astron. Astrophys. 
25,345. 

Filippenko, A. V., 1997, Ann. Rev. Astron. Astrophys. 35, 309. 

Finkbeiner, D. P., 2003, Astrophys. J. Supp. 146, 407. 

Fixsen, D. J., E. S. Cheng, J. M. Gales, J. C. Mather, R. A. Shafer, 

and E. L. Wright, 1996, Astrophys. J. 473, 576. 
Fixsen, D. J., and J. C. Mather, 2002, Astrophys. J. 581, 817. 
Fowler, W A., 1989, New York Academy Sciences Annals 571, 68. 
Freedman, W. L., S. M. Hughes, B. F. Madore, J. R. Mould, M. G. 

Lee, P. Stetson, R. C. Kennicutt, A. Turner, L. Ferrarese, H. Ford, 



J. A. Graham, R. Hill, et al, 1994, Astrophys. J. 427, 628. 
Freedman, W. L., B. F. Madore, B. K. Gibson, L. Ferrarese, D. D. 
Kelson, S. Sakai, J. R. Mould, J. R. C. Kennicutt, H. C. Ford, J. A. 
Graham, J. P. Huchra, S. M. G. Hughes, et al, 2001, Astrophys. 
J. 553, 47. 

Fu, L., E. Semboloni, H. Hoekstra, M. Kilbinger, L. van Waerbeke, 
I. Tereno, Y. Mellier, C. Heymans, J. Coupon, K. Benabed, J. Ben- 
jamin, E. Bertin, et al, 2008, Astron. Astrophys. 479, 9. 

Gamow, G., 1948, Physical Review 74, 505. 

Gaztanaga, E., M. Manera, and T. Multamaki, 2006, Mon. Not. 
Royal Astron. Soc. 365, 171. 

Giannantonio, T, R. G. Crittenden, R. C. Nichol, R. Scranton, G. T. 
Richards, A. D. Myers, R. J. Brunner, A. G. Gray, A. J. Connolly, 
and D. P. Schneider, 2006, Phys. Rev. D 74(6), 063520. 

Gould, A., 1995, Astrophys. J. 452, 189. 

Grundahl, F, D. A. VandenBerg, R. A. Bell, M. I. Andersen, and 

P. B. Stetson, 2000, Astron. J. 120, 1884. 
Gunn, J. E., and J. R. I. Gott, 1972, Astrophys. J. 176, 1. 
Gunn, J. E., and B. A. Peterson, 1965, Astrophys. J. 142, 1633. 
Guth, A. H., and S.-Y. Pi, 1982, Physical Review Letters 49, 1 1 10. 
Halliwell, J. J., and S. W. Hawking, 1985, Phys. Rev. D 31, 1777. 
Hamilton, A. J. S., 1993a, Astrophys. J. Lett. 406, L47. 
Hamilton, A. J. S., 1993b, Astrophys. J. 417, 19. 
Hansen, B. M. S., and J. Liebert, 2003, Ann. Rev. Astron. Astrophys. 

41,465. 

Hansen, B. M. S., H. B. Richer, G. G. Fahlman, P. B. Stetson, 
J. Brewer, T. Currie, B. K. Gibson, R. Ibata, R. M. Rich, and 
M. M. Shara, 2004, Astrophys. J. Supp. 155, 551. 

Harrison, E. R., 1970, Phys. Rev. D 1, 2726. 

Haslam, C. G. T, C. J. Salter, H. Stoffel, and W. E. Wilson, 1982, 
Astron. Astrophys. Supp. 47, 1. 

Hawkins, E., S. Maddox, S. Cole, O. Lahav, D. S. Madgwick, 
P. Norberg, J. A. Peacock, I. K. Baldry, C. M. Baugh, J. Bland- 
Hawthorn, T. Bridges, R. Cannon, et al, 2003, Mon. Not. Royal 
Astron. Soc. 346, 78. 

Herschel, W., 1786, Philosophical Transactions Series I 76, 457. 

Herschel, W., 1789, Philosophical Transactions Series 1 79, 212. 

Hetterscheidt, M., P. Simon, M. Schirmer, H. Hildebrandt, T. Schrab- 
back, T. Erben, and P. Schneider, 2007, Astron. Astrophys. 468, 
859. 

Heymans, C, M. L. Brown, M. Barden, J. A. R. Caldwell, K. Jahnke, 
C. Y. Peng, H.-W. Rix, A. Taylor, S. V. W. Beckwith, E. F. Bell, 
A. Borch, B. HauBler, et al, 2005, Mon. Not. Royal Astron. Soc. 
361, 160. 

Heymans, C, L. Van Waerbeke, D. Bacon, J. Berge, G. Bernstein, 
E. Bertin, S. Bridle, M. L. Brown, D. Clowe, H. Dahle, T. Erben, 
M. Gray, etal, 2006, Mon. Not. Royal Astron. Soc. 368, 1323. 

Hillebrandt, W., and J. C. Niemeyer, 2000, Ann. Rev. Astron. Astro- 
phys. 38, 191. 

Hockney, R. W., and J. W. Eastwood, 1988, Computer simulation 

using particles (Bristol: Hilger, 1988). 
Hoekstra, H., 2004, Mon. Not. Royal Astron. Soc. 347, 1337. 
Hoekstra, H., 2007, Mon. Not. Royal Astron. Soc. 379, 317. 
Hoekstra, H., Y. Mellier, L. van Waerbeke, E. Semboloni, L. Fu, 

M. J. Hudson, L. C. Parker, I. Tereno, and K. Benabed, 2006, 

Astrophys. J. 647, 116. 
Hoekstra, H., L. van Waerbeke, M. D. Gladders, Y. Mellier, and 

H. K. C. Yee, 2002a, Astrophys. J. 577, 604. 
Hoekstra, H, H. K. C. Yee, and M. D. Gladders, 2002b, Astrophys. 

J. 577, 595. 

Holz, D. E., and R. M. Wald, 1998, Phys. Rev. D 58(6), 063501. 

Hu, W., 2002a, Phys. Rev. D 66(8), 083515. 

Hu, W., 2002b, Phys. Rev. D 65(2), 023003. 

Hu, W., and N. Sugiyama, 1996, Astrophys. J. 471, 542. 



49 



Hubble, E„ and M. L. Humason, 1931, Astrophys. J. 74, 43. 
Hui, X., H. C. Ford, R. Ciardullo, and G. H. Jacoby, 1993, Astrophys. 
J. 414, 463. 

Izotov, Y. I., and T. X. Thuan, 2004, Astrophys. J. 602, 200. 

Jacoby, G. H., 1980, Astrophys. J. Supp. 42, 1. 

Jarvis, M., G. M. Bernstein, P. Fischer, D. Smith, B. Jain, J. A. Tyson, 

and D. Wittman, 2003, Astron. J. 125, 1014. 
Jenkins, A., C. S. Frenk, S. D. M. White, J. M. Colberg, S. Cole, 

A. E. Evrard, H. M. P. Couchman, and N. Yoshida, 2001, Mon. 

Not. Royal Astron. Soc. 321, 372. 
Kaiser, N., 1987, Mon. Not. Royal Astron. Soc. 227, 1. 
Kaiser, N., G. Squires, and T. Broadhurst, 1995, Astrophys. J. 449, 

460. 

Kamionkowski, M., and A. H. Jaffe, 2001, International Journal of 

Modern Physics A 16, 1 16. 
Kasen, D., and S. E. Woosley, 2007, Astrophys. J. 656, 661. 
Kilbinger, M., and P. Schneider, 2004, Astron. Astrophys. 413, 465. 
Kinney, W. H., 1998, Phys. Rev. D 58(12), 123506. 
Klessen, R., and A. Burkert, 1996, Mon. Not. Royal Astron. Soc. 

280, 735. 

Kochanek, C. S., 2002, Astrophys. J. 578, 25. 
Kochanek, C. S., 2003, Astrophys. J. 583, 49. 

Komatsu, E., J. Dunkley, M. R. Nolta, C. L. Bennett, B. Gold, 

G. Hinshaw, N. Jarosik, D. Larson, M. Limon, L. Page, D. N. 
Spergel, M. Halpern, et al, 2008, ArXiv e-prints 803. 

Korn, A. J., F. Grundahl, O. Richard, P. S. Barklem, L. Mashonkina, 
R. Collet, N. Piskunov, and B. Gustafsson, 2006, Nature 442, 657. 

Kovac, J. M., E. M. Leitch, C. Pryke, J. E. Carlstrom, N. W. Halver- 
son, and W. L. Holzapfel, 2002, Nature 420, 772. 

Kregel, M., and P. C. van der Kruit, 2004, Mon. Not. Royal Astron. 
Soc. 352, 787. 

Kunz, M., and D. Sapone, 2007, Physical Review Letters 98(12), 
121301. 

Lahav, O., S. L. Bridle, W. J. Percival, J. A. Peacock, G. Efstathiou, 
C. M. Baugh, J. Bland-Hawthorn, T. Bridges, R. Cannon, S. Cole, 
M. Colless, C. Collins, et al, 2002, Mon. Not. Royal Astron. Soc. 
333, 961. 

Lamarre, J. M., J. L. Puget, F. Bouchet, P. A. R. Ade, A. Benoit, J. P. 
Bernard, J. Bock, P. de Bernardis, J. Charra, F. Couchot, J. De- 
labrouille, G. Efstathiou, et al, 2003, New Astronomy Review 
47, 1017. 

Landy, S. D., and A. S. Szalay, 1993, Astrophys. J. 412, 64. 

Lasserre, T., C. Afonso, J. N. Albert, J. Andersen, R. Ansari, 
E. Aubourg, P. Bareyre, F. Bauer, J. P. Beaulieu, G. Blanc, A. Bou- 
quet, S. Char, et al, 2000, Astron. Astrophys. 355, L39. 

Leibundgut, B., 2001, Ann. Rev. Astron. Astrophys. 39, 67. 

Leibundgut, B., and S. Blondin, 2005, Nuclear Physics B Proceed- 
ings Supplements 138, 10. 

Liddle, A. R., and D. H. Lyth, 2000, Cosmological Inflation and 
Large-Scale Structure (Cosmological Inflation and Large-Scale 
Structure, by Andrew R. Liddle and David H. Lyth, pp. 414. ISBN 
052166022X. Cambridge, UK: Cambridge University Press, April 
2000.). 

Madgwick, D. S., O. Lahav, I. K. Baldry, C. M. Baugh, J. Bland- 
Hawthorn, T. Bridges, R. Cannon, S. Cole, M. Colless, C. Collins, 
W. Couch, G. Dalton, et al, 2002, Mon. Not. Royal Astron. Soc. 
333, 133. 

Madore, B. F, and W. L. Freedman, 1991, Pub. Astron. Soc. Pacific 
103, 933. 

Massey, R., C. Heymans, J. Berge, G. Bernstein, S. Bridle, D. Clowe, 

H. Dahle, R. Ellis, T. Erben, M. Hetterscheidt, F. W. High, C. Hi- 
rata, et al, 2007a, Mon. Not. Royal Astron. Soc. 376, 13. 

Massey, R., J. Rhodes, A. Leauthaud, P. Capak, R. Ellis, A. Koeke- 
moer, A. Refregier, N. Scoville, J. E. Taylor, J. Albert, J. Berge, 



C. Heymans, et al, 2007b, Astrophys. J. Supp. 172, 239. 
Mather, J. C, E. S. Cheng, D. A. Cottingham, J. R. E. Eplee, D. J. 

Fixsen, T. Hewagama, R. B. Isaacman, K. A. Jensen, S. S. Meyer, 
P. D. Noerdlinger, S. M. Read, L. P. Rosen, et al, 1994, Astro- 
phys. J. 420, 439. 
McDonald, P., U. Seljak, R. Cen, D. Shih, D. H. Weinberg, S. Buries, 

D. P. Schneider, D. J. Schlegel, N. A. Bahcall, J. W. Briggs, 
J. Brinkmann, M. Fukugita, et al, 2005, Astrophys. J. 635, 761. 

Mellier, Y, 1999, Ann. Rev. Astron. Astrophys. 37, 127. 

Meyer, B. S., and D. N. Schramm, 1986, Astrophys. J. 311, 406. 

Monaghan, J. J., 2005, Reports of Progress in Physics 68, 1703. 

Mukhanov, V. F, andG. V. Chibisov, 1981, Soviet Journal of Exper- 
imental and Theoretical Physics Letters 33, 532. 

Norberg, P., C. M. Baugh, E. Hawkins, S. Maddox, J. A. Peacock, 
S. Cole, C. S. Frenk, J. Bland-Hawthorn, T. Bridges, R. Cannon, 
M. Colless, C. Collins, et al, 2001, Mon. Not. Royal Astron. Soc. 
328, 64. 

Nusser, A., and M. Haehnelt, 1999, Mon. Not. Royal Astron. Soc. 
303, 179. 

Ostman, L., and E. Mortsell, 2005, Journal of Cosmology and Astro- 
Particle Physics 2, 5. 

Oswalt, T. D., J. A. Smith, M. A. Wood, and P. Hintzen, 1996, Nature 
382, 692. 

Padmanabhan, T, 1993, Structure Formation in the Universe 
(Structure Formation in the Universe, by T. Padmanabhan, 
pp. 499. ISBN 0521424860. Cambridge, UK: Cambridge Univer- 
sity Press, June 1993.). 

Page, L., G. Hinshaw, E. Komatsu, M. R. Nolta, D. N. Spergel, C. L. 
Bennett, C. Barnes, R. Bean, O. Dore, J. Dunkley, M. Halpern, 
R. S. Hill, et al, 2007, Astrophys. J. Supp. 170, 335. 

Pagels, H. R., 1983, Astrophys. Lett. 23, 151. 

Panagia, N., R. Gilmozzi, F. Macchetto, H.-M. Adorf, and R. P. Kir- 

shner, 1991, Astrophys. J. Lett. 380, L23. 
Patterson, C, 1956, Geochim. Cosmochim. Acta 10, 230. 
Peacock, J. A., 1999, Cosmological Physics (Cosmological Physics, 

by John A. Peacock, pp. 704. ISBN 052141072X. Cambridge, 

UK: Cambridge University Press, January 1999.). 
Peacock, J. A., and S. J. Dodds, 1996, Mon. Not. Royal Astron. Soc. 

280, L19. 

Peebles, P. J., and B. Ratra, 2003, Reviews of Modern Physics 75, 
559. 

Peebles, P. J. E., 1973, Astrophys. J. 185, 413. 

Peebles, P. J. E., 1982, Astrophys. J. Lett. 263, LI. 

Peebles, P. J. E., 1993, Principles of physical cosmology (Prince- 
ton Series in Physics, Princeton, NJ: Princeton University Press, 
— cl993). 

Peebles, P. J. E., and B. Ratra, 1988, Astrophys. J. Lett. 325, L17. 
Peebles, P. J. E., and J. T. Yu, 1970, Astrophys. J. 162, 815. 
Peimbert, M., V. Luridiana, and A. Peimbert, 2007, Astrophys. J. 
666, 636. 

Penzias, A. A., and R. W. Wilson, 1965, Astrophys. J. 142, 419. 

Percival, W. J., R. C. Nichol, D. J. Eisenstein, J. A. Frieman, 
M. Fukugita, J. Loveday, A. C. Pope, D. P. Schneider, A. S. Sza- 
lay, M. Tegmark, M. S. Vogeley, D. H. Weinberg, et al, 2007, 
Astrophys. J. 657, 645. 

Perlmutter, S., G. Aldering, G. Goldhaber, R. A. Knop, P. Nugent, 
P. G. Castro, S. Deustua, S. Fabbro, A. Goobar, D. E. Groom, I. M. 
Hook, A. G. Kim, et al, 1999, Astrophys. J. 517, 565. 

Pettini, M., and D. V. Bowen, 2001, Astrophys. J. 560, 41. 

Phillips, M. M., P. Lira, N. B. Suntzeff, R. A. Schommer, M. Hamuy, 
and J. Maza, 1999, Astron. J. 118, 1766. 

Pierpaoli, E., S. Borgani, D. Scott, and M. White, 2003, Mon. Not. 
Royal Astron. Soc. 342, 163. 

Pierpaoli, E., D. Scott, and M. White, 2001, Mon. Not. Royal Astron. 



50 



Soc. 325, 77. 

Pinsonneault, M. H., G. Steigman, T. P. Walker, and V. K. 

Narayanan, 2002, Astrophys. J. 574, 398. 
Press, W. H., and P. Schechter, 1974, Astrophys. J. 187, 425. 
Refregier, A., 2003, Ann. Rev. Astron. Astrophys. 41, 645. 
Refsdal, S., 1964, Mon. Not. Royal Astron. Soc. 128, 307. 
Reiprich, T. H., and H. Bohringer, 2002, Astrophys. J. 567, 716. 
Renzini, A., and F. Fusi Pecci, 1988, Ann. Rev. Astron. Astrophys. 

26, 199. 

Rhodes, J., A. Refregier, N. R. Collins, J. P. Gardner, E. J. Groth, 

and R. S. Hill, 2004, Astrophys. J. 605, 29. 
Riess, A. G., A. V. Filippenko, P. Challis, A. Clocchiatti, A. Diercks, 

P. M. Garnavich, R. L. Gilliland, C. J. Hogan, S. Jha, R. P. Kir- 

shner, B. Leibundgut, M. M. Phillips, et ai, 1998, Astron. J. 116, 

1009. 

Riess, A. G., L.-G. Strolger, J. Tonry, S. Casertano, H. C. Fergu- 
son, B. Mobasher, P. Challis, A. V. Filippenko, S. Jha, W. Li, 
R. Chornock, R. P. Kirshner, et al, 2004, Astrophys. J. 607, 665. 

Rix, H.-W., and G. Lake, 1993, Astrophys. J. Lett. 417, L1+. 

Ropke, F. K., W. Hillebrandt, W. Schmidt, J. C. Niemeyer, S. I. Blin- 
nikov, and P. A. Mazzali, 2007, Astrophys. J. 668, 1 132. 

Rosati, P., S. Borgani, and C. Norman, 2002, Ann. Rev. Astron. As- 
trophys. 40, 539. 

Sachs, R. K., and A. M. Wolfe, 1967, Astrophys. J. 147, 73. 

Salaris, M., S. degl'Innocenti, and A. Weiss, 1997a, Astrophys. J. 
484, 986. 

Salaris, M., I. Dominguez, E. Garcia-Berro, M. Hernanz, J. Isern, 
and R. Mochkovitch, 1997b, Astrophys. J. 486, 413. 

Salaris, M., and A. Weiss, 2001, Astron. Astrophys. 376, 955. 

Salpeter, E. E., 1955, Astrophys. J. 121, 161. 

Sandage, A., and C. Cacciari, 1990, Astrophys. J. 350, 645. 

Sandage, A., and G. A. Tammann, 1974, Astrophys. J. 190, 525. 

Sandage, A., G. A. Tammann, A. Saha, B. Reindl, F. D. Macchetto, 
and N. Panagia, 2006, Astrophys. J. 653, 843. 

Schechter, P., 1976, Astrophys. J. 203, 297. 

Schlegel, D. J., D. P. Finkbeiner, and M. Davis, 1998, Astrophys. J. 
500,525. 

Schmidt, B. P., R. P. Kirshner, and R. G. Eastman, 1992, Astrophys. 
J. 395, 366. 

Schneider, P., J. Ehlers, and E. E. Falco, 1992, Gravitational Lenses 
(Gravitational Lenses, XIV, 560 pp. 112 figs.. Springer- Verlag 
Berlin Heidelberg New York. Also Astronomy and Astrophysics 
Library). 

Scranton, R., B. Menard, G. T. Richards, R. C. Nichol, A. D. Myers, 

B. Jain, A. Gray, M. Bartelmann, R. J. Brunner, A. J. Connolly, 
J. E. Gunn, R. K. Sheth, et al. , 2005, Astrophys. J. 633, 589. 

Seljak, U., 2002, Mon. Not. Royal Astron. Soc. 337, 769. 
Seljak, U., and M. Zaldarriaga, 1996, Astrophys. J. 469, 437. 
Seljak, U., and M. Zaldarriaga, 1997, Physical Review Letters 78, 
2054. 

Semboloni, E., Y. Mellier, L. van Waerbeke, H. Hoekstra, I. Tereno, 
K. Benabed, S. D. J. Gwyn, L. Fu, M. J. Hudson, R. Maoli, and 
L. C. Parker, 2006, Astron. Astrophys. 452, 51. 

Sheth, R. K., and G. Tormen, 2002, Mon. Not. Royal Astron. Soc. 
329,61. 

Silk, J., 1968, Astrophys. J. 151, 459. 

Simon, P., M. Hetterscheidt, M. Schirmer, T. Erben, P. Schneider, 

C. Wolf, and K. Meisenheimer, 2007, Astron. Astrophys. 461, 
861. 

Slipher, V. M., 1927, Pub. Astron. Soc. Pacific 39, 143. 

Smith, G. P., A. C. Edge, V. R. Eke, R. C. Nichol, I. Smail, and J.-P. 

Kneib, 2003a, Astrophys. J. Lett. 590, L79. 
Smith, R. E., J. A. Peacock, A. Jenkins, S. D. M. White, C. S. Frenk, 

F. R. Pearce, P. A. Thomas, G. Efstathiou, and H. M. P. Couchman, 



2003b, Mon. Not. Royal Astron. Soc. 341, 1311. 
Sofue, Y, and V. Rubin, 2001, Ann. Rev. Astron. Astrophys. 39, 137. 
Spergel, D. N, R. Bean, O. Dore, M. R. Nolta, C. L. Bennett, 

J. Dunkley, G. Hinshaw, N. Jarosik, E. Komatsu, L. Page, H. V. 

Peiris, L. Verde, et al, 2007, Astrophys. J. Supp. 170, 377. 
Spite, F, andM. Spite, 1982, Astron. Astrophys. 115, 357. 
Springel, V, and L. Hernquist, 2002, Mon. Not. Royal Astron. Soc. 

333, 649. 

Springel, V, S. D. M. White, A. Jenkins, C. S. Frenk, N. Yoshida, 
L. Gao, J. Navarro, R. Thacker, D. Croton, J. Helly, J. A. Peacock, 
S. Cole, et al. , 2005, Nature 435, 629. 

Steigman, G, 2007, Annual Review of Nuclear and Particle Science 
57, 463. 

Steigman, G., D. N. Schramm, and J. E. Gunn, 1977, Physics Letters 
B 66, 202. 

Steigman, G, and T. P. Walker, 1992, Astrophys. J. Lett. 385, L13. 
Straniero, O., and A. Chieffi, 1991, Astrophys. J. Supp. 76, 525. 
Sunyaev, R. A., and Y. B. Zeldovich, 1970, Astrophys. Space Sc. 7, 
3. 

Sunyaev, R. A., and Y. B. Zeldovich, 1972, Comments on Astro- 
physics and Space Physics 4, 173. 

Tauber, J. A., 2004, Advances in Space Research 34, 491. 

Tegmark, M., M. R. Blanton, M. A. Strauss, F. Hoyle, D. Schlegel, 
R. Scoccimarro, M. S. Vogeley, D. H. Weinberg, I. Zehavi, 
A. Berlind, T. Budavari, A. Connolly, et al, 2004, Astrophys. J. 
606,702. 

Toffolatti, L., F. Argueso Gomez, G. de Zotti, P. Mazzei, A. Frances- 

chini, L. Danese, and C. Burigana, 1998, Mon. Not. Royal Astron. 

Soc. 297, 117. 
Tonry, J., and D. P. Schneider, 1988, Astron. J. 96, 807. 
Tonry, J. L., J. P. Blakeslee, E. A. Ajhar, and A. Dressier, 1997, 

Astrophys. J. 475, 399. 
Tonry, J. L., J. P. Blakeslee, E. A. Ajhar, and A. Dressier, 2000, 

Astrophys. J. 530, 625. 
Truran, J. W., W. D. Arnett, and A. G. W. Cameron, 1967, Canadian 

Journal of Physics 45, 2315. 
Tully, R. B., and J. R. Fisher, 1977, Astron. Astrophys. 54, 661. 
Udalski, A., M. Szymanski, M. Kubiak, G. Pietrzynski, I. Soszynski, 

P. Wozniak, and K. Zebrun, 1999, Acta Astronomica 49, 201. 
Uson, J. M., and D. T. Wilkinson, 1984, Astrophys. J. Lett. 277, LI. 
Valenziano, L., M. Sandri, G. Morgante, C. Burigana, M. Bersanelli, 

R. C. Butler, F. Cuttaia, F. Finelli, E. Franceschi, M. Galaverni, 

A. Gruppuso, M. Malaspina, et al, 2007, New Astronomy Review 

51, 287. 

Van Waerbeke, L., Y. Mellier, and H. Hoekstra, 2005, Astron. Astro- 
phys. 429, 75. 

Vandenberg, D. A., 1983, Astrophys. J. Supp. 51, 29. 

VandenBerg, D. A., 2000, Astrophys. J. Supp. 129, 315. 

Viel, M., G. D. Becker, J. S. Bolton, M. G. Haehnelt, M. Rauch, and 
W. L. W. Sargent, 2008, Physical Review Letters 100(4), 041304. 

Viel, M., and M. G. Haehnelt, 2006, Mon. Not. Royal Astron. Soc. 
365,231. 

Viel, M., J. Weller, and M. G. Haehnelt, 2004, Mon. Not. Royal 

Astron. Soc. 355, L23. 
Wambsganss, J., and B. Paczynski, 1992, Astrophys. J. Lett. 397, LI. 
Wang, X., M. Tegmark, B. Jain, and M. Zaldarriaga, 2003, Phys. 

Rev. D 68(12), 123001. 
Wasserburg, G. J., F. Tera, D. A. Papanastassiou, and J. C. Huneke, 

1977, Earth and Planetary Science Letters 35, 294. 
Weinberg, S., 2008, Cosmology (Oxford University Press). 
Wesselink, A. J., 1946, Bull. Astron. Inst. Netherlands 10, 91. 
Wetterich, C, 1988, Nuclear Physics B 302, 668. 
White, S. D. M., C. S. Frenk, M. Davis, and G. Efstathiou, 1987, 

Astrophys. J. 313, 505. 



51 

White, S. D. M., J. F. Navarro, A. E. Evrard, and C. S. Frenk, 1993, 
Nature 366, 429. 

Woosley, S. E., A. Heger, and T. A. Weaver, 2002, Reviews of Mod- 
ern Physics 74, 1015. 

Yang, J., M. S. Turner, D. N. Schramm, G. Steigman, and K. A. 
Olive, 1984, Astrophys. J. 281, 493. 

York, D. G., J. Adelman, J. J. E. Anderson, S. F. Anderson, J. An- 
nis, N. A. Bahcall, J. A. Bakken, R. Barkhouser, S. Bastian, 
E. Berman, W. N. Boroski, S. Bracker, et al, 2000, Astron. J. 
120, 1579. 

Zaldarriaga, M„ 1997, Phys. Rev. D 55, 1822. 

Zaroubi, S., M. Viel, A. Nusser, M. Haehnelt, and T.-S. Kim, 2006, 
Mon. Not. Royal Astron. Soc. 369, 734. 

Zel'dovich, Y. B., 1970, Astron. Astrophys. 5, 84. 

Zeldovich, Y. B., 1972, Mon. Not. Royal Astron. Soc. 160, IP 

Zlatev, I., L. Wang, and P. J. Steinhardt, 1999, Physical Review Let- 
ters 82, 896. 

Zwicky, F, 1933, Helvetica Physica Acta 6, 110. 
Zwicky, F, 1937, Astrophys. J. 86, 217. 



